Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 793,055 C→A 5.5% S349* (TCG→TAG)  modC → molybdate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273793,0550CA 5.5% 1072.1 / 25.5 272S349* (TCG→TAG) modCmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (117/140);  new base A (5/10);  total (122/150)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 281 total aligned reads displayed.

AGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCGGCTTCCAGCATACTG  >  CP009273/792906‑793204
                                                                                                                                                     |                                                                                                                                                     
aGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTCAAAGTGTGTc                                                                                                                                                       >  2:127733/1‑150 (MQ=255)
aGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTc                                                                                                                                                       <  2:2693778/150‑1 (MQ=255)
aGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTc                                                                                                                                                       <  1:2020350/150‑1 (MQ=255)
    gttAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATa                                                                                                                                                   >  1:1407119/1‑150 (MQ=255)
    gttAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTAGATa                                                                                                                                                   <  2:3219578/150‑1 (MQ=255)
       aaTAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAAcc                                                                                                                                                >  2:444066/1‑150 (MQ=255)
          aGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGcc                                                                                                                                             <  1:2137743/150‑1 (MQ=255)
           gTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                                                                                                                            >  1:655441/1‑150 (MQ=255)
            ttATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTg                                                                                                                                           >  1:1102588/1‑150 (MQ=255)
            ttATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTg                                                                                                                                           <  1:5134791/150‑1 (MQ=255)
               tGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGAtt                                                                                                                                        >  2:1422294/1‑150 (MQ=255)
                gacgacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTa                                                                                                                                       <  1:3430743/150‑1 (MQ=255)
                 acgacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTaa                                                                                                                                      <  2:4200767/150‑1 (MQ=255)
                  cgacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTaaa                                                                                                                                     >  1:2670455/1‑150 (MQ=255)
                   gacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAAt                                                                                                                                    >  2:1265119/1‑150 (MQ=255)
                   gacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGGAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAAt                                                                                                                                    <  2:4941112/150‑1 (MQ=255)
                      aaCGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAg                                                                                                                                 <  1:3080114/150‑1 (MQ=255)
                             aGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCta                                                                                                                          <  2:4172396/150‑1 (MQ=255)
                               gtggaagtggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCtata                                                                                                                        >  2:3656770/1‑150 (MQ=255)
                                    agtggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGa                                                                                                                   >  1:4877270/1‑150 (MQ=255)
                                       ggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACt                                                                                                                <  1:4802336/150‑1 (MQ=255)
                                       ggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACt                                                                                                                <  1:4971403/150‑1 (MQ=255)
                                         aaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTgg                                                                                                              >  1:3781891/1‑150 (MQ=255)
                                         aaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTgg                                                                                                              <  2:3009009/150‑1 (MQ=255)
                                         aaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTgg                                                                                                              <  2:2344368/150‑1 (MQ=255)
                                          aCTGGAAGTCGGCGGTAAACCGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTggg                                                                                                             <  2:2996694/150‑1 (MQ=255)
                                              gAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAAt                                                                                                         <  2:3333862/150‑1 (MQ=255)
                                              gAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAAt                                                                                                         >  2:1754354/1‑150 (MQ=255)
                                              gAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAAt                                                                                                         >  1:4186085/1‑150 (MQ=255)
                                               aaGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATg                                                                                                        >  1:4954000/1‑150 (MQ=255)
                                                  tCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTg                                                                                                     >  1:68099/1‑150 (MQ=255)
                                                  tCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTg                                                                                                     <  2:4643489/150‑1 (MQ=255)
                                                  tCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTAGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTg                                                                                                     <  2:2847087/150‑1 (MQ=255)
                                                   cggcggTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGt                                                                                                    <  1:1071134/150‑1 (MQ=255)
                                                     gcggTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCg                                                                                                  <  2:3915395/150‑1 (MQ=255)
                                                       ggTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCggt                                                                                                >  1:905480/1‑150 (MQ=255)
                                                       ggTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCggt                                                                                                >  1:3558908/1‑150 (MQ=255)
                                                          aaaaCGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCggtggt                                                                                             <  2:2073632/150‑1 (MQ=255)
                                                           aaaCGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGgtggtg                                                                                            <  2:2619107/150‑1 (MQ=255)
                                                            aaCGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGgtgt                                                                                           <  2:1420632/150‑1 (MQ=255)
                                                                 tgtgGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTACAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCa                                                                                      <  2:4009785/150‑1 (MQ=255)
                                                                 tgtgGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCa                                                                                      <  1:5007662/150‑1 (MQ=255)
                                                                 tgtgGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCa                                                                                      <  2:1199985/150‑1 (MQ=255)
                                                                    gggCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCa                                                                                   >  1:331118/1‑150 (MQ=255)
                                                                      gcgcGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAAt                                                                                 >  2:2716117/1‑150 (MQ=255)
                                                                      gcgcGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAAt                                                                                 >  1:3698693/1‑150 (MQ=255)
                                                                       cgcgTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTAGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATc                                                                                <  2:3662174/150‑1 (MQ=255)
                                                                          gTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATcaca                                                                             >  2:1953316/1‑150 (MQ=255)
                                                                          gTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATcaca                                                                             <  1:1313217/150‑1 (MQ=255)
                                                                          gTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATcaca                                                                             >  1:3635582/1‑150 (MQ=255)
                                                                             tCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATg                                                                          <  1:119329/150‑1 (MQ=255)
                                                                             tCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATg                                                                          >  1:4041956/1‑150 (MQ=255)
                                                                                 cccGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTgg                                                                      <  2:3800130/150‑1 (MQ=255)
                                                                                 cccGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTgg                                                                      >  1:1150468/1‑150 (MQ=255)
                                                                                       ggCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGgcgcgcg                                                                <  2:3853697/150‑1 (MQ=255)
                                                                                         ccAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCt                                                              >  2:39695/1‑150 (MQ=255)
                                                                                         ccAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCt                                                              >  2:682433/1‑150 (MQ=255)
                                                                                         ccAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCt                                                              >  2:2641212/1‑150 (MQ=255)
                                                                                          cAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCtt                                                             <  1:4787904/150‑1 (MQ=255)
                                                                                           aGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCttt                                                            >  2:1086267/1‑150 (MQ=255)
                                                                                            gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCtttt                                                           >  2:4504805/1‑150 (MQ=255)
                                                                                                 gAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACcgc                                                      <  1:1611398/150‑1 (MQ=255)
                                                                                                 gAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAAATGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACcgc                                                      >  2:3504710/1‑150 (MQ=255)
                                                                                                  aaCTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACcgcg                                                     <  2:2528302/150‑1 (MQ=255)
                                                                                                    cTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTc                                                   <  2:4139283/150‑1 (MQ=255)
                                                                                                    cTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTc                                                   <  2:3666171/150‑1 (MQ=255)
                                                                                                     tGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCa                                                  <  2:607026/150‑1 (MQ=255)
                                                                                                     tGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCa                                                  >  1:467935/1‑150 (MQ=255)
                                                                                                      ggCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCAt                                                 <  1:5044399/150‑1 (MQ=255)
                                                                                                       gCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATc                                                <  2:5061244/150‑1 (MQ=255)
                                                                                                        cGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCg                                               <  1:3338622/150‑1 (MQ=255)
                                                                                                          aTCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGc                                             >  2:1628948/1‑150 (MQ=255)
                                                                                                            cAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGt                                           >  2:748823/1‑150 (MQ=255)
                                                                                                             aaaCCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGtt                                          <  1:2161067/150‑1 (MQ=255)
                                                                                                                ccTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGcc                                       <  1:2489877/150‑1 (MQ=255)
                                                                                                                 cTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGccc                                      <  2:3237939/150‑1 (MQ=255)
                                                                                                                    gCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTGTACCGCGTCATCGGCGTTGCCCATc                                   <  2:4914323/150‑1 (MQ=255)
                                                                                                                                aCGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTa                       >  1:2434698/1‑150 (MQ=255)
                                                                                                                                 cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTAc                      >  2:3420392/1‑150 (MQ=255)
                                                                                                                                  gcgcAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTAcc                     <  2:931490/150‑1 (MQ=255)
                                                                                                                                  gcgcAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTAcc                     <  2:548917/150‑1 (MQ=255)
                                                                                                                                   cgcAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCa                    <  2:4732523/150‑1 (MQ=255)
                                                                                                                                   cgcAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCa                    >  1:1287367/1‑150 (MQ=255)
                                                                                                                                    gcAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAg                   <  1:71475/150‑1 (MQ=255)
                                                                                                                                    gcAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAg                   <  1:2882593/150‑1 (MQ=255)
                                                                                                                                     cAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGc                  >  2:433640/1‑150 (MQ=255)
                                                                                                                                     cAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGc                  <  1:3566797/150‑1 (MQ=255)
                                                                                                                                     cAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGc                  <  1:2969883/150‑1 (MQ=255)
                                                                                                                                     cAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGc                  >  2:920098/1‑150 (MQ=255)
                                                                                                                                      aaaTTAACAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAACTGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCg                 <  2:1607607/150‑1 (MQ=255)
                                                                                                                                       aaTTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCgg                <  2:226707/150‑1 (MQ=255)
                                                                                                                                           aaaaGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCGGCTTc            >  2:2126796/1‑150 (MQ=255)
                                                                                                                                            acaGTGTGTCGATAACCGCCTGAATAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCGGCTTcc           <  1:4535452/148‑1 (MQ=255)
                                                                                                                                             aaGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCGGCTTCCa          >  1:4082285/1‑150 (MQ=255)
                                                                                                                                              aGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCGGCTTCCaa         <  2:2009047/150‑2 (MQ=255)
                                                                                                                                              aGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCGGCTTCCAg         <  2:961052/150‑1 (MQ=255)
                                                                                                                                                 gtgtCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCGGCTTCCAGCAt      >  2:5099268/1‑150 (MQ=255)
                                                                                                                                                  tgtCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCGGCTTCCAGCATa     <  1:4257182/150‑1 (MQ=255)
                                                                                                                                                     cGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCGGCTTCCAGCATACTg  >  1:100415/1‑150 (MQ=255)
                                                                                                                                                     cGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCGGCTTCCAGCATACTg  <  1:444066/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTGTCGGTGGTGTTATCACCAATCACAATGTTGGCGCGCGCTTTTACCGCGTCATCGGCGTTGCCCATCGCCACGCCTGTACCAGCGGCTTCCAGCATACTG  >  CP009273/792906‑793204

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: ‑