Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,396,772 T→C F236F (TTT→TTC)  amiB → N‑acetylmuramoyl‑L‑alanine amidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,396,7720TC96.5% 171.2 / ‑3.8 58F236F (TTT→TTC) amiBN‑acetylmuramoyl‑L‑alanine amidase B
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (25/30);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGC  >  NC_000913/4396525‑4397014
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
gCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCaaa                                                                                                                                                                                                                                                 <  1:107248/251‑1 (MQ=255)
        gcgcgcAATCCGTTTAAAACGGAAAGTAACCGCTCTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACTCGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTTCGCGGGAG‑AAAAATGTCACCATC‑GCCATCGTGCGTAAATTGCGTACTTTTCTCAATGACGATCC‑GATGTTCAAAGGCGtttt                                                                                                                                                                                                                                         <  1:348364/251‑1 (MQ=255)
         cgcgcAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGGTGCCGGACACGGCGGGCAGGACCCCGGCGCTATCGGCCCCGGGGGGAAGCGGGGGGAAAAAGGT‑ACCATCCGCCA‑AGCGCGGAAAATGCGACCGTGGCTCAATTACGGTCC‑GATGTTTAACTGCGttttt                                                                                                                                                                                                                                        >  2:442496/1‑250 (MQ=255)
              aaTCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAAcccc                                                                                                                                                                                                                                   >  1:355173/1‑250 (MQ=255)
              aaTCCGTTTAAAACGGAAAGGAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGAGGTCAGGGACCTGGTGCTAGCGGCCCAGGTGGAACGCGGGAG‑GACAATGTCCCCATC‑GCCATAGCGCGTAAATAGGCTACTTTGCTCAATGAAGGTAC‑GATGTTCAAAAGCATTTTacac                                                                                                                                                                                                                                     >  2:160080/1‑246 (MQ=255)
                                     ccGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCCGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAAAGTCACCCTC‑GCCCTCGCGCGTCAATTTCGGACTTTGCTCAATGACGATCC‑GATGTTCAAAAGCGGTTTAACCCGGGGCGGGGATTACTTTATCTCgg                                                                                                                                                                                                            >  2:37912/1‑251 (MQ=255)
                                                 tgtTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCCGGGGGACGCGGGGG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGCTACGGGGATTACTTTATCTCGGTGGTggggggcc                                                                                                                                                                                                >  1:241371/1‑247 (MQ=255)
                                                 ggtTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCa                                                                                                                                                                                                <  1:171808/250‑1 (MQ=255)
                                                    tATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTTCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGGTTTCTTTATCTCGGGGATGGGGCGCAGTg                                                                                                                                                                                             >  1:99802/1‑251 (MQ=255)
                                                      tCAGAAGTAATACGGTAACGCGTCCGGCATCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCCGGACCCTTGCGCTATAGGCCCCCGGGGTACGCGGGGG‑AAAAATGTCACCCTC‑GCCATCGCGCGGAAATTGCGTACTTTTCTCAATGAAGGTCCCGTTGTC‑AAAGGCGTTTTATCCCGTGACGGGTGTTACTTTTTCTCGGGGATGGGgcgcggggat                                                                                                                                                                                           >  2:141222/1‑245 (MQ=255)
                                                        agcagTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGGG‑AAAAATGTAACCATC‑GACATTGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGGTCAAAGGAGTTTTAAACCGTGCCGGGGGATACTTTATCTCGGGGATGGGGCGCAGTggtgg                                                                                                                                                                                         >  1:436722/1‑247 (MQ=255)
                                                              aaTACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACg                                                                                                                                                                                   <  1:171891/251‑1 (MQ=255)
                                                                       aatggtccggccggccgcgtgACGGCTAAAACTGGCGATAAAATGATCATCGCTATTGATGCCGTACACGGCGGTCAGGACCCTGGCGTTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCGCCAGT‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAAGGACGATCC‑GAGGTTCATAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCaaaa                                                                                                                                                                          <  2:444082/236‑1 (MQ=255)
                                                                           gTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAAAGTCCCCATC‑GCCATCGCGCGTAAATTTCGGACTTTGCTCAATGACGATCC‑GGTGTTCCAAGGCGTTTTAACCCGGTACGGGGATTACTTTATCTCGGTGAAGGGGCGCAGTGATGTgggacctcaggaaaaaccc                                                                                                                                                                      >  2:143346/1‑234 (MQ=255)
                                                                                gCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAAttt                                                                                                                                                                 <  2:236155/251‑1 (MQ=255)
                                                                                 cAGCGCGCGCGACGGCGAACCCTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATTGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTc                                                                                                                                                                <  1:75757/251‑1 (MQ=255)
                                                                                        cgcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCCAAGGCGTTTTAACCCGTGACGGGGATTACTTTTTCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCCATTTCCTCgtgt                                                                                                                                                         >  1:295753/1‑251 (MQ=255)
                                                                                        cgcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCgtgt                                                                                                                                                         <  1:37912/251‑1 (MQ=255)
                                                                                          cgACGGCCAAAACTGGCTATAAAATTATCATCGCTAGTGATGCCGGAGGCGGCGGTCAGGACCCTTGCTCTATCGGCCCCAGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGAGGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCg                                                                                                                                                       <  2:389079/251‑1 (MQ=255)
                                                                                           gACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCgg                                                                                                                                                      >  1:58527/1‑250 (MQ=255)
                                                                                            aCGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGGCCGGACCCTGGCGCTATCGGCCCCCGGGGTACGCGGGGG‑AAAAATTTCACCATC‑GCCCTCGCGCGGAAATTGCGGACTTTGCACAATGACGATAC‑GATGTTCAAAGGCGTTTTAACCCGTGCAGGGGATTACTTTATCTCGGGGATGGGGGGACGTGCGGTGGCACGTAAGCAAAACGCCAATTTCCTCGGGTCGAt                                                                                                                                                     >  2:130314/1‑251 (MQ=255)
                                                                                             aggatAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGAACCTGGCGCTATCGGCCCCGGCGGTACGCGGGGG‑AAAAATGTCAACATC‑GCCATCGCGCGTAAATTGCGGAATTTGCTCAATGACGAGCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGACAAACGCCAATTTCCTCGTGTCGAtt                                                                                                                                                    >  1:259813/5‑251 (MQ=255)
                                                                                              ggCTAACACGGGCGATAAAATTATCATCGCGAGTGATGCCGGACACGGCGGTCAGGACCCGGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTGTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTc                                                                                                                                                   <  2:152775/251‑1 (MQ=255)
                                                                                                      cTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGGG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTa                                                                                                                                                                                                                  >  1:412241/1‑180 (MQ=255)
                                                                                                      cTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTa                                                                                                                                                                                                                  <  2:412241/180‑1 (MQ=255)
                                                                                                               aaaTTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCg                                                                                                                                  <  1:250879/251‑1 (MQ=255)
                                                                                                                     tcatcGCTATTGATGACGGACCCGGCGGTCAGGACCCGGGCGCTATCGGCCCCGGGGGTACGCGGGAG‑AAAAATGGCACCATC‑GCCAGCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCg                                                                                                                             <  2:250388/250‑1 (MQ=255)
                                                                                                                            tATTGATGCCGGACACGGCGGTCAGGGCCCTGGCGCCATTGGCCCCGGTGGGACGAGGGAG‑AAAAATGTCAACATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAAACCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGGGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCCGTGCg                                                                                                                      >  1:172730/1‑250 (MQ=255)
                                                                                                                               tGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGGGTCGCTTCCCGCTGATGCCGCACCGAACCGCAGGGGGACTg                                                                                                                  >  1:18644/1‑251 (MQ=255)
                                                                                                                                 aTGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGgc                                                                                                                <  2:145395/251‑1 (MQ=255)
                                                                                                                                 aTGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGgc                                                                                                                <  1:207771/251‑1 (MQ=255)
                                                                                                                                   gCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGgcg                                                                                                               <  1:48408/250‑1 (MQ=255)
                                                                                                                                     cGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCtt                                                                                                            <  1:54357/251‑1 (MQ=255)
                                                                                                                                            ggcggTCCGGACCGTGGCGCTTGCGGCCCCGGGGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATg                                                                                                     <  2:435199/251‑1 (MQ=255)
                                                                                                                                                     gACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTACCCCGTGACGGGGATTACTTTATCTCGGTGATGGTGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGctctc                                                                                            <  2:455667/251‑1 (MQ=255)
                                                                                                                                                         cTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAAc                                                                                        >  1:161768/1‑251 (MQ=255)
                                                                                                                                                                     gCCCCGGTGGGACGCGGGAG‑AAAAATGGCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAACGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAAc                                                                            <  2:241576/251‑1 (MQ=255)
                                                                                                                                                                           gtggtACGCGGGAG‑AAAAATGTCACCATC‑GGCATCGCGCGTAAATTTCGTACTTTGGTCAATGACGATCC‑GATGTTCAAAGGAGTTTTAACCCGTGACGGGGATTACGTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGaa                                                                      <  1:291027/251‑1 (MQ=255)
                                                                                                                                                                                tacgcgggaccgaaaagTCACCATC‑GCAATCGCGAGTCAATAGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACACGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAACAGCCGATTTCCTCGTGTCGATTCCCGCTGATGACGCACCGAACAGCAGTGCGACTCGCGCGTCCGGATGGGGGCTATCTAACAGTCGCGCCACCAGTgcaaggg                                                                  >  2:195067/17‑245 (MQ=255)
                                                                                                                                                                                      gag‑aAAAATGTCACCATC‑GCCAGCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGTATTACTTTATCGCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAgctg                                                           <  1:443064/251‑1 (MQ=255)
                                                                                                                                                                                             aaTGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAACCGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGCACCGCAGTGCGACTGGCGCTTCCGTATGGGGGCTCTCTACCCGTCGCGCCAACAGAGAGATGGCCCGctggatggg                                                     >  2:8909/1‑250 (MQ=255)
                                                                                                                                                                                                ttCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCCCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCGCGGTTATGGGGCGCAGTGATGTGGCACGGAAGCAAAACGCCAATGTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACGGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGACAACAGTGAAATGGCCAGCTGGCTGGagt                                                   <  2:235257/249‑2 (MQ=255)
                                                                                                                                                                                                 gtccccacctccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGAACC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCGCGGGGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGagc                                                   <  2:189019/240‑1 (MQ=255)
                                                                                                                                                                                                           ccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGACATGGCCAGCTGGCTGGTGCAGCCCGAGaa                                         >  2:77184/1‑250 (MQ=255)
                                                                                                                                                                                                           ccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGaaa                                        <  1:79373/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGCTTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTa                                                                                          >  2:228403/1‑190 (MQ=255)
                                                                                                                                                                                                                      aaaTTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTa                                                                                          <  1:228403/190‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAgc                              <  1:130314/250‑1 (MQ=255)
                                                                                                                                                                                                                        aTTGCGTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCggggatgc                           >  2:211106/1‑245 (MQ=255)
                                                                                                                                                                                                                             gTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCtgggt                      <  1:452875/251‑1 (MQ=255)
                                                                                                                                                                                                                             gTACTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGACGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGCGCTGCtgggg                      >  1:40829/1‑250 (MQ=255)
                                                                                                                                                                                                                               aCTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGCTTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCAAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCACCAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCtgggtgg                    >  1:451952/1‑251 (MQ=255)
                                                                                                                                                                                                                               aCTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCtgggtgg                    >  2:140424/1‑251 (MQ=255)
                                                                                                                                                                                                                               aCTTTGCTCAATGACGATCC‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCtgggtgg                    >  1:68341/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ttGCTCAATGACGATCC‑GATGTGCAAAGGCGTTTTAACCTGTGACGGGGAGTACTTTATCTCGGTGATGGGGCGCAGTGATGGGGCACGTAAGCAAAACGCCAATTTCCTCGTGGCGAGGCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAAGGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTgggg                  <  2:295753/250‑1 (MQ=255)
                                                                                                                                                                                                                                        aaTGACGATCC‑GATGTCCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTCTCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGAAAAACGCCAATTTCCTCGTATCGATTCACTCTGATGCCGCACCGCACCGCGGTGCGACTGGCGCTTCCGTATGGGGGCTCTCTAACCGTCGCGCCAACAGTGAAATGGACAGCTGGCTGGAGCCGAACGAGACACAGTCGGAGCTGTTGGGAGGGGCggtggg           >  2:299852/1‑247 (MQ=255)
                                                                                                                                                                                                                                               gtCC‑GATGTTCAAAGGCGTTTTAACCCGGGACGGGGATTACTTTATCTCGGTGAGGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTg    <  2:18644/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 cc‑GATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGc  <  1:282979/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
GCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAG‑AAAAATGTCACCATC‑GCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCC‑GATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGC  >  NC_000913/4396525‑4397014

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