Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,258,182 A→T 11.5% K281* (AAG→TAG)  setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,258,1820AT11.5% 845.3 / 83.6 261K281* (AAG→TAG) setBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base A (123/107);  new base T (18/12);  total (142/119)
Fisher's exact test for biased strand distribution p-value = 5.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 267 total aligned reads displayed.

GAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCGGGATGCTCTATTT  >  CP009273/2258034‑2258330
                                                                                                                                                    |                                                                                                                                                    
gAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTa                                                                                                                                                      <  1:2687506/149‑1 (MQ=255)
 aaCAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTaa                                                                                                                                                     >  1:4521074/1‑149 (MQ=255)
  aCAGCCTGTACATCATCAACATGCCTCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGAATGCCGGATATTTCGCCAAACGTCTGGGTAAg                                                                                                                                                    >  2:1698447/1‑149 (MQ=255)
   cAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGc                                                                                                                                                   <  2:3186339/149‑1 (MQ=255)
    aGCCTGTACATCATCAACATGCCGCGATGTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGCTATTTCGCCAAACGTCTGGGTCAGCg                                                                                                                                                  >  2:1067764/1‑149 (MQ=255)
      ccTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGtt                                                                                                                                                >  1:1194577/1‑149 (MQ=255)
         gTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTTAGCGTTTCt                                                                                                                                             <  2:3706118/149‑1 (MQ=255)
                atcaACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTTAGCGTTTCTTAATgcg                                                                                                                                      >  1:111841/1‑149 (MQ=255)
                atcaACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATgcg                                                                                                                                      >  1:240923/1‑149 (MQ=255)
                 tcaACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATgcgc                                                                                                                                     >  2:338876/1‑149 (MQ=255)
                 tcaACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATgcgc                                                                                                                                     >  1:1812299/1‑149 (MQ=255)
                  caACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGcgcg                                                                                                                                    >  1:3890671/1‑149 (MQ=255)
                   aaCATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGt                                                                                                                                   <  2:73944/149‑1 (MQ=255)
                    aCATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGtt                                                                                                                                  <  1:3813152/149‑1 (MQ=255)
                    aCATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGtt                                                                                                                                  >  1:2127718/1‑149 (MQ=255)
                    aCATGCCGCTATTTATTATCAACGAACTGCATCTTAACGAGAAACTGGCCGGTGTGATGAAGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCAGGAAATTTCGCCAAACGTCTGGGTAATCGTTTCTTAATGAGCGtt                                                                                                                                  >  2:1685823/1‑149 (MQ=255)
                      aTGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTtgc                                                                                                                                >  1:3584730/1‑149 (MQ=255)
                      aTGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTtgc                                                                                                                                <  2:408936/149‑1 (MQ=255)
                       tGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTtgct                                                                                                                               <  1:1769052/149‑1 (MQ=255)
                           gCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCg                                                                                                                           >  2:799928/1‑149 (MQ=255)
                               tttattATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCgtggg                                                                                                                       >  2:1837040/1‑149 (MQ=255)
                               tttattATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCgtggg                                                                                                                       >  2:4302620/1‑149 (MQ=255)
                               tttattATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCgtggg                                                                                                                       >  2:42059/1‑149 (MQ=255)
                                  attatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCgtgggtgg                                                                                                                    >  1:1693017/1‑149 (MQ=255)
                                   ttatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGc                                                                                                                   <  2:4660613/149‑1 (MQ=255)
                                    tatCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGGGGGTGGCg                                                                                                                  >  2:3287742/1‑149 (MQ=255)
                                      tCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTc                                                                                                                <  2:2031468/149‑1 (MQ=255)
                                       cAACGAACTGCATCGTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCt                                                                                                               <  1:2960789/149‑1 (MQ=255)
                                        aaCGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTTAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTg                                                                                                              >  1:2427061/1‑149 (MQ=255)
                                        aaCGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTg                                                                                                              >  1:3311429/1‑149 (MQ=255)
                                         aCGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGt                                                                                                             <  1:2125971/149‑1 (MQ=255)
                                              cTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTa                                                                                                        <  1:1233579/149‑1 (MQ=255)
                                               tGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTAc                                                                                                       >  1:1977408/1‑149 (MQ=255)
                                                  aTCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCa                                                                                                    >  2:3189394/1‑149 (MQ=255)
                                                    cTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAgg                                                                                                  <  2:2857557/149‑1 (MQ=255)
                                                     ttCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGa                                                                                                 >  1:3412913/1‑149 (MQ=255)
                                                         cGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGc                                                                                             >  2:4127560/1‑149 (MQ=255)
                                                          gagaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTTAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCt                                                                                            <  2:1092401/149‑1 (MQ=255)
                                                          gagaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCt                                                                                            <  1:2170969/149‑1 (MQ=255)
                                                          gagaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCt                                                                                            <  1:575784/149‑1 (MQ=255)
                                                          gagaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCt                                                                                            >  2:1135525/1‑149 (MQ=255)
                                                           agaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTTAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTg                                                                                           <  2:175756/149‑1 (MQ=255)
                                                            gaAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGa                                                                                          >  1:875832/1‑149 (MQ=255)
                                                               aCTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATgg                                                                                       <  2:3584724/149‑1 (MQ=255)
                                                               aCTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATgg                                                                                       <  2:3852909/149‑1 (MQ=255)
                                                               aCTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATgg                                                                                       <  2:2523989/149‑1 (MQ=255)
                                                                cTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGgc                                                                                      <  1:3910024/149‑1 (MQ=255)
                                                                cTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGgc                                                                                      <  2:4451984/149‑1 (MQ=255)
                                                                  ggCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTTAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGgcgc                                                                                    >  2:2983788/1‑149 (MQ=255)
                                                                       gtgtGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCa                                                                               >  2:3492266/1‑149 (MQ=255)
                                                                       gtgtGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCa                                                                               <  2:2104783/149‑1 (MQ=255)
                                                                       gtgtGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCa                                                                               >  2:3005657/1‑149 (MQ=255)
                                                                         gttatgatgGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCAcc                                                                             >  1:4370309/4‑149 (MQ=255)
                                                                             tgatgGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTc                                                                         <  1:2533141/149‑1 (MQ=255)
                                                                             tgatgGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTc                                                                         <  1:229277/149‑1 (MQ=255)
                                                                              gatgGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCa                                                                        >  1:3441117/1‑149 (MQ=255)
                                                                                 ggggACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTc                                                                     >  1:1039218/1‑149 (MQ=255)
                                                                                  gggACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCt                                                                    <  2:221974/149‑1 (MQ=255)
                                                                                     aCCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGtt                                                                 <  1:1568018/149‑1 (MQ=255)
                                                                                       cgccgccgGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTgg                                                               <  2:1919871/149‑1 (MQ=255)
                                                                                                gCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAgc                                                      <  2:1284320/149‑1 (MQ=255)
                                                                                                   ggAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAgctgc                                                   >  1:62757/1‑149 (MQ=255)
                                                                                                    gAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAgctgct                                                  <  1:1837982/149‑1 (MQ=255)
                                                                                                    gAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAgctgct                                                  <  1:4172770/149‑1 (MQ=255)
                                                                                                      aaTCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTaa                                                >  1:2979878/1‑149 (MQ=255)
                                                                                                       aTCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTaaa                                               <  1:2492316/149‑1 (MQ=255)
                                                                                                       aTCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTaaa                                               <  2:3383656/149‑1 (MQ=255)
                                                                                                        tCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAAt                                              >  1:1778044/1‑149 (MQ=255)
                                                                                                        tCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAAt                                              >  2:3349696/1‑149 (MQ=255)
                                                                                                        tCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAAt                                              <  2:4595041/149‑1 (MQ=255)
                                                                                                          ccGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGc                                            <  1:3529175/149‑1 (MQ=255)
                                                                                                           cgacgaTGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCt                                           <  2:4521074/149‑1 (MQ=255)
                                                                                                            gacgaTGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTa                                          >  2:2398437/1‑149 (MQ=255)
                                                                                                             acgaTGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTAt                                         <  1:3913525/149‑1 (MQ=255)
                                                                                                               gaTGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTAttt                                       <  1:4385122/149‑1 (MQ=255)
                                                                                                                    tGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTAtt                                  >  2:1678548/1‑149 (MQ=255)
                                                                                                                      aTTGCCGGATATTTCGCCAAACGTCTGGGTTAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTgg                                >  2:3423167/1‑149 (MQ=255)
                                                                                                                           cGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTc                           >  1:1971348/1‑149 (MQ=255)
                                                                                                                           cGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTc                           >  1:204383/1‑149 (MQ=255)
                                                                                                                             gATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTg                         >  2:4418164/1‑149 (MQ=255)
                                                                                                                              atatTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTgg                        >  2:3807830/1‑149 (MQ=255)
                                                                                                                               tatTTCGCCAAACGTCTGGGTTAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTggg                       >  1:3684891/1‑149 (MQ=255)
                                                                                                                               tatTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTggg                       >  1:3738614/1‑149 (MQ=255)
                                                                                                                                  ttCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGggcgg                    <  2:2746205/149‑1 (MQ=255)
                                                                                                                                   tCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGggcggc                   <  2:3524562/149‑1 (MQ=255)
                                                                                                                                     gCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCAt                 >  1:3043283/1‑149 (MQ=255)
                                                                                                                                      ccAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATc                <  2:3671182/149‑1 (MQ=255)
                                                                                                                                       cAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCg               <  1:4210348/149‑1 (MQ=255)
                                                                                                                                       cAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCg               <  1:406560/149‑1 (MQ=255)
                                                                                                                                        aaaCGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCgg              <  1:1079332/149‑1 (MQ=255)
                                                                                                                                        aaaCGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCgg              >  2:1909256/1‑149 (MQ=255)
                                                                                                                                         aaCGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTTGCATTATGGGCGGCATCggg             >  1:1446468/1‑149 (MQ=255)
                                                                                                                                           cGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCGGGAt           <  2:1437362/149‑1 (MQ=255)
                                                                                                                                            gTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCGGGATg          >  1:826008/1‑149 (MQ=255)
                                                                                                                                             tCTGGGTTAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCGGGATGc         <  1:4248286/149‑1 (MQ=255)
                                                                                                                                             tCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCGGGATGc         >  2:929284/1‑149 (MQ=255)
                                                                                                                                               tGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCGGGATGctc       >  1:542755/1‑149 (MQ=255)
                                                                                                                                                 ggTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCGGGATGCTCTa     >  2:3546194/1‑149 (MQ=255)
                                                                                                                                                  gTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCGGGATGCTCTAt    <  1:3274405/149‑1 (MQ=255)
                                                                                                                                                   tAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCGGGATGCTCTAtt   <  2:761853/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAACAGCCTGTACATCATCAACATGCCGCTATTTATTATCAACGAACTGCATCTTCCCGAGAAACTGGCCGGTGTGATGATGGGGACCGCCGCCGGGCTGGAAATCCCGACGATGTTGATTGCCGGATATTTCGCCAAACGTCTGGGTAAGCGTTTCTTAATGCGCGTTGCTGCCGTGGGTGGCGTCTGTTTTTACGCAGGAATGCTGATGGCGCATTCACCTGTCATTCTGTTGGGCTTGCAGCTGCTAAATGCTATTTTTATTGGCATTCTGGGCGGCATCGGGATGCTCTATTT  >  CP009273/2258034‑2258330

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: