Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,897,492 G→T 8.4% P489Q (CCG→CAG)  bglF ← fused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,897,4920GT 8.4% 1630.1 / 86.7 427P489Q (CCG→CAG) bglFfused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base G (168/223);  new base T (18/18);  total (186/241)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 437 total aligned reads displayed.

TGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCACCGTACCGACAAATGCGCAACC  >  CP009273/3897344‑3897640
                                                                                                                                                    |                                                                                                                                                    
 gACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCgg                                                                                                                                                     >  2:733966/1‑149 (MQ=255)
 gACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCgg                                                                                                                                                     <  2:5067013/149‑1 (MQ=255)
     cAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGcaca                                                                                                                                                 >  2:1726848/1‑149 (MQ=255)
      aaTGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACaa                                                                                                                                                <  2:2076861/149‑1 (MQ=255)
       aTGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTCTCAGCGACGTGAATGAGCGGCACAAt                                                                                                                                               >  1:1296021/1‑149 (MQ=255)
       aTGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAAt                                                                                                                                               >  1:440315/1‑149 (MQ=255)
       aTGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAAt                                                                                                                                               <  1:2091565/149‑1 (MQ=255)
        tGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAAtc                                                                                                                                              >  2:4637444/1‑149 (MQ=255)
                 cGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTc                                                                                                                                     <  1:3222094/149‑1 (MQ=255)
                   tgtAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCAt                                                                                                                                   >  2:564655/1‑149 (MQ=255)
                    gtAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCACCCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCTGCACAAGCTCTCCCGTCATc                                                                                                                                  <  2:4685299/149‑1 (MQ=255)
                    gtAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATc                                                                                                                                  >  1:2937482/1‑149 (MQ=255)
                           ggCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAAc                                                                                                                           <  1:1390473/149‑1 (MQ=255)
                            gCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACa                                                                                                                          <  1:4136434/149‑1 (MQ=255)
                              gAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGa                                                                                                                        >  1:3391792/1‑149 (MQ=255)
                                acaacGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATa                                                                                                                      >  1:60547/1‑149 (MQ=255)
                                      aaGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGccc                                                                                                                <  2:4349483/149‑1 (MQ=255)
                                      aaGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGccc                                                                                                                <  1:724476/149‑1 (MQ=255)
                                      aaGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGccc                                                                                                                <  2:604705/149‑1 (MQ=255)
                                       aGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGTCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCt                                                                                                               >  2:2123047/1‑149 (MQ=255)
                                       aGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCt                                                                                                               >  1:3177901/1‑149 (MQ=255)
                                         cAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGc                                                                                                             >  1:2453169/1‑149 (MQ=255)
                                          aaTTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCTGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCt                                                                                                            >  1:380245/1‑149 (MQ=255)
                                                aCCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAggtg                                                                                                      <  2:272084/149‑1 (MQ=255)
                                                aCCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAggtg                                                                                                      <  2:2015902/149‑1 (MQ=255)
                                                 cccGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAggtgg                                                                                                     <  2:343343/149‑1 (MQ=255)
                                                  ccGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAggtggt                                                                                                    <  1:4900122/149‑1 (MQ=255)
                                                  ccGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAggtggt                                                                                                    <  1:1980609/149‑1 (MQ=255)
                                                   cGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGgtggtg                                                                                                   <  2:1427098/149‑1 (MQ=255)
                                                     cAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGta                                                                                                 <  1:3940707/149‑1 (MQ=255)
                                                      aaCCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCTGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGtaa                                                                                                >  2:4997985/1‑149 (MQ=255)
                                                         cGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGtaataa                                                                                             <  2:3789305/149‑1 (MQ=255)
                                                          ggAGAACGCACGTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCACCAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCTGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAAc                                                                                            <  2:3396339/149‑1 (MQ=255)
                                                           gagaACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAAcc                                                                                           <  2:2386007/149‑1 (MQ=255)
                                                           gagaACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAAcc                                                                                           <  2:4018071/149‑1 (MQ=255)
                                                             gaACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACctc                                                                                         <  1:1252732/149‑1 (MQ=255)
                                                                   aCTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGgtgt                                                                                   >  2:4063121/1‑149 (MQ=255)
                                                                   aCTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGgtgt                                                                                   <  1:4195321/149‑1 (MQ=255)
                                                                   aCTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGgtgt                                                                                   <  2:1904238/149‑1 (MQ=255)
                                                                   aCTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGgtgt                                                                                   >  2:2214534/1‑149 (MQ=255)
                                                                       cACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTtc                                                                               <  2:1852173/149‑1 (MQ=255)
                                                                       cACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTtc                                                                               <  1:4033473/149‑1 (MQ=255)
                                                                        aCCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTtct                                                                              >  1:981934/1‑149 (MQ=255)
                                                                          cAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTctct                                                                            <  2:910504/149‑1 (MQ=255)
                                                                           aaCCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCtt                                                                           <  2:147918/149‑1 (MQ=255)
                                                                           aaCCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCtt                                                                           >  1:4234103/1‑149 (MQ=255)
                                                                            aCCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTg                                                                          <  2:295855/149‑1 (MQ=255)
                                                                            aCCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTg                                                                          <  1:4709313/149‑1 (MQ=255)
                                                                            aCCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTg                                                                          >  1:87126/1‑149 (MQ=255)
                                                                               gAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCgg                                                                       <  1:3034305/149‑1 (MQ=255)
                                                                                aGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCggg                                                                      <  2:488381/149‑1 (MQ=255)
                                                                                aGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCggg                                                                      >  1:266480/1‑149 (MQ=255)
                                                                                  ggCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGc                                                                    <  1:1002261/149‑1 (MQ=255)
                                                                                      gAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCAccc                                                                >  1:2328038/1‑149 (MQ=255)
                                                                                      gAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCAccc                                                                >  1:286319/1‑149 (MQ=255)
                                                                                          ggCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTgcg                                                            >  1:1723900/1‑149 (MQ=255)
                                                                                          ggCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTgcg                                                            >  1:3754731/1‑149 (MQ=255)
                                                                                             aaTACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCgctg                                                         >  2:1059354/1‑149 (MQ=255)
                                                                                               tACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCgctggc                                                       >  1:207966/1‑149 (MQ=255)
                                                                                                  cTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGa                                                    >  2:859975/1‑149 (MQ=255)
                                                                                                  cTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGa                                                    >  1:4852451/1‑149 (MQ=255)
                                                                                                   tttACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGAc                                                   <  2:692591/149‑1 (MQ=255)
                                                                                                   tttACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGAc                                                   >  2:2992431/1‑149 (MQ=255)
                                                                                                      aCCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGtt                                                <  1:2616676/149‑1 (MQ=255)
                                                                                                      aCCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGtt                                                >  1:3624980/1‑149 (MQ=255)
                                                                                                       cccAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCTGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGttt                                               <  2:1168285/149‑1 (MQ=255)
                                                                                                         cAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAg                                             <  2:4499923/149‑1 (MQ=255)
                                                                                                         cAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAg                                             <  2:3630208/149‑1 (MQ=255)
                                                                                                          aaCAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGc                                            >  1:626730/1‑149 (MQ=255)
                                                                                                          aaCAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGc                                            >  1:4591130/1‑149 (MQ=255)
                                                                                                            cAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCgg                                          >  1:3107650/1‑149 (MQ=255)
                                                                                                               gCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGtga                                       <  2:2616104/149‑1 (MQ=255)
                                                                                                               gCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGtga                                       <  2:3502308/149‑1 (MQ=255)
                                                                                                               gCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGtga                                       <  1:694428/149‑1 (MQ=255)
                                                                                                                 cACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGtgatg                                     <  1:1475280/149‑1 (MQ=255)
                                                                                                                  aCTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGtgatga                                    <  2:3289318/149‑1 (MQ=255)
                                                                                                                  aCTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGtgatga                                    <  2:4217009/149‑1 (MQ=255)
                                                                                                                   cTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGatgaa                                   <  1:5117007/149‑1 (MQ=255)
                                                                                                                    tGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGatgaaa                                  <  1:2726982/149‑1 (MQ=255)
                                                                                                                     ggCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGatgaaat                                 >  1:1778381/1‑149 (MQ=255)
                                                                                                                      gCAAACGTGGTATCAGCGACGTGAATGAGCTGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGatgaaatg                                >  1:3410128/1‑149 (MQ=255)
                                                                                                                       cAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGatgaaatga                               <  2:2328038/149‑1 (MQ=255)
                                                                                                                        aaaCGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGatgaaatgaa                              <  2:821666/149‑1 (MQ=255)
                                                                                                                         aaCGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAg                             >  1:1666420/1‑149 (MQ=255)
                                                                                                                         aaCGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAg                             >  1:2368528/1‑149 (MQ=255)
                                                                                                                         aaCGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAg                             >  1:255349/1‑149 (MQ=255)
                                                                                                                          aCGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGc                            >  1:4542009/1‑149 (MQ=255)
                                                                                                                             tGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATc                         >  1:4774276/1‑149 (MQ=255)
                                                                                                                              ggTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCa                        <  2:4682089/149‑1 (MQ=255)
                                                                                                                              ggTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCa                        >  1:557454/1‑149 (MQ=255)
                                                                                                                               gTATCAGCGACGTGAATGAGCTGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCAc                       <  2:380245/149‑1 (MQ=255)
                                                                                                                                    aGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCACCGTAc                  <  2:3232594/149‑1 (MQ=255)
                                                                                                                                    aGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCACCGTAc                  <  1:330109/149‑1 (MQ=255)
                                                                                                                                        aCGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCACCGTACCGAc              <  1:3599526/149‑1 (MQ=255)
                                                                                                                                         cGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCACCGTACCGACa             <  1:2226989/149‑1 (MQ=255)
                                                                                                                                            gAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCACCGTACCGACAAAt          <  1:115278/149‑1 (MQ=255)
                                                                                                                                             aaTGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCACCGTACCGACAAATg         >  1:1564832/1‑149 (MQ=255)
                                                                                                                                                   cGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATATCCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCACCGTACCGACAAATGCGCAAc   <  1:2295922/149‑1 (MQ=255)
                                                                                                                                                   cGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCACCGTACCGACAAATGCGCAAc   <  2:305918/149‑1 (MQ=255)
                                                                                                                                                   cGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCACCGTACCGACAAATGCGCAAc   <  1:3256198/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGACTCAATGCCAATGGCGTGTAATGTGGCGAACAACGAAGCAATTCGACCCGCAACCGGAGAACGCACTTCACCAACCGAGGGCAGAATGGCAATACCTTTACCCAACAGGCCACTGGCAAACGTGGTATCAGCGACGTGAATGAGCGGCACAATCTCTCCCGTCATCGGTGAACAGATACCGCCCTGCTCAGGTGGTGTAATAACCTCTGGTGTTTTCTCTTGCGGGGCACCCTGCGCTGGCTGACGTTTAGCGGTGATGAAATGAAGCATCACCGTACCGACAAATGCGCAACC  >  CP009273/3897344‑3897640

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: