Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,602,236 A→T 17.2% Q249H (CAA→CAT lsrA → autoinducer 2 import ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,602,2360AT17.2% 402.3 / 76.2 169Q249H (CAA→CATlsrAautoinducer 2 import ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (72/68);  new base T (15/14);  total (87/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-01

Only 100 of 178 total aligned reads displayed.

CCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTC  >  NC_000913/1602139‑1602335
                                                                                                 |                                                                                                   
ccATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaaa                                                                                                   <  2:1044943/100‑1 (MQ=255)
 cATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaaaa                                                                                                  >  2:3512607/1‑100 (MQ=255)
 cATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaaaa                                                                                                  >  1:2532487/1‑100 (MQ=255)
 cATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaaaa                                                                                                  >  1:162227/1‑100 (MQ=255)
 cATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaaaa                                                                                                  >  2:4005291/1‑100 (MQ=255)
     gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTAt                                                                                              >  1:935020/1‑100 (MQ=255)
      ccTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACGCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTAtn                                                                                             <  1:3292463/100‑2 (MQ=255)
      ccTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCACTCTCTGCCAGCCAAAAATTATg                                                                                             <  2:280039/100‑1 (MQ=255)
        ttAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGc                                                                                           >  2:2488975/1‑100 (MQ=255)
         tAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCt                                                                                          >  2:3845084/1‑100 (MQ=255)
          aaGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTg                                                                                         <  2:2450792/100‑1 (MQ=255)
           aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTgg                                                                                        >  2:4147647/1‑100 (MQ=255)
            gCGGCAAAACCAGCGAACTGTCTACCGCCGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                       <  1:114139/100‑2 (MQ=255)
            gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                       <  2:1569374/100‑2 (MQ=255)
            gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                       <  2:1773529/100‑2 (MQ=255)
            gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                       <  1:1562412/100‑2 (MQ=255)
            gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGa                                                                                       >  2:4247647/1‑100 (MQ=255)
            gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGa                                                                                       <  2:2511866/100‑1 (MQ=255)
             cagCAAAACCAGCGAACTGTCTACCGATGACATTATTCAGGACATCAACAAAAAAAAAAAAAAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                       <  1:1347569/97‑2 (MQ=255)
             cGGCAAAGCGAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAAGCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                                                                      <  2:1461404/100‑1 (MQ=255)
             cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTTGAg                                                                                      >  1:3365855/1‑100 (MQ=255)
             cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                                                                      <  2:2356677/100‑1 (MQ=255)
             cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                                                                      <  2:1754313/100‑1 (MQ=255)
              ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                     <  2:3563252/100‑4 (MQ=255)
              ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                     <  1:3892469/100‑4 (MQ=255)
              ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                     <  1:709313/100‑4 (MQ=255)
               gCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGtt                                                                                    >  1:1830013/1‑100 (MQ=255)
                cAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTa                                                                                   >  2:312942/1‑100 (MQ=255)
                cAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTa                                                                                   >  2:4296355/1‑100 (MQ=255)
                     ccAGCGAACTGTCTACCGACGACATTATTGAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTg                                                                               >  1:449063/1‑99 (MQ=255)
                     ccAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCAAAAAATAAAAAAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTgg                                                                              <  1:752465/100‑1 (MQ=255)
                      cAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGt                                                                             >  2:275101/1‑100 (MQ=255)
                       aGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcagagagc                                                                                <  2:636736/96‑9 (MQ=37)
                       aGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTa                                                                            <  1:1044511/100‑1 (MQ=255)
                        gCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTaa                                                                           <  2:570316/100‑1 (MQ=255)
                         tgAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAAc                                                                          >  1:2249826/2‑100 (MQ=255)
                          gAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAAcc                                                                         >  1:2282315/1‑100 (MQ=255)
                             cTGTCTACCGACGACATTATTCAGGCCATCACCCCTGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGcc                                                                      <  1:3450685/100‑1 (MQ=255)
                             cTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGcc                                                                      <  2:1357560/100‑1 (MQ=255)
                              tGTCTACCGGCAACATAATTCATGCCATCACCCGAGCGGGGCGGGAAAAATCACTCTCTCCCAGCCAAAAATAATGGCTAGAGCTACCTGGTAACCGcc                                                                      >  1:587187/1‑99 (MQ=255)
                               gTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCca                                                                    >  1:637234/1‑100 (MQ=255)
                                 cTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCcaca                                                                  <  1:4267351/100‑1 (MQ=255)
                                 cTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAAGCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCca                                                                    <  2:1143449/98‑1 (MQ=255)
                                  tACCGACGACATTATTGAGGCCATCAGACCAGCGGTAGGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGAAGTTAACCCGTAACCGCCCacaa                                                                 >  1:4082488/1‑100 (MQ=255)
                                  tACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaa                                                                 >  2:3363189/1‑100 (MQ=255)
                                  tACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaa                                                                 >  1:3622064/1‑100 (MQ=255)
                                   aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                                >  2:2548414/1‑100 (MQ=255)
                                   aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                                <  1:835176/100‑1 (MQ=255)
                                   aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                                >  1:4210226/1‑100 (MQ=255)
                                    ccGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaaca                                                               >  2:3533209/1‑100 (MQ=255)
                                     cgacgaCATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACAt                                                              <  2:911309/100‑1 (MQ=255)
                                     cgacgaCATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACAt                                                              <  1:916118/100‑1 (MQ=255)
                                        cgacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgcc                                                           <  1:3261397/100‑1 (MQ=255)
                                         gacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccg                                                          <  2:1862095/100‑1 (MQ=255)
                                          acATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccg                                                          >  1:3203696/1‑99 (MQ=255)
                                           cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccgcc                                                        >  2:318430/1‑100 (MQ=255)
                                           cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccgcc                                                        <  2:4352158/100‑1 (MQ=255)
                                           cATTAGTCGGCCCATCAGCCGAGGGGTACGGGAAAACTCTCTCTGTGCCAGCCAAAACTTGTGGCTGGAGTTGCCTGGTAACCGTCCACAACATgccgcc                                                        <  1:1382432/100‑1 (MQ=255)
                                            attattCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGccgccg                                                       <  1:1029545/100‑1 (MQ=255)
                                              tattCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAATAATTATGGCAGGAGTGACCTGGTAGCCGCCCACAACATGCCGCCGGa                                                     >  1:2969079/1‑100 (MQ=255)
                                               attCAGGCCATCACCCCATCGTTACGGGAAAAATCGCGCTCTACCAACCAAAAATTATGGCTCCAGTTCCCCCGTAACCGACCACAACCTGccgccg                                                       <  1:2128903/97‑1 (MQ=255)
                                               attCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGaa                                                    >  1:1101226/1‑100 (MQ=255)
                                                 tCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACg                                                  <  2:3757064/100‑1 (MQ=255)
                                                  cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGc                                                 >  2:4210048/1‑100 (MQ=255)
                                                   aGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGcc                                                >  1:2911866/1‑100 (MQ=255)
                                                    ggCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCg                                               >  2:2603501/1‑100 (MQ=255)
                                                     gCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCgg                                              <  2:2072953/100‑1 (MQ=255)
                                                     gCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCgg                                              <  2:763696/100‑1 (MQ=255)
                                                            ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGac                                       >  1:4002029/1‑100 (MQ=255)
                                                                   ggTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaa                                <  1:1529104/100‑1 (MQ=255)
                                                                   ggTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaa                                >  1:352782/1‑100 (MQ=255)
                                                                       cGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                            >  2:3991323/1‑100 (MQ=255)
                                                                       cGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                            <  1:838509/100‑1 (MQ=255)
                                                                        gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                           >  2:2914737/1‑100 (MQ=255)
                                                                        gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                           <  1:2267198/100‑1 (MQ=255)
                                                                        gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                           >  1:3787752/1‑100 (MQ=255)
                                                                         ggAAAAATCGTTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTg                          >  1:4354427/1‑100 (MQ=255)
                                                                         ggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTg                          <  2:1709827/100‑1 (MQ=255)
                                                                            aaaaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGAcc                       <  1:794763/100‑1 (MQ=255)
                                                                                tCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCg                   <  2:759947/100‑1 (MQ=255)
                                                                                  gCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGaa                 <  1:2738126/100‑1 (MQ=255)
                                                                                   ctctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAg                <  2:3066525/100‑1 (MQ=255)
                                                                                    tctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAgg               <  2:3460355/100‑1 (MQ=255)
                                                                                       nctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGttt            >  2:686140/4‑100 (MQ=255)
                                                                                        ctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTc           >  2:3381056/3‑100 (MQ=255)
                                                                                        ctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTc           >  2:2910619/3‑100 (MQ=255)
                                                                                        ctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTc           >  1:373225/3‑100 (MQ=255)
                                                                                         tccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa          >  1:29942/2‑100 (MQ=255)
                                                                                         tccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa          >  2:1924269/2‑100 (MQ=255)
                                                                                         gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa          >  2:2933674/1‑100 (MQ=255)
                                                                                         gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa          >  2:3998154/1‑100 (MQ=255)
                                                                                          ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg         >  2:683257/1‑100 (MQ=255)
                                                                                          ccagccaAAAATTATGGCTGGAGTTACCTGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAgg        >  1:3772812/1‑100 (MQ=255)
                                                                                             gccaAAAATTATGGCTGGAGTTACCTGGCAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg         <  1:2835243/97‑1 (MQ=255)
                                                                                               naaaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg    >  2:2030963/2‑100 (MQ=255)
                                                                                               caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg    >  1:3429168/1‑100 (MQ=255)
                                                                                               caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg    >  2:2702740/1‑100 (MQ=255)
                                                                                               caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg    >  2:161464/1‑100 (MQ=255)
                                                                                               caAAAATTATGGCTGGAGTCACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg    >  1:1420054/1‑100 (MQ=255)
                                                                                                 aaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACGTGAAAATCTGACCGGCGAAGGTTTCAGGAATGTc  >  1:1966053/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
CCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTC  >  NC_000913/1602139‑1602335

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