Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,031,324	A→C	Q55P (CAA→CCA)	pepQ →	proline dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,031,324	0	A	C	100.0%	217.5 / NA	59	Q55P (CAA→CCA)	pepQ	proline dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (35/24);  total (35/24)

ATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAG‑CTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCC  >  NC_000913/4031093‑4031552                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     accgcGCTATGAATGACAGATGTGATAATAGGCGGATACTTAAGGTCAAAAAATCAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAACCATATAG‑CTACCTTACAAAAACGGACCCGCGATGGCGGGGCGCGCTTCAAGCTGGAGGCGTTCCTTATCCCCCCCGGCGAGCTGTTCAAGGTTTTCCTCGAGGATCACCCCAACCCGTTTAAA‑GGAAA‑CCCGCCAtt                                                                                                                                                                                                                                   <  2:442490/231‑1 (MQ=255)          tGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAG‑CTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTg                                                                                                                                                                                                                          >  2:567320/1‑235 (MQ=255)                    aTGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAG‑CTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCCAAGCGTAGCTGCCGGTaa                                                                                                                                                                                                               >  2:122404/1‑236 (MQ=255)                              aTGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAG‑CTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGcc                                                                                                                                                                                                     <  1:567320/236‑1 (MQ=255)                                 cGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAG‑CTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGc                                                                                                                                                                                                      >  1:743909/1‑232 (MQ=255)                                   gagaCTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAG‑CTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCaa                                                                                                                                                                                                   >  1:448948/1‑233 (MQ=255)                                          aaGGTCAAAAAAACAGAAGGATAAAACATGAAATCACGGCCCTCGCTCTATAAAAATCATATAG‑CTACTTTACAAGAACGGACTCGCGATGCGCTGGCGCGTTTCAAGCTGGATGCTTTACTTATCCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAC‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACtgctggt                                                                                                                                                                                          <  2:50112/235‑1 (MQ=255)                                                aaaaaaaCAGAAGGGTAAAACATGGAATCACTGGCCTCGCTCTATAAAAATCATATAG‑CTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACtgctggttgctggt                                                                                                                                                                                   >  1:34288/1‑236 (MQ=255)                                                                aaaaatTGAAATCACTGGCCTCTCTATAAAAAACTCATATAC‑CTCCCGAACAAGAACTGCCCCCCTATGCGCTGGCGCGCTTCATGCTGATTGCTTTCCTTATCCACTCGGGCGAGCTGTTCAATTTTTTCCTCGACGATCACCCCTACCCGTTTAAA‑GTAAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGaa                                                                                                                                                                       <  2:67631/232‑1 (MQ=255)                                                                       ggAATCACTGGCCTCGCTCTATAAAAATCATATAG‑CTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGcc                                                                                                                                                                 >  1:244751/1‑231 (MQ=255)                                                                                 cccTCACTCTA‑AAAAATCTAATAGCCTACCTTACAAGAGCGGCCTCAGGAGGCGCGGGCGGTCTTAAAGCTGGATGCTTTATTTTTTCCCCCCGGCGAGCTGTTCATTGTTTTCCTGGACGATCACCCCTATCCGTTTAAA‑GTGAC‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTa                                                                                                                                                  <  2:522964/235‑1 (MQ=255)                                                                                  cctcgtctaaaaaaacTCATTTGC‑CTACCTTACAAGAACGGACCCGGGGTGCCGTGGCGCGCTTCAAGCTGGATGCTTTCCTTATTCCCCCCGGCGAGTTGTTCAATTTTTTTCCCGCCGATCATCCCTATCCGTTAAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTAt                                                                                                                                                 <  2:34288/225‑1 (MQ=255)                                                                                                 aTCATATAG‑CTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTg                                                                                                                                       >  1:738687/1‑231 (MQ=255)                                                                                                      ataG‑CTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGCTTACTGGc                                                                                                                              >  1:146190/1‑235 (MQ=255)                                                                                                             aCCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCtggtgg                                                                                                                                                                                 <  1:617250/178‑1 (MQ=255)                                                                                                             aCCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCtggtgg                                                                                                                                                                                 >  2:617250/1‑178 (MQ=255)                                                                                                             aCCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGt                                                                                                                       >  2:597618/1‑236 (MQ=255)                                                                                                              ccTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGt                                                                                                                       >  1:422506/1‑235 (MQ=255)                                                                                                               cTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGt                                                                                                                       <  1:597618/234‑1 (MQ=255)                                                                                                                 tACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGa                                                                                                                    >  1:512917/1‑235 (MQ=255)                                                                                                                 tACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGa                                                                                                                    >  1:574937/1‑235 (MQ=255)                                                                                                                   cAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGa                                                                                                                    >  1:457982/1‑233 (MQ=255)                                                                                                                             ctgcgagggcgcTGGCTGGTTTCATGCGGGAGGCTTTACATTTCCACTCGGGCGAGCTTTAAATTGTTTTCCTCGACATTCACCCCTACCCGTTAAAA‑GTAAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGcc                                                                                                           <  2:146190/224‑1 (MQ=255)                                                                                                                                   aTGCGGTGGCGCGCTTAAGGTTAGAGGCTTTACTTTTTCACACCGGGGAGCTGTTCAATTTTTTCCGCGACGACCCCCCCTACCCGTTTAAG‑GAGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTc                                                                                                 <  2:512917/236‑1 (MQ=255)                                                                                                                                                cTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCt                                                                                                >  1:591808/1‑224 (MQ=255)                                                                                                                                                    aagtgggATGCGTTTCTTTTTCACCCGGGCGGGCTTTTCAATGTTTTTCTGGACGATCACCCCTACCCGTTTAAA‑GGAAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGcaca                                                                                                                           <  2:435687/188‑1 (MQ=255)                                                                                                                                                          gATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACtgctggt                                                                                                                                                                                          >  2:389073/1‑124 (MQ=255)                                                                                                                                                          gATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACtgctggt                                                                                                                                                                                          <  1:389073/124‑1 (MQ=255)                                                                                                                                                           aTGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGa                                                                                                                    >  1:435687/1‑193 (MQ=255)                                                                                                                                                           aTGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCAAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGAt                                                                          >  1:81653/1‑235 (MQ=255)                                                                                                                                                                taccTTATTCACTCCGGCGAGTTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCgcgc                                                                     <  1:715684/232‑1 (MQ=255)                                                                                                                                                                             ccGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCCCTTTCTGGACTGAAGATGCGGAAGGGATGCGGCTGCCGTAAGCCg                                                         >  1:679788/1‑234 (MQ=255)                                                                                                                                                                              cGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATg                                                      >  2:758425/1‑236 (MQ=255)                                                                                                                                                                               ggggAGCTGTGCAGTGTTTTTCGCGACGATCATCCCTACCCGTTTAAA‑GTGAC‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATgg                                                     <  1:653099/233‑1 (MQ=255)                                                                                                                                                                               ggCGAGCTTTTAAGTGTTTTTCTCGACGTTCATCCATATCCGTTTAAA‑GTGAC‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATgg                                                     <  1:242094/236‑1 (MQ=255)                                                                                                                                                                                          cAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcc                                          >  2:220315/1‑236 (MQ=255)                                                                                                                                                                                                  ttCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGa                                                                                         >  2:101762/1‑181 (MQ=255)                                                                                                                                                                                                  ttCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGa                                                                                         <  1:101762/181‑1 (MQ=255)                                                                                                                                                                                                  ttCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctg                                     >  2:127154/1‑233 (MQ=255)                                                                                                                                                                                                  ttCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctg                                     >  1:389356/1‑233 (MQ=255)                                                                                                                                                                                                       gacgaTCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctg                                        >  2:462493/1‑225 (MQ=255)                                                                                                                                                                                                       gacgaTCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCctgctgcctgctgc                             >  2:121965/1‑236 (MQ=255)                                                                                                                                                                                                         cgaTCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGCAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTgcgc                           >  1:660323/1‑236 (MQ=255)                                                                                                                                                                                                         caatcatcCTTATCCGTTAAAA‑GTGAT‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTgcgc                           <  1:61517/234‑1 (MQ=255)                                                                                                                                                                                                          gaTCATCCCTATCCGTTAAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTgcgcg                          <  1:518571/236‑1 (MQ=255)                                                                                                                                                                                                              atcCCTATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGCCCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGCAAGCCGCTGGACTTGGTAGCGTGCTGCTTGCTgcgcgagcg                      >  2:527884/1‑232 (MQ=255)                                                                                                                                                                                                                  ttATCCGTTAAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAata                  <  1:471555/235‑1 (MQ=255)                                                                                                                                                                                                                   tATCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCTGGCTGCCTGATGCGGCCGGCa                     >  2:471555/1‑232 (MQ=255)                                                                                                                                                                                                                   cgtCCCTTTAAA‑GACAACCCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCAAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAata                  <  2:81653/234‑1 (MQ=255)                                                                                                                                                                                                                     tCCGTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGACCAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGCCGAACCGCCTCCGACCTCTTTCTGGACTGAAGATGTGTAAAGGATCGCGCTGCCGAAAGCCGATGGCATTGGTGACCTCCTACCCGCTGCCCTCGGCAAAATCg               >  2:567495/1‑236 (MQ=255)                                                                                                                                                                                                                       cctttAAAGGTGAC‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCgg              <  1:474829/234‑1 (MQ=255)                                                                                                                                                                                                                        gTTTAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGtt            >  1:407034/1‑236 (MQ=255)                                                                                                                                                                                                                          ttAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGTTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGTTGCCTGCTCCGTGCGGCAATATCCGTtaa          >  2:692054/1‑235 (MQ=255)                                                                                                                                                                                                                          ttAAA‑GTGAA‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGGCCGGCAATATCCGttca          >  2:8195/1‑234 (MQ=255)                                                                                                                                                                                                                             aa‑GTGAA‑TCCGCCATTCAAAGCGTGGGTGCCGGTACCTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAGGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTatat        <  1:567495/235‑1 (MQ=255)                                                                                                                                                                                                                                  aaa‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTc   <  1:122404/235‑1 (MQ=255)                                                                                                                                                                                                                                    a‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGa                                                        >  2:750841/1‑181 (MQ=255)                                                                                                                                                                                                                                    a‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGa                                                        <  1:750841/181‑1 (MQ=255)                                                                                                                                                                                                                                    a‑CCCGCCATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGTGCGCGGCAATATCGGTTTAATCGGTcc  >  2:176223/1‑235 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     ATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAG‑CTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAA‑GTGAA‑CCCGCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCC  >  NC_000913/4031093‑4031552

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑