Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,378,778 G→T 100% S137* (TCG→TAG)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,378,7780GT95.6% 533.1 / 15.3 180S137* (TCG→TAG) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base T (98/74);  total (102/78)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 180 total aligned reads displayed.

CGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTG  >  CP009273/4378633‑4378924
                                                                                                                                                 |                                                                                                                                                  
cgcCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTaa                                                                                                                                                  <  2:1558420/148‑1 (MQ=255)
cgcCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTaa                                                                                                                                                  <  2:3321398/148‑1 (MQ=255)
  ccTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATc                                                                                                                                                >  2:2926723/1‑148 (MQ=255)
   cTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCg                                                                                                                                               <  2:738054/148‑1 (MQ=255)
    ttAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCgg                                                                                                                                              >  1:2183004/1‑148 (MQ=255)
    ttAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCgg                                                                                                                                              >  1:1891548/1‑148 (MQ=255)
     tAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTGGCTGTTGCGGCAGTTGATTCAGCTAATCGGc                                                                                                                                             >  1:952438/1‑148 (MQ=255)
      aaGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCa                                                                                                                                            <  2:1389962/148‑1 (MQ=255)
        gACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAAt                                                                                                                                          >  1:638076/1‑148 (MQ=255)
        gACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAAt                                                                                                                                          >  2:892453/1‑148 (MQ=255)
         aCGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAAtc                                                                                                                                         >  2:2561676/1‑148 (MQ=255)
            tGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTcg                                                                                                                                      <  1:2264884/148‑1 (MQ=255)
            tGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTcg                                                                                                                                      <  1:922549/148‑1 (MQ=255)
                     aGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGgcgcgc                                                                                                                             <  2:347484/148‑1 (MQ=255)
                          gCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTg                                                                                                                        >  1:1357926/1‑148 (MQ=255)
                            aaCGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGct                                                                                                                      >  2:117295/1‑148 (MQ=255)
                                cgAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTg                                                                                                                  <  1:957269/148‑1 (MQ=255)
                                cgAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTg                                                                                                                  <  1:1767436/148‑1 (MQ=255)
                                cgAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTg                                                                                                                  <  2:31923/148‑1 (MQ=255)
                                  aaaTTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCt                                                                                                                <  1:81161/148‑1 (MQ=255)
                                     ttttGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTggg                                                                                                             <  2:795871/148‑1 (MQ=255)
                                      tttGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGTCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCTGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGc                                                                                                            >  2:3377344/1‑148 (MQ=255)
                                      gttGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGc                                                                                                            <  2:2120039/147‑1 (MQ=255)
                                          tgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGa                                                                                                        >  2:1633457/1‑148 (MQ=255)
                                           gCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGAc                                                                                                       >  1:200592/1‑148 (MQ=255)
                                            ccTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACg                                                                                                      >  2:3071391/1‑148 (MQ=255)
                                             cTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACgg                                                                                                     <  2:2793512/148‑1 (MQ=255)
                                               gATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCt                                                                                                   >  1:2802794/1‑148 (MQ=255)
                                               gATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCt                                                                                                   >  2:2294940/1‑148 (MQ=255)
                                                 ttGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCttt                                                                                                 <  2:563954/148‑1 (MQ=255)
                                                  tGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCtttt                                                                                                >  2:302133/1‑148 (MQ=255)
                                                  tGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCtttt                                                                                                >  2:1112441/1‑148 (MQ=255)
                                                  tGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTAAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCtttt                                                                                                >  2:302130/1‑148 (MQ=255)
                                                       ggAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCa                                                                                           >  1:2135144/1‑148 (MQ=255)
                                                        gAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAg                                                                                          <  2:3132773/148‑1 (MQ=255)
                                                            aTTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAAc                                                                                      >  2:524042/1‑148 (MQ=255)
                                                            aTTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTGTATCCAGCacc                                                                                      >  2:2560890/1‑146 (MQ=255)
                                                                  ggTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgctg                                                                                <  1:2436729/148‑1 (MQ=255)
                                                                    tAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgatgct                                                                              >  2:2304471/1‑148 (MQ=255)
                                                                     aaGCGTTCCCAGGCGGCGCTCTATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGctgctg                                                                             <  2:2429359/148‑1 (MQ=255)
                                                                     aaGCGGTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGctgctg                                                                             <  2:1914778/148‑1 (MQ=255)
                                                                      aGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGa                                                                            >  1:3160347/1‑148 (MQ=255)
                                                                      aGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGa                                                                            >  1:2066262/1‑148 (MQ=255)
                                                                           tCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTTCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                       >  1:3501263/1‑148 (MQ=255)
                                                                           tCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                       <  1:1013436/148‑1 (MQ=255)
                                                                           tCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                       >  2:185826/1‑148 (MQ=255)
                                                                           tCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                       >  2:2157309/1‑148 (MQ=255)
                                                                           tCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                       <  1:1010993/148‑1 (MQ=255)
                                                                           tCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                       >  1:332995/1‑148 (MQ=255)
                                                                           tCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                       >  1:671728/1‑148 (MQ=255)
                                                                            cccAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGga                                                                      >  1:1348469/1‑148 (MQ=255)
                                                                            cccAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGga                                                                      <  2:2103204/148‑1 (MQ=255)
                                                                            cccAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGga                                                                      <  2:1729547/148‑1 (MQ=255)
                                                                               aGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGaa                                                                   >  2:1553275/1‑148 (MQ=255)
                                                                               aGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGaa                                                                   <  2:1003760/148‑1 (MQ=255)
                                                                                ggcggcGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAAt                                                                  <  1:124163/148‑1 (MQ=255)
                                                                                 gcggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAAtt                                                                 >  1:3113632/1‑148 (MQ=255)
                                                                                      gcTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTg                                                            >  2:2163569/1‑148 (MQ=255)
                                                                                        tCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAt                                                          >  1:1925397/1‑148 (MQ=255)
                                                                                        tCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAt                                                          >  2:75674/1‑148 (MQ=255)
                                                                                         cGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAtt                                                         >  1:2434223/1‑148 (MQ=255)
                                                                                            tctcATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAg                                                      >  1:1738876/1‑148 (MQ=255)
                                                                                            tctcATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAg                                                      >  1:1099666/1‑148 (MQ=255)
                                                                                             ctcATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgc                                                     >  2:2288575/1‑148 (MQ=255)
                                                                                               cATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgcgc                                                   >  2:2309763/1‑148 (MQ=255)
                                                                                               cATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgcgc                                                   >  2:3490241/1‑148 (MQ=255)
                                                                                                aTTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGcgcg                                                  <  1:205668/148‑1 (MQ=255)
                                                                                                 tttAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTTATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGt                                                 >  1:2364083/1‑148 (MQ=255)
                                                                                                 tttAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGt                                                 <  2:2183004/148‑1 (MQ=255)
                                                                                                 tttAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGt                                                 <  1:3401781/148‑1 (MQ=255)
                                                                                                    aaCTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCgg                                              >  1:244334/1‑148 (MQ=255)
                                                                                                    aaCTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCgg                                              >  1:2040064/1‑148 (MQ=255)
                                                                                                     aCTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGt                                             >  1:2680675/1‑148 (MQ=255)
                                                                                                      cTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTa                                            >  2:2016110/1‑148 (MQ=255)
                                                                                                       tGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGtcg                                           >  2:475195/1‑146 (MQ=255)
                                                                                                       tGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAg                                           >  2:475192/1‑148 (MQ=255)
                                                                                                         aCGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGAc                                         <  1:739036/148‑1 (MQ=255)
                                                                                                          cggcggGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACa                                        <  1:139505/148‑1 (MQ=255)
                                                                                                          cggcggGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACa                                        <  1:133244/148‑1 (MQ=255)
                                                                                                                    cGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGa                              <  2:2398383/148‑1 (MQ=255)
                                                                                                                    cGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGa                              >  1:484987/1‑148 (MQ=255)
                                                                                                                     ggTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGac                             <  2:1526249/148‑1 (MQ=255)
                                                                                                                        ttGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACg                          <  2:1160065/148‑1 (MQ=255)
                                                                                                                        ttGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACg                          <  1:2594014/148‑1 (MQ=255)
                                                                                                                        ttGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACg                          <  1:1257759/148‑1 (MQ=255)
                                                                                                                           cTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTg                       >  2:1791035/1‑148 (MQ=255)
                                                                                                                           cTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTg                       >  2:1951831/1‑148 (MQ=255)
                                                                                                                              ttGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTgcgc                    >  2:883339/1‑148 (MQ=255)
                                                                                                                               tGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGcgcg                   >  1:1455870/1‑148 (MQ=255)
                                                                                                                                gCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCgg                  <  1:1479873/148‑1 (MQ=255)
                                                                                                                                  ggCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCt                <  1:892402/148‑1 (MQ=255)
                                                                                                                                      gTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTc            >  2:1913652/1‑148 (MQ=255)
                                                                                                                                      gTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTc            <  1:2561676/148‑1 (MQ=255)
                                                                                                                                       ttGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCa           <  2:638076/148‑1 (MQ=255)
                                                                                                                                       ttGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCa           <  2:2608832/148‑1 (MQ=255)
                                                                                                                                           ttCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCa       >  2:3232280/1‑148 (MQ=255)
                                                                                                                                           ttCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCa       >  2:1156171/1‑148 (MQ=255)
                                                                                                                                           ttCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCa       <  1:2288575/148‑1 (MQ=255)
                                                                                                                                              aGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCAtgt    <  1:2664655/148‑1 (MQ=255)
                                                                                                                                              aGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCAtgt    <  1:1978767/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTG  >  CP009273/4378633‑4378924

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: