Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,602,236 A→T 16.8% 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,2360AT16.8% 342.8 / 60.7 142Q249H (CAA→CATlsrAautoinducer 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 = 2.13e-01

Only 100 of 148 total aligned reads displayed.

ACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTC  >  NC_000913/1602138‑1602335
                                                                                                  |                                                                                                   
aCCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaa                                                                                                    <  2:2911973/100‑1 (MQ=255)
  cATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCaaaaa                                                                                                  >  1:3551829/1‑100 (MQ=255)
   aTCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAAt                                                                                                 >  1:3308193/1‑100 (MQ=255)
   aTCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAAt                                                                                                 <  1:2928376/100‑1 (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:365774/100‑1 (MQ=255)
      gCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTAt                                                                                              >  1:3369466/1‑100 (MQ=255)
       ccTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATg                                                                                             <  2:524778/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)
            aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTgg                                                                                        <  1:1290170/100‑1 (MQ=255)
             gCGGCAAAACCAGTGAACTGTCTACCGACGACATAATGCGGGCTATCACCCCAGCGGTGCGGGAAAAATCGCTCTCTGCCAGCAATAATTTTTGGCTggc                                                                                       <  2:3201804/100‑2 (MQ=255)
             gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                       <  1:881504/100‑2 (MQ=255)
              cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggca                                                                                      <  1:3504093/100‑3 (MQ=255)
              cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggca                                                                                      <  1:892157/100‑3 (MQ=255)
              cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggca                                                                                      <  1:2270710/100‑3 (MQ=255)
              cGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggca                                                                                      <  1:2581212/100‑3 (MQ=255)
               ggCAAAACCAGCGAACTGTCTACGGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGt                                                                                     >  1:1275811/1‑100 (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)
                   aaaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                                                                 <  1:3240627/100‑1 (MQ=255)
                   aaaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                                                                 >  1:1663037/1‑100 (MQ=255)
                   aaaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                                                                 <  1:1781867/100‑1 (MQ=255)
                   aaaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTAcc                                                                                 >  1:1479352/1‑100 (MQ=255)
                    aaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcagagagc                                                                                <  1:2753069/100‑9 (MQ=37)
                    aaCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAACTATGGCTGGAGTTACCt                                                                                <  1:3358435/100‑1 (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                                                                 >  1:1207357/1‑100 (MQ=255)
                                   tACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaa                                                                 >  2:3468816/1‑100 (MQ=255)
                                    aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                                >  2:946331/1‑100 (MQ=255)
                                    aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                                >  2:231935/1‑100 (MQ=255)
                                    aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                                >  2:1801572/1‑100 (MQ=255)
                                     ccgatctcATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaaca                                                               <  1:27950/93‑1 (MQ=255)
                                     ccGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaaca                                                               >  2:2216265/1‑100 (MQ=255)
                                          gacATTATTCACGCCATCCCCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccg                                                          <  2:3617628/100‑1 (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)
                                                attCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGaa                                                    >  2:3233576/1‑100 (MQ=255)
                                                 ttCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAAc                                                   >  1:1870078/1‑100 (MQ=255)
                                                  tCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACg                                                  <  1:2429171/100‑1 (MQ=255)
                                                   cAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGc                                                 >  1:2409228/1‑100 (MQ=255)
                                                     ggCCATCACCCCAGCGGTACGGGAACAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCg                                               >  2:282908/1‑100 (MQ=255)
                                                     ggCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCg                                               >  1:2600293/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)
                                                       ccATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGt                                             >  1:1262669/1‑100 (MQ=255)
                                                        cATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTg                                            >  2:2234954/1‑100 (MQ=255)
                                                          tCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCt                                          <  2:189796/100‑1 (MQ=255)
                                                          tCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCt                                          <  2:1851877/100‑1 (MQ=255)
                                                             ccccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGac                                       >  1:3284254/1‑100 (MQ=255)
                                                                    ggTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGaa                                >  1:686170/1‑100 (MQ=255)
                                                                       aCGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAAt                             <  2:1816061/100‑1 (MQ=255)
                                                                       aCGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAAt                             >  1:216600/1‑100 (MQ=255)
                                                                        cGGGAAGAATCGCTTTCTGCCAGCCAAAAATTATGGCTGGTGTTACCTGGTAACCGCCCACATCTTGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                            <  2:1110977/100‑1 (MQ=255)
                                                                         gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                           >  1:983497/1‑100 (MQ=255)
                                                                         gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                           >  1:1366490/1‑100 (MQ=255)
                                                                          ggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTg                          <  1:2374216/100‑1 (MQ=255)
                                                                          ggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTg                          >  1:1109847/1‑100 (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)
                                                                                  cGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGa                  <  1:290016/100‑1 (MQ=255)
                                                                                  cGCTCTCTGCCAGCCAAAAATTATGGCGGGAGTTACCGGGGAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGa                  <  2:2008277/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)
                                                                                          gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTc           >  1:850072/1‑99 (MQ=255)
                                                                                          gccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCa          <  2:1181088/100‑1 (MQ=255)
                                                                                           ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg         >  1:1202145/1‑100 (MQ=255)
                                                                                           ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg         >  1:1465776/1‑100 (MQ=255)
                                                                                            cagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAgg        <  1:3167489/100‑1 (MQ=255)
                                                                                            cagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAgg        <  2:2693954/100‑1 (MQ=255)
                                                                                               ccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAAt     >  1:1548092/1‑100 (MQ=255)
                                                                                               ccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAAt     >  1:2205158/1‑100 (MQ=255)
                                                                                                caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg    >  2:2790966/1‑100 (MQ=255)
                                                                                                caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg    >  1:2990970/1‑100 (MQ=255)
                                                                                                caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg    >  1:1387101/1‑100 (MQ=255)
                                                                                                  aaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCTGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTc  <  2:3431583/100‑1 (MQ=255)
                                                                                                  aaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTc  <  2:2024450/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
ACCATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGTC  >  NC_000913/1602138‑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: