Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,602,239 A→T 16.8% 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,2390AT16.8% 343.3 / 58.9 142K250N (AAA→AATlsrAautoinducer 2 import ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (72/46);  new base T (8/16);  total (80/62)
Fisher's exact test for biased strand distribution p-value = 2.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

Only 100 of 147 total aligned reads displayed.

CATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCA  >  NC_000913/1602140‑1602336
                                                                                                   |                                                                                                 
 aTCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAAt                                                                                                  <  1:2928376/100‑1 (MQ=255)
 aTCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAAt                                                                                                  >  1:3308193/1‑100 (MQ=255)
  tCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCCGCTGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAAtt                                                                                                 >  2:39075/1‑100 (MQ=255)
   cGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTa                                                                                                <  1:2821660/100‑1 (MQ=255)
    gCCTTAAGCGGCAAAACCAGCGCACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTAt                                                                                               >  2:863872/1‑100 (MQ=255)
    gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTAt                                                                                               >  1:1157075/1‑100 (MQ=255)
    gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTAt                                                                                               >  1:3369466/1‑100 (MQ=255)
    gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTAt                                                                                               <  1:365774/100‑1 (MQ=255)
    gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATGAt                                                                                               >  2:1833491/1‑100 (MQ=255)
     ccTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATg                                                                                              <  2:524778/100‑1 (MQ=255)
      cTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAGTCGCTCTCTGCCAGCCAAAAATTATgg                                                                                             <  1:1222201/100‑1 (MQ=255)
      cTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATgg                                                                                             >  1:3555005/1‑100 (MQ=255)
       ttAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATg                                                                                              >  2:1579701/1‑98 (MQ=255)
       ttAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGc                                                                                            <  2:556829/100‑1 (MQ=255)
         aaGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTg                                                                                          <  1:1062100/100‑1 (MQ=255)
           gaggCAAAACCAGCGAACTGTCTACCGACGACATTATCCAGGCCATCACCCCAGAAGAAAGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                        <  2:1993397/98‑2 (MQ=255)
           gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                        <  1:881504/100‑2 (MQ=255)
           gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                        <  1:127486/100‑2 (MQ=255)
            cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggca                                                                                       <  1:892157/100‑3 (MQ=255)
            cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggca                                                                                       <  1:3504093/100‑3 (MQ=255)
            cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggca                                                                                       <  1:2581212/100‑3 (MQ=255)
            cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggca                                                                                       <  1:2270710/100‑3 (MQ=255)
             ggCAAAACCAGCGAACTGTCTACGGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGt                                                                                      >  1:1275811/1‑100 (MQ=255)
             ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCATGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                        <  2:223972/98‑2 (MQ=255)
             ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                      <  1:1009070/100‑4 (MQ=255)
             ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                      <  1:2709062/100‑4 (MQ=255)
             ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAAGCCACCACCCCAGCGGTACCGGAAAAATCCCTCACTGCCAGCCAAAAATTCTGGCTGGAg                                                                                       >  2:764712/1‑99 (MQ=255)
              gCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcaga                                                                                     <  1:3318019/100‑5 (MQ=255)
               cAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcagag                                                                                    <  1:3534005/100‑6 (MQ=38)
                 aaaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                                                                  <  1:3240627/100‑1 (MQ=255)
                 aaaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                                                                  >  1:1663037/1‑100 (MQ=255)
                 aaaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                                                                  <  1:1781867/100‑1 (MQ=255)
                  aaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAACTATGGCTGGAGTTACCt                                                                                 <  1:3358435/100‑1 (MQ=255)
                    cccccGAACTGTCCGCCGACGACATTATTCCGGTCATGGACCCAGCCGGACGTGTAAAATCGCTCGCTACCAGCCCACAATTTTCGCTGGAGTTcccgg                                                                                <  2:1887858/95‑3 (MQ=255)
                    ccAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTgg                                                                               <  1:593633/100‑1 (MQ=255)
                     cAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGACAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGGTACCTgg                                                                               >  1:3108909/1‑99 (MQ=255)
                     cAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGt                                                                              <  1:3342330/100‑1 (MQ=255)
                      aGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTa                                                                             <  2:1448952/100‑1 (MQ=255)
                       gCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTaa                                                                            >  1:1637528/1‑100 (MQ=255)
                              gTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCca                                                                     >  2:3221203/1‑100 (MQ=255)
                               tCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCcac                                                                    <  1:3011852/100‑1 (MQ=255)
                                cTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCcaca                                                                   >  2:187903/1‑100 (MQ=255)
                                 tACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaa                                                                  >  2:3468816/1‑100 (MQ=255)
                                 tACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaa                                                                  >  1:1310969/1‑100 (MQ=255)
                                  aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                                 >  2:1801572/1‑100 (MQ=255)
                                  aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                                 >  2:946331/1‑100 (MQ=255)
                                   ccgatctcATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaaca                                                                <  1:27950/93‑1 (MQ=255)
                                        gacATTATTCACGCCATCCCCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccg                                                           <  2:3617628/100‑1 (MQ=255)
                                         acATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccac                                                          >  2:2201026/1‑98 (MQ=255)
                                          cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaacaagatcgg                                                         >  2:27950/1‑93 (MQ=255)
                                          cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccgcc                                                         >  1:2030823/1‑100 (MQ=255)
                                           attattCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGccgccg                                                        <  1:1658422/100‑1 (MQ=255)
                                            ttattCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGccgccg                                                        >  2:3359003/1‑99 (MQ=255)
                                            ttattCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCgg                                                       >  2:3167466/1‑100 (MQ=255)
                                             tattCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGa                                                      >  1:1635674/1‑100 (MQ=255)
                                             tattCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGa                                                      <  1:1544675/100‑1 (MQ=255)
                                               ttCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAAc                                                    >  1:1870078/1‑100 (MQ=255)
                                                 cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGc                                                  >  1:2409228/1‑100 (MQ=255)
                                                   ggCCATCACCCCAGCGGTACGGGAACAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCg                                                >  2:282908/1‑100 (MQ=255)
                                                    gCCATCACCCCAGCGGTACGGGCAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTAGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCgg                                               >  2:2328157/1‑100 (MQ=255)
                                                    gCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCgg                                               >  1:779454/1‑100 (MQ=255)
                                                     ccATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGt                                              >  2:2426108/1‑100 (MQ=255)
                                                      cATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTg                                             >  2:2234954/1‑100 (MQ=255)
                                                       aTCACCCCAGCGGTACGGGAAAAATCGCTCTCTGACAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGc                                            <  2:3204528/100‑1 (MQ=255)
                                                           ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGac                                        >  1:3284254/1‑100 (MQ=255)
                                                                  ggTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaa                                 >  1:686170/1‑100 (MQ=255)
                                                                   gTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaaa                                >  1:1432893/1‑100 (MQ=255)
                                                                     aCGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACCCTGGAAAAt                              >  1:2753333/1‑100 (MQ=255)
                                                                     aCGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAAt                              >  1:216600/1‑100 (MQ=255)
                                                                     aCGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAAt                              <  2:1816061/100‑1 (MQ=255)
                                                                      cGGGAAGAATCGCTTTCTGCCAGCCAAAAATTATGGCTGGTGTTACCTGGTAACCGCCCACATCTTGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                             <  2:1110977/100‑1 (MQ=255)
                                                                       gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                            >  1:1366490/1‑100 (MQ=255)
                                                                       gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                            >  1:983497/1‑100 (MQ=255)
                                                                        ggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTg                           <  1:2374216/100‑1 (MQ=255)
                                                                        cgAACAATCCATCTCTGCCACCCAAAAATTATCCCTGGACTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                             <  2:340395/97‑1 (MQ=255)
                                                                          aacaaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCAACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGAc                         <  2:1207349/97‑1 (MQ=255)
                                                                          aaaaaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGAc                         >  1:1232262/1‑100 (MQ=255)
                                                                            aaaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCg                       >  1:3207858/1‑100 (MQ=255)
                                                                             aaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCgg                      <  1:2475166/100‑1 (MQ=255)
                                                                               tCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCg                    >  1:1481390/1‑100 (MQ=255)
                                                                                cGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGa                   <  1:290016/100‑1 (MQ=255)
                                                                                  ctctctGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTGACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAg                 <  2:2745714/100‑1 (MQ=255)
                                                                                    taactccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAgg                >  2:763081/6‑99 (MQ=38)
                                                                                      ctGCGAGCTAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGGCGGTGCTGACACTGGAACATCTGACCGGCGAAg                 >  2:1958357/1‑96 (MQ=255)
                                                                                       tGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGttt             >  1:3210548/1‑99 (MQ=255)
                                                                                       ctccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGgtt             >  2:3038766/3‑99 (MQ=255)
                                                                                        gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTc            >  1:850072/1‑99 (MQ=255)
                                                                                        gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa           <  2:1181088/100‑1 (MQ=255)
                                                                                         ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg          >  1:1631569/1‑100 (MQ=255)
                                                                                         ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg          >  1:1465776/1‑100 (MQ=255)
                                                                                         ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg          >  1:1202145/1‑100 (MQ=255)
                                                                                          cagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAgg         <  1:3167489/100‑1 (MQ=255)
                                                                                          cagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAgg         <  2:2693954/100‑1 (MQ=255)
                                                                                             ccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAAt      >  1:2205158/1‑100 (MQ=255)
                                                                                              caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg     >  1:2990970/1‑100 (MQ=255)
                                                                                              caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg     >  2:2790966/1‑100 (MQ=255)
                                                                                              caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg     >  1:2619095/1‑100 (MQ=255)
                                                                                                aaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCTGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTc   <  2:3431583/100‑1 (MQ=255)
                                                                                                aaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTc   <  2:2024450/100‑1 (MQ=255)
                                                                                                 aaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCa  >  1:1091510/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
CATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTCA  >  NC_000913/1602140‑1602336

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: