Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,233,142 A→T 9.8% V229E (GTG→GAG)  mglB ← methyl‑galactoside transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,233,1420AT 9.8% 1022.3 / 78.7 312V229E (GTG→GAG) mglBmethyl‑galactoside transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (114/166);  new base T (15/16);  total (129/183)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 320 total aligned reads displayed.

TCAGTACGGTGCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCT  >  CP009273/2232993‑2233291
                                                                                                                                                     |                                                                                                                                                     
tCAGTACGGTGCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCa                                                                                                                                                       >  2:1525660/1‑150 (MQ=255)
   gTACGGTGCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACtt                                                                                                                                                    >  2:2933745/1‑150 (MQ=255)
    tACGGTGCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTc                                                                                                                                                   <  1:746721/150‑1 (MQ=255)
     aCGGTGCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCg                                                                                                                                                  <  2:1765855/150‑1 (MQ=255)
     aCGGTGCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCg                                                                                                                                                  >  1:2525917/1‑150 (MQ=255)
         tGCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGAttt                                                                                                                                              >  2:3284953/1‑150 (MQ=255)
         tGCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGAttt                                                                                                                                              >  1:1240797/1‑150 (MQ=255)
          gCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCTCTTCGAtttt                                                                                                                                             >  2:3780003/1‑150 (MQ=255)
          gCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCTCTTCGAtttt                                                                                                                                             >  1:3154442/1‑150 (MQ=255)
          gCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCTCTTCGAtttt                                                                                                                                             >  1:2568736/1‑150 (MQ=255)
             cGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGCTAACCACTTCGATTttgtt                                                                                                                                          >  2:3077201/1‑150 (MQ=255)
              gCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTttgttg                                                                                                                                         <  2:3703182/150‑1 (MQ=255)
              gCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTttgttg                                                                                                                                         <  1:4239725/150‑1 (MQ=255)
               ccAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTgg                                                                                                                                        <  1:2830188/150‑1 (MQ=255)
                 aGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCTCTTCGATTTTGTTGGCg                                                                                                                                      <  1:2216193/150‑1 (MQ=255)
                 aGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCg                                                                                                                                      <  1:3538188/150‑1 (MQ=255)
                   tGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGtt                                                                                                                                    <  2:4125957/150‑1 (MQ=255)
                     cACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCTCTTCGATTTTGTTGGCGTTCg                                                                                                                                  <  1:3066264/150‑1 (MQ=255)
                      aCCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCgg                                                                                                                                 >  1:4153205/1‑150 (MQ=255)
                      aCCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCgg                                                                                                                                 <  1:18923/150‑1 (MQ=255)
                       ccGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCggg                                                                                                                                >  1:2453062/1‑150 (MQ=255)
                        cGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGc                                                                                                                               <  1:4386104/150‑1 (MQ=255)
                            tttCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGa                                                                                                                           >  2:1265894/1‑150 (MQ=255)
                             ttCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGAc                                                                                                                          <  2:2651135/150‑1 (MQ=255)
                                aCCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGc                                                                                                                       <  2:164955/150‑1 (MQ=255)
                                 ccagcgccagcgCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGcc                                                                                                                      >  2:3614673/1‑150 (MQ=255)
                                 ccagcgccagcgCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGcc                                                                                                                      <  1:4632975/150‑1 (MQ=255)
                                    gcgccagcgcTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCTCTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAgg                                                                                                                   <  1:4783997/150‑1 (MQ=255)
                                    gcgccagcgcTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAgg                                                                                                                   >  1:1169590/1‑150 (MQ=255)
                                        cagcgcTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTc                                                                                                               >  1:2497976/1‑150 (MQ=255)
                                        cagcgcTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGAGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTc                                                                                                               >  2:4560322/1‑150 (MQ=255)
                                         agcgcTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTcc                                                                                                              >  1:3525610/1‑150 (MQ=255)
                                           cgcTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCat                                                                                                            <  1:2621901/150‑1 (MQ=255)
                                             cTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCatct                                                                                                          >  2:2452366/1‑150 (MQ=255)
                                                  ggCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCatcttatct                                                                                                     <  2:2803466/150‑1 (MQ=255)
                                                          aTCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTg                                                                                             >  2:3025691/1‑150 (MQ=255)
                                                            cGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAg                                                                                           >  2:2365054/1‑150 (MQ=255)
                                                                gCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCTCTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTCTCTTTCGCCTGAGCGgt                                                                                       <  2:4232237/150‑1 (MQ=255)
                                                                 ccAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGgtg                                                                                      <  2:3227293/150‑1 (MQ=255)
                                                                 ccAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGgtg                                                                                      <  2:840657/150‑1 (MQ=255)
                                                                     acacCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTccc                                                                                  >  2:3319345/1‑150 (MQ=255)
                                                                      cacCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCca                                                                                 <  1:114210/150‑1 (MQ=255)
                                                                       acCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCCCTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTGATCTTTCGCCTGAGCGGTGTCCcac                                                                                <  2:220495/150‑1 (MQ=255)
                                                                        ccGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCcaca                                                                               <  1:3019487/150‑1 (MQ=255)
                                                                          ggAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACAtt                                                                             <  1:2104533/150‑1 (MQ=255)
                                                                          ggAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACAtt                                                                             <  1:270605/150‑1 (MQ=255)
                                                                          ggAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACAtt                                                                             <  1:1688748/150‑1 (MQ=255)
                                                                           gAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTg                                                                            <  2:2409263/150‑1 (MQ=255)
                                                                           gAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATCGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTg                                                                            <  2:3908892/150‑1 (MQ=255)
                                                                            aaTGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGc                                                                           <  2:2102072/150‑1 (MQ=255)
                                                                                cTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTa                                                                       >  2:642567/1‑150 (MQ=255)
                                                                                 tGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCTGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTAt                                                                      <  2:147370/150‑1 (MQ=255)
                                                                                 tGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTAt                                                                      <  1:1265955/150‑1 (MQ=255)
                                                                                  ggACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATc                                                                     <  2:691284/150‑1 (MQ=255)
                                                                                       tgttgtGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCtaact                                                                <  2:1227144/150‑1 (MQ=255)
                                                                                        gttgtGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGTTAACGACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGGCTGAGCGGTGTCCCACATTGCGGTATCtaactg                                                               <  1:2936610/150‑1 (MQ=255)
                                                                                         ttgtGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCtaactgt                                                              <  2:3462013/150‑1 (MQ=255)
                                                                                          tgtgtgCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCtaactgta                                                             >  2:757399/1‑150 (MQ=255)
                                                                                              tgCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACAGTGCGGTATCtaactgtaactg                                                         >  2:4097303/1‑150 (MQ=255)
                                                                                                 tttCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCTCTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTc                                                      >  2:2988296/1‑150 (MQ=255)
                                                                                                  ttCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCa                                                     <  2:2717997/150‑1 (MQ=255)
                                                                                                   tCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAg                                                    >  2:1842746/1‑150 (MQ=255)
                                                                                                    cAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGGGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGt                                                   >  2:1928029/1‑150 (MQ=255)
                                                                                                     aGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGtt                                                  <  2:4388155/150‑1 (MQ=255)
                                                                                                        gcTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCTCTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTg                                               >  1:1806432/1‑150 (MQ=255)
                                                                                                         cTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGa                                              >  2:1674619/1‑150 (MQ=255)
                                                                                                          ttCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGAt                                             <  2:3985686/150‑1 (MQ=255)
                                                                                                          ttCAACCGCGCCCAGTGCCATCGCATCGTTGTGGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGAt                                             <  2:2831616/150‑1 (MQ=255)
                                                                                                              aCCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCt                                         <  1:1993302/150‑1 (MQ=255)
                                                                                                                    cccATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCg                                   >  2:954026/1‑150 (MQ=255)
                                                                                                                    cccATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCg                                   <  1:3978620/150‑1 (MQ=255)
                                                                                                                    cccATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCg                                   >  1:864430/1‑150 (MQ=255)
                                                                                                                     ccATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGt                                  <  2:1181105/150‑1 (MQ=255)
                                                                                                                     ccATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGt                                  >  1:4107599/1‑150 (MQ=255)
                                                                                                                      cATTGCCATCGCATCGTTGTTGGCGATAACCTCTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGtt                                 <  1:3446411/150‑1 (MQ=255)
                                                                                                                      cATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGtt                                 <  1:321033/150‑1 (MQ=255)
                                                                                                                      cATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGtt                                 <  2:2508050/150‑1 (MQ=255)
                                                                                                                          gCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAAt                             <  2:1411165/150‑1 (MQ=255)
                                                                                                                             atcgcatcgTTGTTGGCGATAACCTCTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCt                          >  1:683218/1‑150 (MQ=255)
                                                                                                                                gcatcgTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTa                       <  2:4652405/150‑1 (MQ=255)
                                                                                                                                  atcgTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAAt                     <  1:2341747/150‑1 (MQ=255)
                                                                                                                                  atcgTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAAt                     <  1:4304765/150‑1 (MQ=255)
                                                                                                                                    cgTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCa                   <  1:1018526/150‑1 (MQ=255)
                                                                                                                                    cgTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCa                   >  2:2852085/1‑150 (MQ=255)
                                                                                                                                    cgTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCa                   <  1:1525559/150‑1 (MQ=255)
                                                                                                                                     gttgttGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCAc                  >  2:3125097/1‑150 (MQ=255)
                                                                                                                                     gttgttGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCAc                  >  2:1099260/1‑150 (MQ=255)
                                                                                                                                      ttgttgGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACg                 <  1:3932735/150‑1 (MQ=255)
                                                                                                                                             cGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTgg          >  1:3254614/1‑150 (MQ=255)
                                                                                                                                               aTAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTa        >  2:4074450/1‑150 (MQ=255)
                                                                                                                                               aTAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTa        <  2:428298/150‑1 (MQ=255)
                                                                                                                                                tAACCTCTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTAc       <  1:1631936/150‑1 (MQ=255)
                                                                                                                                                tAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTAc       <  2:1882619/150‑1 (MQ=255)
                                                                                                                                                 aaCCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACg      >  1:3323902/1‑150 (MQ=255)
                                                                                                                                                   ccACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTg    <  2:1074478/150‑1 (MQ=255)
                                                                                                                                                   ccACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTg    <  2:159339/150‑1 (MQ=255)
                                                                                                                                                    cACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGc   <  2:2331225/150‑1 (MQ=255)
                                                                                                                                                     aCTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTATCTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCt  <  1:3827670/150‑1 (MQ=255)
                                                                                                                                                     aCTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCt  <  1:1418227/150‑1 (MQ=255)
                                                                                                                                                     aCTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCt  <  1:4750482/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAGTACGGTGCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTTGTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCT  >  CP009273/2232993‑2233291

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: