Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	715,762	G→A	7.7%	K117K (AAG→AAA)	ybfP →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	715,762	0	G	A	7.7%	218.9 / 7.4	78	K117K (AAG→AAA)	ybfP	lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (42/30);  new base A (3/3);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

ACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATA  >  NC_000913/715615‑715910                                                                                                                                                    |                                                                                                                                                     aCCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGTGCAAGc                                                                                                                                                     >  2:1124722/1‑149 (MQ=255)   cTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGccc                                                                                                                                                   >  1:327555/1‑149 (MQ=255)    tGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAACAGAACTTACTTATCTTCATATTATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGCTTGCCCTGGCAAACCCg                                                                                                                                                  <  1:63150/149‑1 (MQ=255)    tGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCg                                                                                                                                                  >  2:733781/1‑149 (MQ=255)      ataataCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGCTTGCCCTGGCAAACCCgca                                                                                                                                                <  2:323941/149‑3 (MQ=255)           aCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCg                                                                                                                                           >  1:689931/1‑149 (MQ=255)            cccGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGa                                                                                                                                          >  1:1131040/1‑149 (MQ=255)            cccGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGa                                                                                                                                          <  2:1000511/149‑1 (MQ=255)             ccgccgCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCCCAAGAAAAGGTCGGTTTAAATCCCGGATGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCTGCTCTCGAc                                                                                                                                         >  1:767041/1‑149 (MQ=255)             ccgccgCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGAc                                                                                                                                         >  2:4678/1‑149 (MQ=255)                      gCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCa                                                                                                                                >  1:322099/1‑149 (MQ=255)                        aGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAATAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTGTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAAt                                                                                                                              <  2:737019/149‑1 (MQ=255)                        aGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTAAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGCGTTTGCCAGGGCAAGCCCGGCAATCGACCGAAAGTCaaa                                                                                                                              >  1:1057916/1‑148 (MQ=255)                         gTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGGTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCCAGCACGGCACTCGACCGATAGTCAAtt                                                                                                                             >  1:481513/1‑148 (MQ=255)                            gTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAgg                                                                                                                          >  2:137625/1‑149 (MQ=255)                               aGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGCCCGATAGTCAATGAGGAtt                                                                                                                       >  2:151248/1‑149 (MQ=255)                                 gCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTAt                                                                                                                     >  2:45950/1‑149 (MQ=255)                                  cGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGGCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATg                                                                                                                    <  1:527168/149‑1 (MQ=255)                                         gAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTAAACCTCATTTTTTATGCGGGTTTGCCATTGCAAGCACGGAACTCGACCGATATTCAATGAGGATTATGGCAGCGa                                                                                                             >  2:840288/1‑149 (MQ=255)                                           aaatGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGaaa                                                                                                           <  1:817085/149‑1 (MQ=255)                                           aaaaGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGaaa                                                                                                           <  1:1055267/149‑1 (MQ=255)                                            aaaGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAg                                                                                                          >  2:243046/1‑149 (MQ=255)                                             aaGAACTTACTTATCATCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGc                                                                                                         <  1:21475/149‑1 (MQ=255)                                                       ttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGTTTTGCCAGGGCAAGCCCGGCCCTCGACCGATAGTCAATGAGGATTATGGAAGCGAAAGCGGAATTTTTg                                                                                               >  2:58676/1‑149 (MQ=255)                                                       ttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCAGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTg                                                                                               >  2:739823/1‑149 (MQ=255)                                                             tCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCga                                                                                         <  1:1170364/149‑1 (MQ=255)                                                             tCATAATATGTCAGCACAAGAAAAGGGCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCga                                                                                         <  2:462297/149‑1 (MQ=255)                                                               ataataTGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTAAGCGAAAGCGGATTTTTTGGCTCgagt                                                                                       >  1:573629/1‑147 (MQ=255)                                                                taataTGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgatga                                                                                      >  1:200721/1‑149 (MQ=255)                                                                 aataTGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgatgat                                                                                     <  1:733781/149‑1 (MQ=255)                                                                  atatGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTTGCTCGatgatg                                                                                    >  2:1098103/1‑149 (MQ=255)                                                                     tGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCAt                                                                                 >  1:411821/1‑149 (MQ=255)                                                                       tCAGCTCAAGAAAAGGTCGGTTTAAATCCAGGTTGGACATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGTCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGt                                                                               <  1:207893/149‑1 (MQ=255)                                                                          gCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCgg                                                                            >  2:1015340/1‑149 (MQ=255)                                                                             caAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGat                                                                         <  1:341577/149‑1 (MQ=255)                                                                             caAGAAAACGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTATGCCAGCGCAAGCCAGGCACTCGACCGAAAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGAAGATGCATGTCGGGat                                                                         >  1:352477/1‑149 (MQ=255)                                                                                gAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGattat                                                                      >  1:1122896/1‑149 (MQ=255)                                                                                gAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGAAGCATGTCGGGattat                                                                      >  1:348554/1‑149 (MQ=255)                                                                                  aaaGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTg                                                                    <  1:603119/149‑1 (MQ=255)                                                                                   aaGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGa                                                                   >  1:532620/1‑149 (MQ=255)                                                                                              aaaTCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCt                                                        >  2:526975/1‑149 (MQ=255)                                                                                                  cccGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCa                                                    >  2:77356/1‑149 (MQ=255)                                                                                                    cGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGt                                                  <  1:1124722/149‑1 (MQ=255)                                                                                                    cGGTTGGCAATGCTACACCTCATTTTATATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGt                                                  <  2:406789/149‑1 (MQ=255)                                                                                                      gTTGGCAATGCTACACCTCATTTTTTATTCGGGTATGCCAGGGAAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTTTTGAGTTACGACGCTGCCAGTCa                                                <  1:694336/149‑1 (MQ=255)                                                                                                       ttGGCAATGCTACACCTCATTTTTTATGCGGGTTTGACATGGCAAGCCCGGGACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAg                                               <  1:155191/149‑1 (MQ=255)                                                                                                            aaTGCTACACCTCATTTTTTATGCGGGTTTGACATTGCAAGCCCGGCACTCGACAGATAGACAATGATGATTATGTAAGCTAAAGCGGATTTTTTGGATCGATGATGCATGTCGGGAATATTGAGTTACGACGCTGCCAGTCAGAAAAt                                          >  2:697981/1‑149 (MQ=255)                                                                                                             aTGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAAtt                                         <  1:85508/149‑1 (MQ=255)                                                                                                               gCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAACGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGc                                       >  2:1058852/1‑149 (MQ=255)                                                                                                                 tACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCa                                     <  2:1131040/149‑1 (MQ=255)                                                                                                                 tACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCa                                     >  2:302053/1‑149 (MQ=255)                                                                                                                     ccTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCaa                                 <  1:137625/149‑1 (MQ=255)                                                                                                                       tCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGCCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGa                               >  1:1008889/1‑149 (MQ=255)                                                                                                                       tCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGa                               <  2:689931/149‑1 (MQ=255)                                                                                                                       tCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGa                               >  2:980782/1‑149 (MQ=255)                                                                                                                       tCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGa                               >  1:491935/1‑149 (MQ=255)                                                                                                                       tCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAAGGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGGCGGGAGTTTTGGGTTACGCCGCGGCCAGTCAGCAAATTGccggaatca                               >  2:217884/1‑142 (MQ=255)                                                                                                                        cATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCTATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGaa                              <  2:926595/149‑1 (MQ=255)                                                                                                                        cATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGaa                              <  1:77356/149‑1 (MQ=255)                                                                                                                           tttttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTg                           <  1:45950/149‑1 (MQ=255)                                                                                                                             tttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaa                         <  2:322099/149‑1 (MQ=255)                                                                                                                              ttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaaa                        >  1:1005352/1‑149 (MQ=255)                                                                                                                              ttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaaa                        <  2:532620/149‑1 (MQ=255)                                                                                                                                tgccGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGCATATTTTGGCACGATGATGCATGTCGGGATTATTGAGTTACGACGCTGACAGTCAGAAAAATGCCAGCAAGAATTGAAAGc                      >  2:490918/4‑149 (MQ=255)                                                                                                                                 tGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCg                     <  1:753870/149‑1 (MQ=255)                                                                                                                                  ccGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGa                    >  1:481823/2‑149 (MQ=255)                                                                                                                                   cGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAAtt                            <  2:940374/140‑1 (MQ=255)                                                                                                                                   cGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAAtt                            >  1:940374/1‑140 (MQ=255)                                                                                                                                    gggTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATa                  >  1:271669/1‑149 (MQ=255)                                                                                                                                          gCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATaca            >  1:675617/1‑149 (MQ=255)                                                                                                                                           ccAGGGCAAGCCCGGGACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTAGGGATTATTGAGTTACGACGCTGCGAGTCAGAAATTTGCCAGCAAGAATTGAAAGCGATAAATacac           <  1:227889/149‑1 (MQ=255)                                                                                                                                           ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATacac           <  2:200721/149‑1 (MQ=255)                                                                                                                                           ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATacac           >  1:1144214/1‑149 (MQ=255)                                                                                                                                           ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATacac           >  1:1175151/1‑149 (MQ=255)                                                                                                                                             cgggCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACAt         <  2:78415/148‑1 (MQ=255)                                                                                                                                               ggCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACAtaa       >  1:29506/1‑147 (MQ=255)                                                                                                                                               ggCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTCTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTa       >  1:105339/1‑149 (MQ=255)                                                                                                                                                 cAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTaaa     <  1:243046/149‑1 (MQ=255)                                                                                                                                                    gCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGCCGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAAACACATTAAAATa  >  2:73542/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATA  >  NC_000913/715615‑715910

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑