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,7780GT99.7% 1040.3 / ‑4.5 317S137* (TCG→TAG) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (185/131);  total (186/131)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 317 total aligned reads displayed.

AGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTT  >  CP009273/4378631‑4378923
                                                                                                                                                   |                                                                                                                                                 
 gcgcCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAAt                                                                                                                                                <  2:1494436/150‑1 (MQ=255)
  cgcCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATc                                                                                                                                               <  2:502621/150‑1 (MQ=255)
    ccTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCgg                                                                                                                                             >  2:308363/1‑150 (MQ=255)
    ccTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCgg                                                                                                                                             <  1:1180975/150‑1 (MQ=255)
    ccTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCgg                                                                                                                                             >  2:3855307/1‑150 (MQ=255)
     cTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGc                                                                                                                                            >  2:3195695/1‑150 (MQ=255)
       tAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCaa                                                                                                                                          >  2:654123/1‑150 (MQ=255)
        aaGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAAt                                                                                                                                         >  2:2349556/1‑150 (MQ=255)
          gACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAAtct                                                                                                                                       >  2:5091242/1‑150 (MQ=255)
          gACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAAtca                                                                                                                                       <  1:1290652/150‑2 (MQ=255)
            cGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGTGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTcg                                                                                                                                     <  1:5109084/150‑1 (MQ=255)
             gTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTcgc                                                                                                                                    >  2:130312/1‑150 (MQ=255)
             gTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTcgc                                                                                                                                    <  1:3813540/150‑1 (MQ=255)
                   agTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGgcgc                                                                                                                              >  1:4637106/1‑150 (MQ=255)
                    gTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGgcgcg                                                                                                                             >  2:1347930/1‑150 (MQ=255)
                     tCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGgcgcgc                                                                                                                            >  2:2723759/1‑150 (MQ=255)
                       aGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTc                                                                                                                          >  2:57898/1‑150 (MQ=255)
                        gACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTcc                                                                                                                         >  1:3424454/1‑150 (MQ=255)
                         aCTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCt                                                                                                                        <  1:886687/150‑1 (MQ=255)
                         aCTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCt                                                                                                                        <  1:4137223/150‑1 (MQ=255)
                         aCTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCt                                                                                                                        <  2:1682927/150‑1 (MQ=255)
                          cTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTg                                                                                                                       >  2:3997082/1‑150 (MQ=255)
                            gCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGct                                                                                                                     >  1:5021212/1‑150 (MQ=255)
                             cAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGGGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCATTCTCGCGGGCGCGCTCCTGctc                                                                                                                    <  1:1638479/150‑1 (MQ=255)
                               aCGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCtt                                                                                                                  <  2:2963467/150‑1 (MQ=255)
                                cgcgAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTg                                                                                                                 <  1:846489/150‑1 (MQ=255)
                                  cgAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCt                                                                                                               <  1:3722416/150‑1 (MQ=255)
                                     aaTTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCATGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTggg                                                                                                            >  1:811918/1‑150 (MQ=255)
                                       ttttGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGcc                                                                                                          <  1:698157/150‑1 (MQ=255)
                                           gtgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGAc                                                                                                      >  2:4727127/1‑150 (MQ=255)
                                           gtgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGAc                                                                                                      <  2:2296666/150‑1 (MQ=255)
                                           gtgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGAc                                                                                                      <  1:549774/150‑1 (MQ=255)
                                                    tGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTAt                                                                                             >  2:5110417/1‑150 (MQ=255)
                                                    tGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTAt                                                                                             >  2:2032803/1‑150 (MQ=255)
                                                     gAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGCTCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATc                                                                                            <  2:1029741/150‑1 (MQ=255)
                                                     gAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATc                                                                                            >  2:1439371/1‑150 (MQ=255)
                                                         ggAGTAGTGCCGGTAAGCGGTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGc                                                                                        <  2:5268789/150‑1 (MQ=255)
                                                             tATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACt                                                                                    >  1:4422306/1‑150 (MQ=255)
                                                               ttGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTgg                                                                                  >  2:1647232/1‑150 (MQ=255)
                                                                 gCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgct                                                                                <  1:3262569/150‑1 (MQ=255)
                                                                  ccGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgctg                                                                               >  1:20257/1‑150 (MQ=255)
                                                                    ggTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgctgct                                                                             <  1:4911398/150‑1 (MQ=255)
                                                                    ggTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGgctgct                                                                             <  1:2266491/150‑1 (MQ=255)
                                                                      tAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGa                                                                           >  2:3184311/1‑150 (MQ=255)
                                                                      tAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGa                                                                           >  1:2169107/1‑150 (MQ=255)
                                                                       aaGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGAc                                                                          >  2:3516751/1‑150 (MQ=255)
                                                                          cGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTg                                                                       <  1:3074671/150‑1 (MQ=255)
                                                                           gTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                      >  2:1838375/1‑150 (MQ=255)
                                                                           gTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                      >  2:878275/1‑150 (MQ=255)
                                                                           gTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                      >  2:2448644/1‑150 (MQ=255)
                                                                           gTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTgg                                                                      >  1:3660509/1‑150 (MQ=255)
                                                                           gTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACAGGCTGCTGACCTgg                                                                      >  1:4983979/1‑150 (MQ=255)
                                                                                cAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAAt                                                                 >  1:199063/1‑150 (MQ=255)
                                                                                 aGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAAtt                                                                >  2:3615909/1‑150 (MQ=255)
                                                                                  ggcggcGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAAttt                                                               <  2:4625282/150‑1 (MQ=255)
                                                                                     ggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCt                                                            >  1:3668739/1‑150 (MQ=255)
                                                                                     ggcgCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCt                                                            >  1:158780/1‑150 (MQ=255)
                                                                                      gcgcTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTg                                                           <  2:1186055/150‑1 (MQ=255)
                                                                                        gcTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAt                                                         >  1:4482869/1‑150 (MQ=255)
                                                                                        gcTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAt                                                         >  1:1662628/1‑150 (MQ=255)
                                                                                        gcTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATAACCTGAt                                                         >  2:4893755/1‑150 (MQ=255)
                                                                                         cTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAtt                                                        >  1:4335789/1‑150 (MQ=255)
                                                                                         cTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAtt                                                        >  1:1543414/1‑150 (MQ=255)
                                                                                         cTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGAtt                                                        <  1:5050017/150‑1 (MQ=255)
                                                                                           cGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCa                                                      >  1:4569836/1‑150 (MQ=255)
                                                                                           cGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCa                                                      <  1:1059870/150‑1 (MQ=255)
                                                                                            gATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAg                                                     >  1:1008049/1‑150 (MQ=255)
                                                                                            gATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAg                                                     >  2:5152357/1‑150 (MQ=255)
                                                                                            gATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAg                                                     <  1:4516368/150‑1 (MQ=255)
                                                                                             aTCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgc                                                    <  1:2779030/150‑1 (MQ=255)
                                                                                              tctcATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgcg                                                   >  1:2970285/1‑150 (MQ=255)
                                                                                              tctcATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAgcg                                                   >  1:2569648/1‑150 (MQ=255)
                                                                                                tcATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGcgcg                                                 <  2:2844755/150‑1 (MQ=255)
                                                                                                 cATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCGTTTATCCAGCACCTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGt                                                >  2:4362791/1‑150 (MQ=255)
                                                                                                  aTTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTc                                               <  1:1246470/150‑1 (MQ=255)
                                                                                                  aTTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTc                                               >  1:173985/1‑150 (MQ=255)
                                                                                                     tAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGt                                            >  1:4070517/1‑150 (MQ=255)
                                                                                                          gACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACa                                       >  1:2081546/1‑150 (MQ=255)
                                                                                                            cggcggGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACAtt                                     <  1:3989738/150‑1 (MQ=255)
                                                                                                                    gTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGa                             >  1:1456054/1‑150 (MQ=255)
                                                                                                                    gTCGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGa                             >  2:4205472/1‑150 (MQ=255)
                                                                                                                      cGGTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGaca                           >  2:1546238/1‑150 (MQ=255)
                                                                                                                       ggTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGacac                          >  2:4481451/1‑150 (MQ=255)
                                                                                                                       ggTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGacac                          >  2:2634006/1‑150 (MQ=255)
                                                                                                                        gTTTGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACg                         <  2:779779/150‑1 (MQ=255)
                                                                                                                          ttGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCt                       >  2:3799961/1‑150 (MQ=255)
                                                                                                                          ttGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCt                       <  1:3466583/150‑1 (MQ=255)
                                                                                                                           tGCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTg                      <  2:2529000/150‑1 (MQ=255)
                                                                                                                            gCTGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTgc                     >  1:4669913/1‑150 (MQ=255)
                                                                                                                              tGTTGCGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTgcgc                   >  2:3401192/1‑150 (MQ=255)
                                                                                                                                   cGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTc              <  1:1905827/150‑1 (MQ=255)
                                                                                                                                   cGGCAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTc              >  1:1117727/1‑150 (MQ=255)
                                                                                                                                      cAGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTc           >  2:4314529/1‑150 (MQ=255)
                                                                                                                                       aGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCa          <  2:3878823/150‑1 (MQ=255)
                                                                                                                                       aGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCa          >  2:2192355/1‑150 (MQ=255)
                                                                                                                                       aGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCa          >  1:2693884/1‑150 (MQ=255)
                                                                                                                                       aGTTGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCa          >  1:2908397/1‑150 (MQ=255)
                                                                                                                                          tGATTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGc       >  1:4814291/1‑150 (MQ=255)
                                                                                                                                            aTTCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCAt     <  2:1073724/150‑1 (MQ=255)
                                                                                                                                             ttCAGCTAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCAtg    <  1:5091242/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGCGCCTTAAGACGTGCAGAGTCAGACTGCAACGCGAAATTTTGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTT  >  CP009273/4378631‑4378923

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: