Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,902,432 G→A P265L (CCG→CTG)  bglF ← PTS system beta‑glucosides‑specific enzyme IIABC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,902,4320GA99.3% 1533.6 / ‑3.9 421P265L (CCG→CTG) bglFPTS system beta‑glucosides‑specific enzyme IIABC component
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (250/168);  total (251/170)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 440 total aligned reads displayed.

CCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCA  >  NC_000913/3902342‑3902521
                                                                                          |                                                                                         
ccgccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATagcagc                                                                                         >  1:1189644/1‑93 (MQ=255)
ccgccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATagcagc                                                                                         <  2:2678827/93‑1 (MQ=255)
 cgccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGcc                                                                                        <  2:7578659/93‑1 (MQ=255)
  gccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGccc                                                                                       >  1:10040738/1‑93 (MQ=255)
  gccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGccc                                                                                       >  2:949999/1‑93 (MQ=255)
   cccATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGcccc                                                                                      >  2:9798660/1‑93 (MQ=255)
     cATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccac                                                                                    >  1:8495881/1‑93 (MQ=255)
     cATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccac                                                                                    <  2:9123845/93‑1 (MQ=255)
     cATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccac                                                                                    <  1:5825142/93‑1 (MQ=255)
     cATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccac                                                                                    >  1:18728/1‑93 (MQ=255)
      aTTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccacc                                                                                   >  2:8086712/1‑93 (MQ=255)
       ttACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccacca                                                                                  >  2:1592308/1‑93 (MQ=255)
       ttACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccacca                                                                                  >  2:7196266/1‑93 (MQ=255)
       ttACCGCGCCCGCAAATCCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCccacca                                                                                  >  2:8402806/1‑93 (MQ=255)
        tACCGCGCCAGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACcag                                                                                 >  1:4546030/1‑93 (MQ=255)
            gcgcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGa                                                                             <  1:7340439/93‑1 (MQ=255)
            gcgcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGa                                                                             >  1:1337262/1‑93 (MQ=255)
             cgcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGaa                                                                            <  1:3403300/93‑1 (MQ=255)
             cgcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGaa                                                                            >  1:2930859/1‑93 (MQ=255)
              gcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGaaa                                                                           <  2:3813770/93‑1 (MQ=255)
              gcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGaaa                                                                           <  1:5300234/93‑1 (MQ=255)
               cccGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAg                                                                          <  2:5531910/93‑1 (MQ=255)
               cccGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAg                                                                          >  1:9660239/1‑93 (MQ=255)
                ccGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAgg                                                                         >  1:8886947/1‑93 (MQ=255)
                 cGCAGATGCAGGATCCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGt                                                                        <  1:7782604/93‑1 (MQ=255)
                 cGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGt                                                                        <  2:4094323/93‑1 (MQ=255)
                  gCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTg                                                                       >  1:6499679/1‑93 (MQ=255)
                   cAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGa                                                                      <  2:1245295/93‑1 (MQ=255)
                    aaaTGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGAc                                                                     <  1:9781280/93‑1 (MQ=255)
                    aaaTGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGAc                                                                     <  1:7589720/93‑1 (MQ=255)
                      aTGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACgg                                                                   <  1:7447588/93‑1 (MQ=255)
                       tGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACggg                                                                  >  1:5010354/1‑93 (MQ=255)
                       tGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACggg                                                                  >  2:1771346/1‑93 (MQ=255)
                          aGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGgtgt                                                               <  2:8901018/93‑1 (MQ=255)
                          aGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGgtgt                                                               <  1:2223784/93‑1 (MQ=255)
                           ggAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGtgtg                                                              >  1:995871/1‑93 (MQ=255)
                           ggAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGtgtg                                                              >  1:3878079/1‑93 (MQ=255)
                            gAACCGGCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGa                                                             <  1:793638/93‑1 (MQ=255)
                               ccGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATaa                                                          >  2:6050259/1‑93 (MQ=255)
                                  ccTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCa                                                       <  2:4905630/93‑1 (MQ=255)
                                  ccTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCa                                                       >  2:8532313/1‑93 (MQ=255)
                                    tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     >  1:6776832/1‑93 (MQ=255)
                                    tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     >  1:164025/1‑93 (MQ=255)
                                    tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     <  2:2868816/93‑1 (MQ=255)
                                    tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     >  1:8673015/1‑93 (MQ=255)
                                    tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     >  1:241893/1‑93 (MQ=255)
                                    tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATc                                                     >  2:6298502/1‑93 (MQ=255)
                                     gATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCa                                                    >  2:3558054/1‑93 (MQ=255)
                                      aTAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAg                                                   >  1:4219073/1‑93 (MQ=255)
                                       tAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGa                                                  <  2:7430293/93‑1 (MQ=255)
                                         aaGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACa                                                <  1:1922333/93‑1 (MQ=255)
                                          aGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACAt                                               <  2:612075/93‑1 (MQ=255)
                                           gCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATa                                              >  1:3565905/1‑93 (MQ=255)
                                           accAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATa                                              <  1:6011352/92‑1 (MQ=255)
                                            ccAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAg                                             <  1:6811154/93‑1 (MQ=255)
                                             cAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGc                                            >  2:9057225/1‑93 (MQ=255)
                                             cAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGc                                            >  2:5681446/1‑93 (MQ=255)
                                             cAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGc                                            <  1:6254357/93‑1 (MQ=255)
                                              agagATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCa                                           >  2:2391753/1‑93 (MQ=255)
                                              agagATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCa                                           <  1:6659211/93‑1 (MQ=255)
                                               gagaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCaa                                          <  2:6962400/93‑1 (MQ=255)
                                               gagaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCaa                                          >  1:2523828/1‑93 (MQ=255)
                                                agaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCaag                                         >  1:2006825/1‑92 (MQ=255)
                                                 gaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATg                                        >  2:9210443/1‑93 (MQ=255)
                                                 gaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATg                                        <  2:5336342/93‑1 (MQ=255)
                                                  aTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATgg                                       >  2:3342309/1‑93 (MQ=255)
                                                   tAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgt                                      >  2:6586198/1‑93 (MQ=255)
                                                   tAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgt                                      >  1:9358899/1‑93 (MQ=255)
                                                       cggcggCAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAACGGTGTga                                  <  1:308651/93‑1 (MQ=255)
                                                         gcggcAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaag                                >  2:6736539/1‑93 (MQ=255)
                                                          cggcAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaaga                               <  1:4473965/93‑1 (MQ=255)
                                                          cggcAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaaga                               <  2:2488386/93‑1 (MQ=255)
                                                           ggcAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaagaa                              >  1:3682081/1‑93 (MQ=255)
                                                            gcAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGaaa                             >  1:3936591/1‑93 (MQ=255)
                                                            gcAATCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGaaa                             <  1:5703485/93‑1 (MQ=255)
                                                              aaTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAtt                           >  1:5099568/1‑93 (MQ=255)
                                                               aTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAttt                          >  2:2044192/1‑93 (MQ=255)
                                                               aTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAttt                          >  1:9411993/1‑93 (MQ=255)
                                                               aTCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAttt                          >  2:3060655/1‑93 (MQ=255)
                                                                tCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAtttt                         >  2:5043084/1‑93 (MQ=255)
                                                                tCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAtttt                         >  2:3988729/1‑93 (MQ=255)
                                                                tCAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAtttt                         >  1:4078303/1‑93 (MQ=255)
                                                                 cAGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAttttt                        <  2:2277506/93‑1 (MQ=255)
                                                                  aGTTCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTg                       >  2:3220708/1‑93 (MQ=255)
                                                                    ttCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCANTGGTGTGAAGAAATTTTTGAt                     >  1:7144861/1‑93 (MQ=255)
                                                                    ttCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGAt                     <  2:7027551/93‑1 (MQ=255)
                                                                     tCGCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGAtt                    <  1:8170114/93‑1 (MQ=255)
                                                                       gCTTATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGc                  >  1:8345424/1‑93 (MQ=255)
                                                                          tATCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGa               >  1:5520288/1‑93 (MQ=255)
                                                                           aTCCAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGAc              >  2:4526359/1‑93 (MQ=255)
                                                                             ccAGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACgg            <  2:1157104/93‑1 (MQ=255)
                                                                               aGGTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTa          <  1:47955/93‑1 (MQ=255)
                                                                                ggTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCATTGGTGTGAAGAAATTTTTGATTGCCGACGGTac         >  2:8113954/1‑92 (MQ=255)
                                                                                ggTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTaa         <  2:6510743/93‑1 (MQ=255)
                                                                                ggTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTaa         <  2:2111017/93‑1 (MQ=255)
                                                                                 gTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAAc        >  1:8786239/1‑93 (MQ=255)
                                                                                 gTTGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAAc        >  1:6965036/1‑93 (MQ=255)
                                                                                  ttGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAAcc       <  1:2881950/93‑1 (MQ=255)
                                                                                  ttGATAGCAGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAAcc       >  2:1964872/1‑93 (MQ=255)
                                                                                       agcagcCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCa  <  2:7943124/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCA  >  NC_000913/3902342‑3902521

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