Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,072,207 C→T 11.5% R52R (CGG→CGA Y75_RS10480 ← IS3‑like element IS2 family transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077792,072,2070CT11.5% 247.5 / 19.8 86R52R (CGG→CGAY75_RS10480IS3‑like element IS2 family transposase
Reads supporting (aligned to +/- strand):  ref base C (41/35);  new base T (7/3);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

TTGCCGAGCAGGCGCTGGAGTTCTTTAATCTGCTTCATGGCGGCAGCAAGTTCAGAGGCAGGAACAACCTGTTCTCCGGCGGCGACAGCAGTAAGACTTCCTTCCTGGTATTGCTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCAT  >  NC_007779/2072052‑2072442
                                                                                                                                                           |                                                                                                                                                                                                                                           
tgcccggggatggggagtaagcttagtcagccgttgtctccgggtgtcgcgtgcgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTATCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCg                                                                                                                                                          >  2:1730406‑M2/144‑239 (MQ=255)
tgcccggggatggggagtaagcttagtcagccgttgtctccgggtgtcgcgtgcgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTATCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCg                                                                                                                                                          >  2:2046344‑M2/144‑239 (MQ=255)
        gatggggagtaagcttagtcagccgttgtctccgggtgtcgcgtgcgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTATCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTcc                                                                                                                                                   <  1:1730406‑M2/103‑1 (MQ=255)
               agtaagcttagtcagccgttgtctccgggtgtcgcgtgcgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTa                                                                                                                                            <  2:2110616‑M2/110‑1 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAatcat                                                                                                                               >  2:337849‑M2/91‑213 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAatcat                                                                                                                               <  1:337849‑M2/123‑1 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGt                                                                                                                           <  2:410267‑M2/127‑1 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGt                                                                                                                           >  2:2139879‑M2/91‑217 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGt                                                                                                                           >  1:410267‑M2/91‑217 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGt                                                                                                                           <  2:1503740‑M2/127‑1 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGt                                                                                                                           <  2:2073420‑M2/127‑1 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGt                                                                                                                           >  1:1503740‑M2/91‑217 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGt                                                                                                                           <  1:2139879‑M2/127‑1 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGt                                                                                                                           >  1:2073420‑M2/91‑217 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTtct                                                                                                                        <  1:1034356‑M2/130‑1 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTtct                                                                                                                        >  2:1034356‑M2/91‑220 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTaaaaa                                                                                                     >  2:201141‑M2/91‑239 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTaaaa                                                                                                      >  1:1809992‑M2/91‑238 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTaaaa                                                                                                      >  2:1363800‑M2/91‑238 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTaaaa                                                                                                      >  1:1804259‑M2/91‑238 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTaaaa                                                                                                      >  1:724276‑M2/91‑238 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTaaaa                                                                                                      >  1:1538836‑M2/91‑238 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTATCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCCATCATCTGTTCTCCAATGACTAGCCTaaaa                                                                                                      >  2:1297740‑M2/91‑238 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTATCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTaaaa                                                                                                      >  1:627116‑M2/91‑238 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTATCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTaaaa                                                                                                      >  1:1879666‑M2/91‑238 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgggcttctccgtaTACACCATGTTATCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTaaaa                                                                                                      >  1:1059452‑M2/91‑238 (MQ=255)
                                                     cgcagttgcacgtcattctcagacgaaccgatgactggatggatggccgccgcagccgccacactgatgatacggatgagcttctccgtaTCCACCATGTTATCGGGCAACGAGGGAGACCGTAATCCCCGGTTCAAAGCTCTGCTGAACAACTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCAGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTaaaa                                                                                                      >  1:178364‑M2/91‑238 (MQ=255)
                                                                      ctcagacgtaccgatgactggctggatggtctccgcagtcgtcacactgatgctacggatgtgcttctccgtaTACACCATGTTATCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATc                                                                                    <  2:1059452‑M2/166‑1 (MQ=255)
                                                                      ctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATc                                                                                    <  2:724276‑M2/166‑1 (MQ=255)
                                                                      ctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATc                                                                                    <  1:1363800‑M2/166‑1 (MQ=255)
                                                                      ctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATc                                                                                    <  1:201141‑M2/166‑1 (MQ=255)
                                                                      ctcagacgaaccgatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTATCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATc                                                                                    <  2:1879666‑M2/166‑1 (MQ=255)
                                                                                  gatgactggatggatggccgccgcagtcgtcacactgatgatacggatgtgcttctccgtaTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGt                                                                        <  2:1538836‑M2/178‑1 (MQ=255)
                                                                                                                           agaAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGc                               >  1:447276/1‑239 (MQ=31)
                                                                                                                                      gctggctgCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCAtcct                     >  2:1217833/1‑238 (MQ=35)
                                                                                                                                        tggctgCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTACGACATCAATCATCTGTTCTCCAATGCCTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGTGGCCAGTCTAGTGGCGAAGCATCCTccc                  >  2:2166930/1‑239 (MQ=255)
                                                                                                                                        tggctgCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCAtcctcc                   <  1:213546/238‑1 (MQ=37)
                                                                                                                                         ggctgCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCg                 <  1:539227/239‑1 (MQ=38)
                                                                                                                                         ggctgCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCg                 >  1:527473/1‑239 (MQ=38)
                                                                                                                                         ggctgCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCg                 <  1:1217833/239‑1 (MQ=38)
                                                                                                                                         cgctggtACACCATGTTGCCAGGCATCGAGGGAGCCCGTCATCCCCGGTTCGAAGCTCCGCTGAACGATTGCGAGCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTGCGGCCCTAAGACAGCAATCATCTGTTCTCCAATGACTAGTCTATAATCTAGGATTAGGACTATCCCTTATTTAGGTGATATTCGTTGTCTGGAGAGTCACGGGTCCCATCTAGTGGCGGAGTGTCCTTCCg                 <  2:1104117/238‑1 (MQ=11)
                                                                                                                                             gctACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGtt             <  2:1230894/239‑1 (MQ=40)
                                                                                                                                                   CCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAGCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAAtt       <  1:527504‑M2/239‑1 (MQ=255)
                                                                                                                                                    CATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTa                                                                                               >  2:2137158‑M2/1‑150 (MQ=255)
                                                                                                                                                    CATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTa                                                                                               <  1:2137158‑M2/150‑1 (MQ=255)
                                                                                                                                                    CATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTAtt                                                                                             <  2:280051‑M2/152‑1 (MQ=255)
                                                                                                                                                    CATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTAtt                                                                                             >  1:280051‑M2/1‑152 (MQ=255)
                                                                                                                                                     ATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATtc      >  2:2017266‑M2/1‑238 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGTAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGTCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGCGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCAGTGTTTTAattcgc    >  2:2144220‑M2/1‑237 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTg                                                                                                                                                                                         <  2:553513‑M2/58‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTg                                                                                                                                                                                         >  1:553513‑M2/1‑58 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACatcaat                                                                                                                                  >  2:240618‑M2/1‑113 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACatcaat                                                                                                                                  <  1:240618‑M2/113‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGtt                                                                                                                          <  1:61960‑M2/121‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGtt                                                                                                                          >  2:61960‑M2/1‑121 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTcc                                                                                                                      <  2:238214‑M2/125‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTcc                                                                                                                      >  1:238214‑M2/1‑125 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGAc                                                                                                                <  1:1371572‑M2/131‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGAc                                                                                                                >  2:1371572‑M2/1‑131 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTa                                                                                               >  1:1235606‑M2/1‑148 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTa                                                                                               <  2:1235606‑M2/148‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTAtt                                                                                             >  1:223763‑M2/1‑150 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTAtt                                                                                             <  2:223763‑M2/150‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCCCTTATTTAAGTGATATTGGTTGTCTGGCGATTCAGGGGGCCAGTCTAGTGGCGAAGCATACTCCCGTGTTTTAAttca     >  2:1374740‑M2/1‑237 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTaa                                                                          >  1:2022632‑M2/1‑169 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTaa                                                                          <  2:2022632‑M2/169‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGc                                           <  1:1038846‑M2/200‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGc                                           >  2:804264‑M2/1‑200 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGc                                           >  2:1038846‑M2/1‑200 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGc                                           >  1:1364617‑M2/1‑200 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGc                                           <  1:159386‑M2/200‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGc                                           >  2:159386‑M2/1‑200 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGc                                           <  2:1364617‑M2/200‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGc                                           <  1:804264‑M2/200‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTa                                    >  1:1019329‑M2/1‑207 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTa                                    <  2:1019329‑M2/207‑1 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATtctc    >  1:848257‑M2/1‑239 (MQ=255)
                                                                                                                                                      TGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATtctc    >  1:564216‑M2/1‑239 (MQ=255)
                                                                                                                                                        TTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAatcat                                                                                                                               >  1:782561‑M2/1‑114 (MQ=255)
                                                                                                                                                        TTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAatcat                                                                                                                               <  2:782561‑M2/114‑1 (MQ=255)
                                                                                                                                                        TTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGc                                           <  1:990737‑M2/198‑1 (MQ=255)
                                                                                                                                                        TTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGc                                           >  2:499927‑M2/1‑198 (MQ=255)
                                                                                                                                                        TTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGc                                           <  1:499927‑M2/198‑1 (MQ=255)
                                                                                                                                                        TTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGc                                           >  2:990737‑M2/1‑198 (MQ=255)
                                                                                                                                                        TTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCAt  >  1:1570250‑M2/1‑239 (MQ=255)
                                                                                                                                                        TTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCAt  >  1:1126648‑M2/1‑239 (MQ=255)
                                                                                                                                                           |                                                                                                                                                                                                                                           
TTGCCGAGCAGGCGCTGGAGTTCTTTAATCTGCTTCATGGCGGCAGCAAGTTCAGAGGCAGGAACAACCTGTTCTCCGGCGGCGACAGCAGTAAGACTTCCTTCCTGGTATTGCTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCAT  >  NC_007779/2072052‑2072442

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