Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,319,582 C→G 13.5% V99L (GTG→CTG) ‡ rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,319,5820CG13.5% 215.7 / 23.7 81V99L (GTG→CTG) ‡rcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base C (41/28);  new base G (8/3);  total (49/32)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

GATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGAC  >  NC_000913/2319433‑2319731
                                                                                                                                                     |                                                                                                                                                     
gATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGcac                                                                                                                                                       >  2:1643447/1‑150 (MQ=255)
 aTTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGcaca                                                                                                                                                      >  1:1138427/1‑150 (MQ=255)
     aaaaTACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCg                                                                                                                                                  >  2:962736/1‑150 (MQ=255)
     aaaaTACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTAAAACGCAGGCACATCCg                                                                                                                                                  >  2:87885/1‑150 (MQ=255)
           cccGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCg                                                                                                                                            <  1:743963/150‑1 (MQ=255)
            ccGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCGTAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGt                                                                                                                                           <  2:961942/150‑1 (MQ=255)
               gTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtc                                                                                                                                        <  2:1579901/150‑1 (MQ=255)
                tttAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtct                                                                                                                                       >  1:1044767/1‑150 (MQ=255)
                  tAGATCTTAAGACGCATAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCATTTACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCg                                                                                                                                     <  1:1362200/150‑1 (MQ=255)
                  tAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCg                                                                                                                                     <  1:1730000/150‑1 (MQ=255)
                   aGATCGTACGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGt                                                                                                                                    <  1:1253336/150‑1 (MQ=255)
                   aGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGtttgaacc                                                                                                                                    <  1:1609972/150‑6 (MQ=255)
                   aGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACATCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCGGGGACATCCGCCTGCGTGTCTCGt                                                                                                                                    >  1:205323/1‑150 (MQ=255)
                        gTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcg                                                                                                                               >  2:902562/1‑150 (MQ=255)
                              cGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAGCCTCGCCAGGTGTTACTCATAGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATACGCCTGCGTTTCTCGTCCACGCGGGGa                                                                                                                         >  1:1534403/1‑150 (MQ=255)
                               gCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCCCGCGGGGaa                                                                                                                        >  1:955862/1‑150 (MQ=255)
                                cAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGACCAATCGGAGTCGGCACACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGacg                                                                                                                       >  1:208515/1‑148 (MQ=255)
                                 aGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAg                                                                                                                      >  1:411207/1‑150 (MQ=255)
                                      aaTTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCacac                                                                                                                 <  1:585620/150‑1 (MQ=255)
                                         tATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGt                                                                                                              >  1:1728964/1‑150 (MQ=255)
                                             gcgcCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTc                                                                                                          >  1:1116591/1‑150 (MQ=255)
                                              cgcCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACTCGGGGAAAGCAAACCGTTTTCt                                                                                                         >  1:142584/1‑150 (MQ=255)
                                               gcCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTg                                                                                                        >  2:543447/1‑150 (MQ=255)
                                                  aGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGAGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCCCATCCGCCTGCGTTTCTCGTCCCCGCGGGGAAAGCACACCGTTTTCTGCCg                                                                                                     >  2:1101816/1‑150 (MQ=255)
                                                            aaaCCACGCCAATGCCTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGACTGCGGTTCTAGTCCACGAGGGGAAAGCACCCCGTTTTATGCAGATAAGCGAtt                                                                                           >  2:1305427/1‑150 (MQ=255)
                                                            aaaCCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGAtt                                                                                           >  2:657612/1‑150 (MQ=255)
                                                             aaCCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGAttt                                                                                          >  2:203841/1‑150 (MQ=255)
                                                             aaCCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGAGAttt                                                                                          >  1:1439519/1‑150 (MQ=255)
                                                                 aCGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCgg                                                                                      >  1:1371730/1‑150 (MQ=255)
                                                                  cGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCCCACCGTTTTCTGCCGATAAGCGATTTTCGGc                                                                                     >  1:1318666/1‑150 (MQ=255)
                                                                      aaTGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGAAAATCGCTTATCGGcagaa                                                                                 <  2:122995/150‑5 (MQ=255)
                                                                        tGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTa                                                                               <  1:936998/150‑1 (MQ=255)
                                                                          aCTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACt                                                                             >  2:1702488/1‑150 (MQ=255)
                                                                              cAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGAAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAg                                                                         >  1:392625/1‑150 (MQ=255)
                                                                               agagAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGc                                                                        >  2:670649/1‑150 (MQ=255)
                                                                                gagaACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGct                                                                       <  2:1559812/150‑1 (MQ=255)
                                                                                   aaCCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCtt                                                                    <  2:174321/150‑1 (MQ=255)
                                                                                    aCCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGGTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCttt                                                                   >  1:1273927/1‑150 (MQ=255)
                                                                                      cTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGGTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTtca                                                                 <  1:1215851/150‑1 (MQ=255)
                                                                                        cGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGAGATTTTCGGCGATGTACTTCATCTCTTtcatc                                                               >  2:1265901/1‑150 (MQ=255)
                                                                                          ccAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTcatcac                                                             <  2:1264347/150‑1 (MQ=255)
                                                                                               tgtTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCag                                                        <  2:191506/150‑1 (MQ=255)
                                                                                                  tACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGTCGATGTACTTCAGCTCTTTCATCACATCAGAgg                                                     <  1:1565709/150‑1 (MQ=255)
                                                                                                     tCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGc                                                  <  2:1116591/150‑1 (MQ=255)
                                                                                                       aTTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCCCATCCGCCTGCGGTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGAGGTCCTTCAGCTCGTTCATCACATCAGAGGTGCgg                                                >  2:141935/1‑149 (MQ=255)
                                                                                                         tGCGGAACAATCGGAGTCGGCAAACCGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTg                                              <  2:840721/150‑1 (MQ=255)
                                                                                                         tGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTg                                              <  2:547088/150‑1 (MQ=255)
                                                                                                         tGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTg                                              <  1:203841/150‑1 (MQ=255)
                                                                                                          gCGGAACAATCGGTGTCGGCAATCAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGGCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGa                                             >  2:224272/1‑150 (MQ=255)
                                                                                                          gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTgt                                             >  1:1675896/1‑149 (MQ=255)
                                                                                                          gCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGa                                             >  2:656897/1‑150 (MQ=255)
                                                                                                           cGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATCTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGaa                                            >  1:889821/1‑150 (MQ=255)
                                                                                                            ggAACAAGCGGAGGCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGGTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAAt                                           <  2:848590/150‑1 (MQ=255)
                                                                                                               aCAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATaaa                                        <  1:902562/150‑1 (MQ=255)
                                                                                                                cAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAg                                       <  2:913219/150‑1 (MQ=255)
                                                                                                                 aaTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTATCGTCCACGCGGGGAAAGCACACCGTTTTATGCCGATAGGCGATTTTCGGCGATGTACTTCCGCTCTTTAATCACATCAGAGGTGCGTTGAATAACgc                                      >  1:913496/1‑150 (MQ=255)
                                                                                                                  aTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGTGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCATCTCTTTCATCACATCAGAGGTGCGTTGAATAAcgcg                                     >  2:843528/1‑150 (MQ=255)
                                                                                                                  aTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAgcg                                     >  2:41967/1‑150 (MQ=255)
                                                                                                                    ctgAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCTTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCTTTGAATAAAGCGCt                                   <  2:1228927/148‑1 (MQ=255)
                                                                                                                     ggAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTg                                  >  2:1268899/1‑150 (MQ=255)
                                                                                                                        gTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCCCGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCCGCTCTTTCATCACATCAGAGGTGCGTTGAATAACGCGCTgcgc                               >  2:1226381/1‑150 (MQ=255)
                                                                                                                        gTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGc                               <  2:1728964/150‑1 (MQ=255)
                                                                                                                          cGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCt                             >  1:1412322/1‑150 (MQ=255)
                                                                                                                           ggCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTg                            >  1:943473/1‑150 (MQ=255)
                                                                                                                            gCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGa                           >  1:769230/1‑150 (MQ=255)
                                                                                                                                 cAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACCCCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCGTTCATCACATCAGAGGTGCGTTGAAGAAAGCGCTGAGCCTGAGCGGa                      >  1:476525/1‑150 (MQ=255)
                                                                                                                                 cAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGa                      <  1:543447/150‑1 (MQ=255)
                                                                                                                                         gaaaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGAt              >  2:393332/2‑150 (MQ=255)
                                                                                                                                         gaaaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGAt              >  1:566756/2‑150 (MQ=255)
                                                                                                                                         gaaaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTACCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGAt              >  2:474828/2‑150 (MQ=255)
                                                                                                                                         cAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGAt              <  1:962736/150‑1 (MQ=255)
                                                                                                                                          aaaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTGTCTGCGGATGAGCGCTGTTCGGCGATGTACGTCAGATATTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGACATCAGAtt             >  1:281383/1‑150 (MQ=255)
                                                                                                                                          aaaCGCAGGCACATCCGCGTGAGTTTCTCGTCTACGAGGGGATAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTTCTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGAtt             <  1:166227/150‑1 (MQ=255)
                                                                                                                                          aaaCGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCTATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGAtt             <  1:1265901/150‑1 (MQ=255)
                                                                                                                                           aaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCCCCCCGTTTACTGCCGATAAGCGATTTTCGGCGATGTACTTCAGATGTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTa            >  2:437226/1‑150 (MQ=255)
                                                                                                                                           aaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTa            >  1:736563/1‑150 (MQ=255)
                                                                                                                                           aaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTATTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTa            >  2:1127858/1‑150 (MQ=255)
                                                                                                                                           aaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACAGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTa            >  1:376690/1‑150 (MQ=255)
                                                                                                                                           aaCGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTa            <  1:22032/150‑1 (MQ=255)
                                                                                                                                           aaCGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTa            >  1:1215456/1‑150 (MQ=255)
                                                                                                                                           aaCGCAGCCACATCCGCCTGCGTTTCTCGTCCACGCCGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCTATGTACTTCAGCTCTTTCATCACATCACAGGTGCGTTGAATAAAGCTCTGAGCCTGATCGGAACTCAGAtaa            >  2:1619305/1‑148 (MQ=255)
                                                                                                                                              gCAGGCAAATCCGCCTGCGTTTCTCGTCCACGCGGTGAAAGCACACCGTTTTCTGCCGATAAGCTATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAAt         <  2:1534403/150‑1 (MQ=255)
                                                                                                                                               cAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAAtt        <  1:656897/150‑1 (MQ=255)
                                                                                                                                                     caTCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGAc  >  2:438425/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGAC  >  NC_000913/2319433‑2319731

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: