Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,602,239 A→T 17.3% K250N (AAA→AAT 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,2390AT17.3% 399.2 / 75.4 169K250N (AAA→AATlsrAautoinducer 2 import ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (72/66);  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 = 2.21e-01

Only 100 of 178 total aligned reads displayed.

CATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAG  >  NC_000913/1602140‑1602337
                                                                                                   |                                                                                                  
cATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaaaa                                                                                                    >  2:3512607/1‑100 (MQ=255)
cATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaaaa                                                                                                    >  1:162227/1‑100 (MQ=255)
cATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaaaa                                                                                                    >  2:4005291/1‑100 (MQ=255)
    gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTAt                                                                                                >  1:2409120/1‑100 (MQ=255)
     ccTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATg                                                                                               <  1:3273271/100‑1 (MQ=255)
     ccTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCACTCTCTGCCAGCCAAAAATTATg                                                                                               <  2:280039/100‑1 (MQ=255)
      cTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATg                                                                                               >  1:4064959/1‑99 (MQ=255)
       ttAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGc                                                                                             >  2:2488975/1‑100 (MQ=255)
        tAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCt                                                                                            >  1:921328/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)
          aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTgg                                                                                          >  1:763694/1‑100 (MQ=255)
           gCGGCAAAACCAGCGAACTGTCTACCGCCGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                         <  1:114139/100‑2 (MQ=255)
           gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                         <  2:1773529/100‑2 (MQ=255)
           gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                         <  2:1569374/100‑2 (MQ=255)
           gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                         <  1:4351157/100‑2 (MQ=255)
           gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                         <  1:445631/100‑2 (MQ=255)
           gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGa                                                                                         <  2:2511866/100‑1 (MQ=255)
           gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGa                                                                                         >  2:4247647/1‑100 (MQ=255)
            cagCAAAACCAGCGAACTGTCTACCGATGACATTATTCAGGACATCAACAAAAAAAAAAAAAAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                         <  1:1347569/97‑2 (MQ=255)
            cGGCAAAGCGAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAAGCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                                                                        <  2:1461404/100‑1 (MQ=255)
            cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                                                                        <  2:2356677/100‑1 (MQ=255)
            cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                                                                        <  2:1754313/100‑1 (MQ=255)
             ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                       <  2:3563252/100‑4 (MQ=255)
               cAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTa                                                                                     >  2:312942/1‑100 (MQ=255)
               cAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTa                                                                                     >  2:4296355/1‑100 (MQ=255)
                 aaaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                                                                   <  1:693475/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)
                            cTGTCTACCGACGACATTATTCAGGCCATCACCCCTGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGcc                                                                        <  1:3450685/100‑1 (MQ=255)
                            cTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGcc                                                                        <  2:1357560/100‑1 (MQ=255)
                             tGTCTACCGGCAACATAATTCATGCCATCACCCGAGCGGGGCGGGAAAAATCACTCTCTCCCAGCCAAAAATAATGGCTAGAGCTACCTGGTAACCGcc                                                                        >  1:587187/1‑99 (MQ=255)
                             tGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGccc                                                                       <  1:3445367/100‑1 (MQ=255)
                               tCTACCGACGACATTATTCAGGCCATCTCCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCcac                                                                     <  1:150417/100‑1 (MQ=255)
                               tCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAACTTATGGCTGGAGTTACCTGGTAACCGCCcac                                                                     <  1:462748/100‑1 (MQ=255)
                                cTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAAGCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCca                                                                      <  2:1143449/98‑1 (MQ=255)
                                 tACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaa                                                                   >  1:1576973/1‑100 (MQ=255)
                                 tACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaa                                                                   >  2:3363189/1‑100 (MQ=255)
                                  aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                                  >  1:4210226/1‑100 (MQ=255)
                                  aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                                  <  1:835176/100‑1 (MQ=255)
                                  aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                                  >  2:2548414/1‑100 (MQ=255)
                                   ccGACGTCATTATTCAGGCCGTCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                                  <  1:835142/99‑1 (MQ=255)
                                   ccGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaaca                                                                 >  2:3533209/1‑100 (MQ=255)
                                    cgacgaCATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACAt                                                                <  2:911309/100‑1 (MQ=255)
                                     gacgacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATg                                                               >  1:3691496/1‑100 (MQ=255)
                                     gacgacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTCCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATg                                                               >  1:4206068/1‑100 (MQ=255)
                                        gacATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccg                                                            <  2:1862095/100‑1 (MQ=255)
                                         acATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccg                                                            >  1:3203696/1‑99 (MQ=255)
                                          cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccgcc                                                          >  2:318430/1‑100 (MQ=255)
                                          cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccgcc                                                          <  1:2593836/100‑1 (MQ=255)
                                          cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccgcc                                                          <  2:4352158/100‑1 (MQ=255)
                                          cATTAGTCGGCCCATCAGCCGAGGGGTACGGGAAAACTCTCTCTGTGCCAGCCAAAACTTGTGGCTGGAGTTGCCTGGTAACCGTCCACAACATgccgcc                                                          <  1:1382432/100‑1 (MQ=255)
                                           attattCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGccgccg                                                         <  1:1029545/100‑1 (MQ=255)
                                              attCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGaa                                                      >  1:1101226/1‑100 (MQ=255)
                                                tCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACg                                                    <  2:3757064/100‑1 (MQ=255)
                                                 cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGc                                                   >  1:1843450/1‑100 (MQ=255)
                                                 cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGc                                                   >  2:4210048/1‑100 (MQ=255)
                                                  aGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGcc                                                  >  1:2916692/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)
                                                    gCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCg                                                 >  1:2276106/1‑99 (MQ=255)
                                                      cATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTg                                              <  1:2772950/100‑1 (MQ=255)
                                                       aTCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTg                                              <  1:3569519/99‑1 (MQ=255)
                                                             ccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGacac                                       <  1:2190439/100‑1 (MQ=255)
                                                             ccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGacac                                       <  1:846556/100‑1 (MQ=255)
                                                                gCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTgg                                    >  1:3294374/1‑100 (MQ=255)
                                                                   gTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaaa                                 <  1:3199586/100‑1 (MQ=255)
                                                                      cGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                              >  2:3991323/1‑100 (MQ=255)
                                                                       gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                             >  2:2914737/1‑100 (MQ=255)
                                                                        ggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTg                            <  2:1709827/100‑1 (MQ=255)
                                                                              aTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGc                      <  2:1030231/100‑1 (MQ=255)
                                                                               tCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCg                     <  2:759947/100‑1 (MQ=255)
                                                                                  ctctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAg                  <  2:3066525/100‑1 (MQ=255)
                                                                                   tctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAgg                 <  2:3460355/100‑1 (MQ=255)
                                                                                     tctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAgg                 >  1:3864897/1‑98 (MQ=255)
                                                                                      nctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGttt              >  2:686140/4‑100 (MQ=255)
                                                                                      ctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGttt              >  1:3369736/1‑100 (MQ=255)
                                                                                      actccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGtt               >  1:1557223/4‑99 (MQ=255)
                                                                                       ctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTc             >  2:3381056/3‑100 (MQ=255)
                                                                                       ctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTc             >  2:2910619/3‑100 (MQ=255)
                                                                                       ctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGATGACACTGGAAAATCTGACTGGCGACGGTTTc             >  1:3772369/3‑100 (MQ=255)
                                                                                        tccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa            >  2:1924269/2‑100 (MQ=255)
                                                                                        gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCTCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa            <  1:2969128/100‑1 (MQ=255)
                                                                                        gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa            >  2:3998154/1‑100 (MQ=255)
                                                                                        gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa            >  2:2933674/1‑100 (MQ=255)
                                                                                         ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg           >  1:2012778/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)
                                                                                             ccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGACGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGagt       >  1:397030/1‑98 (MQ=255)
                                                                                              naaaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg      >  2:2030963/2‑100 (MQ=255)
                                                                                              caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg      >  2:2702740/1‑100 (MQ=255)
                                                                                              caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg      >  2:161464/1‑100 (MQ=255)
                                                                                              caAAAATTATGGCTGGAGTCACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg      >  1:1420054/1‑100 (MQ=255)
                                                                                                aaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTc    >  1:747377/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCAG  >  NC_000913/1602140‑1602337

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: