Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,904,858 G→A 24.1% T247I (ACA→ATA)  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
*NC_0009133,904,8580GA24.1% 324.9 / 98.4 178T247I (ACA→ATA) bglFfused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base G (72/63);  new base A (24/19);  total (96/82)
Fisher's exact test for biased strand distribution p-value = 8.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

Only 100 of 185 total aligned reads displayed.

CAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACG  >  NC_000913/3904759‑3904957
                                                                                                   |                                                                                                   
cAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgtg                                                                                                     <  2:3825831/100‑1 (MQ=255)
cAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgtg                                                                                                     <  1:2198068/100‑1 (MQ=255)
   agaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTga                                                                                                  <  1:3558227/100‑1 (MQ=255)
    gaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaa                                                                                                 >  1:282070/1‑100 (MQ=255)
     aTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGTTGTGATAACCATCAGACATAGCAATGGTGTgaag                                                                                                <  2:601144/100‑1 (MQ=255)
      tAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaaga                                                                                               <  1:3253893/100‑1 (MQ=255)
      tAACCGGCGGCAATCAGTTCGCTTATACAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaaga                                                                                               >  2:1598827/1‑100 (MQ=255)
       aaCCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaagaa                                                                                              >  1:439884/1‑100 (MQ=255)
        aCCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGaaa                                                                                             >  1:280136/1‑100 (MQ=255)
        aCCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTATGAAGaaa                                                                                             >  1:1300014/1‑100 (MQ=255)
            gcggcAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAtttt                                                                                         >  1:3536493/1‑100 (MQ=255)
             cgggAAGCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAttttt                                                                                        <  2:3087421/100‑1 (MQ=255)
               gcAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGa                                                                                      <  2:1117630/100‑1 (MQ=255)
               gcAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTATGAAGAAATTTTTGa                                                                                      <  2:2178485/100‑1 (MQ=255)
                cAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGAt                                                                                     >  2:2489157/1‑100 (MQ=255)
                cAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGAt                                                                                     >  1:223151/1‑100 (MQ=255)
                cAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGAt                                                                                     >  1:341063/1‑100 (MQ=255)
                 aaTCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAAACATCAGACATAGCAATGGTATGAAGAAATTTTTGAtt                                                                                    >  1:957507/1‑100 (MQ=255)
                  aTCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTg                                                                                   <  2:3536454/100‑1 (MQ=255)
                   tCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTATGAAGAAATTTTTGATTGc                                                                                  <  1:138855/100‑1 (MQ=255)
                    cAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGcc                                                                                 >  2:1026790/1‑100 (MQ=255)
                    cAGTTCGCTTATCCAGGTTGATAGCGGACCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCCACGGTGTGAAGAAATTTTTGATTAcc                                                                                 >  2:3173294/1‑100 (MQ=255)
                      gTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGa                                                                               >  1:2101815/1‑100 (MQ=255)
                       ttCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGAc                                                                              <  2:1408593/100‑1 (MQ=255)
                       ttCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTATGAAGAAATTTTTGATTGCCGAc                                                                              <  2:2396972/100‑1 (MQ=255)
                        tCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACg                                                                             >  2:698357/1‑100 (MQ=255)
                        tCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACg                                                                             <  2:1455738/100‑1 (MQ=255)
                         cGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTATGAAGAAATTTTTGATTGCCGACgg                                                                            >  1:1874528/1‑100 (MQ=255)
                          gCTTATCCAGGTTGATAGCGGCCTCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTATGAAGAAATTTTTGATTGCCGACgg                                                                            <  1:2449326/99‑1 (MQ=255)
                           cTTATCCAGGTTGTTAGCGGCGGCACCAGCAGAAAGGTGACGGGTGTGATAAGTGTCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGt                                                                           >  2:1950203/1‑99 (MQ=255)
                           cTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTa                                                                          >  1:1819230/1‑100 (MQ=255)
                           cTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTa                                                                          >  1:3098705/1‑100 (MQ=255)
                             tATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAAc                                                                        <  2:163866/100‑1 (MQ=255)
                               tCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCa                                                                      >  1:3163013/1‑100 (MQ=255)
                                ccAGGTTGATAGCGGCCCCACCAGCCGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCAc                                                                     >  1:2548440/1‑100 (MQ=255)
                                ccAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCAc                                                                     >  2:1871455/1‑100 (MQ=255)
                                ccAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCAc                                                                     >  2:480246/1‑100 (MQ=255)
                                ccAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCAc                                                                     >  2:2347155/1‑100 (MQ=255)
                                 cAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACg                                                                    >  2:3261695/1‑100 (MQ=255)
                                      tGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGTCGACGGTAACCACGCATTa                                                               <  2:630457/100‑1 (MQ=255)
                                      tGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTa                                                               <  1:661546/100‑1 (MQ=255)
                                      tGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTATGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTa                                                               >  2:1902023/1‑100 (MQ=255)
                                          aGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTc                                                           >  1:2283986/1‑100 (MQ=255)
                                            cGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCgg                                                         >  2:1996096/1‑100 (MQ=255)
                                             ggCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGATATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGc                                                        <  1:480253/100‑1 (MQ=255)
                                              gCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCg                                                       <  2:212913/100‑1 (MQ=255)
                                               ccccACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGAAATAGCAATGGTATGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGt                                                      <  1:2198443/100‑1 (MQ=255)
                                                 ccaccaGCAGAAATGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGtt                                                     <  1:750289/99‑1 (MQ=255)
                                                  caccaGCAGAAAGGTGACGGGTGTGATAAGCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTcc                                                   <  2:1141464/100‑1 (MQ=255)
                                                  caccaGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTcc                                                   <  2:2457781/100‑1 (MQ=255)
                                                     cagcagAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGa                                                >  2:2297413/1‑100 (MQ=255)
                                                       gcagAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAAt                                              <  1:698365/100‑1 (MQ=255)
                                                       gcagAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACAGTAACCACGCATTAAGTCGGCGTTCCAGAAt                                              >  1:1463459/1‑100 (MQ=255)
                                                          gAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTATGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAAtgct                                           <  2:766766/100‑1 (MQ=255)
                                                            aaGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTATGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAAtgctgc                                         <  2:1127850/100‑1 (MQ=255)
                                                            aaGGCGACGGGTGTGATAACCATCAGCCATAGCAATGGTGTGAAGAAATTTCTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAAtgctgc                                         <  2:3737602/100‑1 (MQ=255)
                                                             aGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGca                                        <  1:1983308/100‑1 (MQ=255)
                                                             aGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGca                                        <  2:809712/100‑1 (MQ=255)
                                                             aGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGca                                        >  2:2894638/1‑100 (MQ=255)
                                                             aGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTATGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGca                                        >  2:3720776/1‑100 (MQ=255)
                                                              ggTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGca                                        >  2:1748084/1‑99 (MQ=255)
                                                              ggTGACGGGTGTGATAACCATCAGACATAGCAATGGTATGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGcac                                       >  1:2361394/1‑100 (MQ=255)
                                                               gTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGcaca                                      <  2:1244788/100‑1 (MQ=255)
                                                                tGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACaa                                     <  2:2062142/100‑1 (MQ=255)
                                                                 gACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAAc                                    >  2:659537/1‑100 (MQ=255)
                                                                 gACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAAc                                    <  1:2696415/100‑1 (MQ=255)
                                                                  aCGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAAcc                                   <  2:2624221/100‑1 (MQ=255)
                                                                     ggTGTGACAACAATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAgg                                >  2:2054831/1‑100 (MQ=255)
                                                                      gtgtGATAACCATCAGACATAGTAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGc                               <  1:813394/100‑1 (MQ=255)
                                                                       tgtgATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAgg                                <  1:213062/98‑1 (MQ=255)
                                                                         tgATAACCCTCAGACATAGCAATGGTATGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTACAGAATGCTGCACAACCAGGCAGa                            >  2:733215/1‑100 (MQ=255)
                                                                         tgATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTTTTTGTCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTTCACAACCAGGCAGa                            >  2:3217355/1‑100 (MQ=255)
                                                                         tgATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGa                            >  2:3688708/1‑100 (MQ=255)
                                                                           aTAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGaaa                          >  1:115764/1‑100 (MQ=255)
                                                                               ccATCAGACATAGCAATGGTGTGAAGAAATTTTTTATTTCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAaat                      >  2:3796565/1‑100 (MQ=255)
                                                                               ccATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACAACCCATCAAGTCGGAGTTCCAGAATGCTGCACAACCAGGCAGaaaaa                        >  2:2690075/1‑98 (MQ=255)
                                                                               ccAGCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAaat                      >  2:386271/1‑100 (MQ=255)
                                                                                cATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCAGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAaata                     >  1:2461403/1‑100 (MQ=255)
                                                                                cATCAGACATAGCAATGGTATGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAaata                     >  2:828173/1‑100 (MQ=255)
                                                                                 aTCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAaataa                    >  1:3547099/1‑100 (MQ=255)
                                                                                  tCAGACATAGTAATGGTGTGAAGAAATTTTTGATTGTCGTCGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAaataat                   <  1:2217529/100‑1 (MQ=255)
                                                                                    aGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATgg                 >  1:743852/1‑100 (MQ=255)
                                                                                    aGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATgg                 <  1:3860926/100‑1 (MQ=255)
                                                                                    aGACATAGCAATGGTATGAAGAAATTTTTGATTGTCGATGTTAACCACGTATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATgg                 <  1:184033/100‑1 (MQ=255)
                                                                                     gACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATggg                >  1:2711854/1‑100 (MQ=255)
                                                                                     gACATAGCAATGGTATGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATggg                >  1:1684038/1‑100 (MQ=255)
                                                                                     gACATAGCAATGGTATGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATggg                <  2:2517255/100‑1 (MQ=255)
                                                                                       cATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGaa              >  2:3300357/1‑100 (MQ=255)
                                                                                           gCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAAc          >  1:3121522/1‑100 (MQ=255)
                                                                                           gCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAAc          >  1:1734367/1‑100 (MQ=255)
                                                                                             aaTTGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCAAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAAc          <  2:2689953/98‑1 (MQ=255)
                                                                                             aaTGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAAAAAAAATAATGGGAATAAcc         <  2:3135642/99‑1 (MQ=255)
                                                                                              aTGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGa       <  1:2862790/100‑1 (MQ=255)
                                                                                               tGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAAc          <  2:3426523/96‑1 (MQ=255)
                                                                                               tGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGAt      <  2:2433771/100‑1 (MQ=255)
                                                                                               tGGTATGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCTGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGAt      <  1:2039004/100‑1 (MQ=255)
                                                                                                ggTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATg     <  1:660450/100‑1 (MQ=255)
                                                                                                ggTATGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATg     >  1:248406/1‑100 (MQ=255)
                                                                                                  tgtgAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGAc   >  2:2557663/1‑100 (MQ=255)
                                                                                                   gtgAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACg  <  2:3015856/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACG  >  NC_000913/3904759‑3904957

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: