Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,442,595 A→C 6.0% K125Q (AAG→CAG)  PP_3060 → tail sheath protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,442,5950AC 6.0% 639.1 / 21.9 184K125Q (AAG→CAG) PP_3060tail sheath protein
Reads supporting (aligned to +/- strand):  ref base A (80/93);  new base C (4/7);  total (84/100)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 188 total aligned reads displayed.

TGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTCCCC  >  NC_002947/3442449‑3442744
                                                                                                                                                  |                                                                                                                                                     
tGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGc                                                                                                                                                    >  2:1496231/1‑150 (MQ=255)
 ggCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCa                                                                                                                                                   >  1:3270333/1‑150 (MQ=255)
  gCTGACGGCAACGCCTTCAACGGCGTATTCTGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCAc                                                                                                                                                  >  1:2407356/1‑150 (MQ=255)
    tGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACgg                                                                                                                                                >  1:2479478/1‑150 (MQ=255)
     gACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGc                                                                                                                                               <  1:691413/150‑1 (MQ=255)
     gACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGc                                                                                                                                               <  1:343299/150‑1 (MQ=255)
        ggCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGTGCTTCAGCTCGGCCtt                                                                                                                                            <  2:2784890/150‑1 (MQ=255)
        ggCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCAtt                                                                                                                                            <  1:127191/150‑1 (MQ=255)
        ggCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCAtt                                                                                                                                            <  2:2975884/150‑1 (MQ=255)
         gCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGTGCTTCAGCTCGGCCTTg                                                                                                                                           <  2:629530/150‑1 (MQ=255)
           aaCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                         >  2:4287072/1‑150 (MQ=255)
           aaCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                         >  1:735349/1‑150 (MQ=255)
           aaCGCCTTCAACGGCGTATTCCGAGGCTCGGTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                         <  2:1981181/150‑1 (MQ=255)
           aaCGCCTTCAACGGCGTATTCCGAGGCTCGGTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                         <  2:3548369/150‑1 (MQ=255)
              gCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGGCAATGCCGTGCTTCAGCTCGGCCTTGTCgcc                                                                                                                                      <  1:1988104/150‑3 (MQ=255)
                cTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCAc                                                                                                                                    <  1:2931290/150‑1 (MQ=255)
                cTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCAc                                                                                                                                    <  1:2299077/150‑1 (MQ=255)
                  tCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACct                                                                                                                                  <  2:4512208/150‑1 (MQ=255)
                    aaCGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACctac                                                                                                                                <  1:1430645/150‑1 (MQ=255)
                       ggCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCtactac                                                                                                                             <  2:410338/150‑1 (MQ=255)
                           tATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGc                                                                                                                         <  2:3186130/150‑1 (MQ=255)
                            aTTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCt                                                                                                                        <  1:1008313/150‑1 (MQ=255)
                             ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                       <  1:3546598/150‑1 (MQ=255)
                             ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                       <  1:801025/150‑1 (MQ=255)
                             ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                       <  1:3126810/150‑1 (MQ=255)
                             ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                       <  1:744612/150‑1 (MQ=255)
                             ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                       <  1:2767545/150‑1 (MQ=255)
                               ccGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGa                                                                                                                     <  2:4181332/150‑1 (MQ=255)
                                cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                    <  1:1205210/150‑1 (MQ=255)
                                cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                    <  1:528410/150‑1 (MQ=255)
                                cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                    >  2:4169035/1‑150 (MQ=255)
                                cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                    <  1:3374467/150‑1 (MQ=255)
                                cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                    <  2:238081/150‑1 (MQ=255)
                                  aGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGt                                                                                                                  >  1:2119588/1‑150 (MQ=255)
                                      tCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTcggc                                                                                                              >  1:4353405/1‑150 (MQ=255)
                                             aGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAgg                                                                                                       >  2:4340520/1‑150 (MQ=255)
                                                tCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCa                                                                                                    >  1:3446014/1‑150 (MQ=255)
                                                   aGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGTGCTTCAGCTCgg                                                                                                                                                <  2:2739366/103‑1 (MQ=255)
                                                   aGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGTGCTTCAGCTCgg                                                                                                                                                >  1:2739366/1‑103 (MQ=255)
                                                       gggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                             <  2:1797973/150‑1 (MQ=255)
                                                       gggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                             <  2:4485473/150‑1 (MQ=255)
                                                       gggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                             <  1:1107977/150‑1 (MQ=255)
                                                           gAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGAtt                                                                                         <  1:1407320/150‑1 (MQ=255)
                                                                ccTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGAccc                                                                                    <  1:1750785/150‑1 (MQ=255)
                                                                  ttGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCgg                                                                                  >  1:3991701/1‑150 (MQ=255)
                                                                     cGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCg                                                                               >  1:1946024/1‑150 (MQ=255)
                                                                       gtcgtcgtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGc                                                                             <  1:2032983/150‑1 (MQ=255)
                                                                        tcgtcgtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCa                                                                            >  1:3441057/1‑150 (MQ=255)
                                                                             gtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAg                                                                       >  2:527109/1‑150 (MQ=255)
                                                                               cACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCg                                                                     >  1:449978/1‑150 (MQ=255)
                                                                                aCCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCgt                                                                    >  2:829856/1‑150 (MQ=255)
                                                                                  ccTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCgtgt                                                                  <  2:1282389/150‑1 (MQ=255)
                                                                                   cTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGtc                                                                 <  2:3285816/150‑1 (MQ=255)
                                                                                        cgGCACCCAGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAATCTAGAAGTCGGAGGCGAGGTCATATACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCaa                                                            <  1:922038/150‑1 (MQ=255)
                                                                                         ggCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAAc                                                           >  1:2056508/1‑150 (MQ=255)
                                                                                            aCCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                    >  2:3406191/1‑90 (MQ=255)
                                                                                              ccTGAAAGGGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCAc                                                      <  1:49376/150‑1 (MQ=255)
                                                                                              ccTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATTAAGCGTGTCATCAACGGCAc                                                      <  1:2792115/150‑1 (MQ=255)
                                                                                                 gAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACTGa                                                   >  1:3701962/1‑150 (MQ=255)
                                                                                                 gAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGa                                                   <  2:2123474/150‑1 (MQ=255)
                                                                                                       gCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGCCATGAAGCGTGTCATCAACGGCACCGACCAACt                                             >  2:3014743/1‑150 (MQ=255)
                                                                                                         tGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTgg                                           >  1:3194182/1‑150 (MQ=255)
                                                                                                         tGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTgg                                           >  1:3810943/1‑150 (MQ=255)
                                                                                                          ggACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGc                                          <  2:2876627/150‑1 (MQ=255)
                                                                                                           gACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGcc                                         >  1:1493128/1‑150 (MQ=255)
                                                                                                           gACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGcc                                         <  1:1772209/150‑1 (MQ=255)
                                                                                                            aCATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCa                                        <  2:2423059/150‑1 (MQ=255)
                                                                                                            aCATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCa                                        <  1:2689569/150‑1 (MQ=255)
                                                                                                            aCATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCa                                        <  1:3977054/150‑1 (MQ=255)
                                                                                                            aCATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCa                                        <  2:3775825/150‑1 (MQ=255)
                                                                                                             cATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAg                                       >  1:1516681/1‑150 (MQ=255)
                                                                                                                gggCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCa                                    >  2:1733599/1‑150 (MQ=255)
                                                                                                                gggCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCa                                    >  2:2155125/1‑150 (MQ=255)
                                                                                                                gggCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCa                                    >  1:1533470/1‑150 (MQ=255)
                                                                                                                  gCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAgc                                  <  1:2936607/150‑1 (MQ=255)
                                                                                                                  gCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAgc                                  <  2:1946024/150‑1 (MQ=255)
                                                                                                                      cTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAgcgcgc                              <  2:3782309/150‑1 (MQ=255)
                                                                                                                      cTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAgcgcgc                              <  1:3824319/150‑1 (MQ=255)
                                                                                                                       tGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGcc                             >  2:1338442/1‑150 (MQ=255)
                                                                                                                       tGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGcc                             <  2:1236700/150‑1 (MQ=255)
                                                                                                                         gAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGa                           >  2:4289671/1‑150 (MQ=255)
                                                                                                                          aaGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCCCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGAc                          <  1:1789658/150‑1 (MQ=255)
                                                                                                                             gCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTc                       <  2:1716837/150‑1 (MQ=255)
                                                                                                                               cGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCgg                     >  2:2966000/1‑150 (MQ=255)
                                                                                                                                ggCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGc                    <  2:1127949/150‑1 (MQ=255)
                                                                                                                                ggCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGc                    <  1:642988/150‑1 (MQ=255)
                                                                                                                                   ggcAATGCCGCGCTTCAGCTCGGCCTTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTg                 >  2:2613609/3‑150 (MQ=255)
                                                                                                                                   gACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTg                 >  1:2183234/1‑150 (MQ=255)
                                                                                                                                    aCAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGt                >  2:4392244/1‑150 (MQ=255)
                                                                                                                                     cAATGCCGTGCTTCAGCTCGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTa               >  1:3911868/1‑150 (MQ=255)
                                                                                                                                      aaGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTaa              >  2:3974025/1‑150 (MQ=255)
                                                                                                                                      aaGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTaa              <  1:3042061/150‑1 (MQ=255)
                                                                                                                                         gCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAccc           >  1:1604712/1‑150 (MQ=255)
                                                                                                                                           cGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAccccc         >  1:3590329/1‑150 (MQ=255)
                                                                                                                                           cGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAccccc         <  2:4380234/150‑1 (MQ=255)
                                                                                                                                              gCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCttt      <  1:4169035/150‑1 (MQ=255)
                                                                                                                                               cTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTc     >  2:369901/1‑150 (MQ=255)
                                                                                                                                                tGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGAAACCCCCTATcc    >  1:3602163/1‑150 (MQ=255)
                                                                                                                                                 gAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTccc   <  2:4113721/150‑1 (MQ=255)
                                                                                                                                                  aaGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTcccc  >  2:3099235/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTCCCC  >  NC_002947/3442449‑3442744

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