Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,401,582 A→T 6.1% L584Q (CTG→CAG)  csrD ← regulator of CsrB and CsrC decay

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,401,5820AT 6.1% 296.6 / 4.3 93L584Q (CTG→CAG) csrDregulator of CsrB and CsrC decay
Reads supporting (aligned to +/- strand):  ref base A (42/45);  new base T (3/3);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTCACTA  >  NC_000913/3401433‑3401729
                                                                                                                                                     |                                                                                                                                                   
cAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagca                                                                                                                                                     <  2:926652/150‑1 (MQ=255)
cAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagca                                                                                                                                                     <  1:786886/150‑1 (MQ=255)
cAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagca                                                                                                                                                     <  1:1592372/150‑1 (MQ=255)
   gTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCt                                                                                                                                                  >  1:41090/1‑150 (MQ=255)
        tgtgAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTtc                                                                                                                                             <  1:363745/150‑1 (MQ=255)
          tgAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTtctc                                                                                                                                           >  1:1264932/1‑150 (MQ=255)
           gAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTcc                                                                                                                                          <  2:1608310/150‑1 (MQ=255)
             ggACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCTGCTGGTtcaaagc                                                                                                                                        <  1:901322/150‑6 (MQ=255)
               aCGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTc                                                                                                                                      >  2:1530979/1‑150 (MQ=255)
               aCGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTc                                                                                                                                      <  1:783908/150‑1 (MQ=255)
                 gcgAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGc                                                                                                                                    >  2:524761/1‑150 (MQ=255)
                  cgAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCt                                                                                                                                   >  2:1661330/1‑150 (MQ=255)
                  cgAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCt                                                                                                                                   >  2:1937722/1‑150 (MQ=255)
                  cgAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCt                                                                                                                                   >  1:989251/1‑150 (MQ=255)
                   gAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAAGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCtt                                                                                                                                  >  2:184104/1‑150 (MQ=255)
                     aaaaaTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTtct                                                                                                                                <  1:944775/150‑1 (MQ=255)
                      aaaaTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGGCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTtctc                                                                                                                               <  2:140902/150‑1 (MQ=255)
                       aaaTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCTGCTGGTtc                                                                                                                                             >  2:1183196/1‑135 (MQ=255)
                       aaaTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCTGCTGGTtc                                                                                                                                             <  1:1183196/135‑1 (MQ=255)
                        aaTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCaa                                                                                                                             >  1:1484678/1‑150 (MQ=255)
                        aaTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCaa                                                                                                                             <  1:1306038/150‑1 (MQ=255)
                           ccccTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGt                                                                                                                          >  2:158988/1‑150 (MQ=255)
                             ccTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGttt                                                                                                                        >  1:915077/1‑150 (MQ=255)
                                 gCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGa                                                                                                                    <  2:691432/150‑1 (MQ=255)
                                   ccGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGAcc                                                                                                                  <  1:1495196/150‑1 (MQ=255)
                                     gCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAg                                                                                                                <  2:1482675/150‑1 (MQ=255)
                                     gCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAg                                                                                                                >  2:1934197/1‑150 (MQ=255)
                                             cacCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGAt                                                                                                        >  1:830973/1‑150 (MQ=255)
                                              acCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATg                                                                                                       >  1:1582437/1‑150 (MQ=255)
                                                  cgcTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGc                                                                                                   <  1:1347896/150‑1 (MQ=255)
                                                  cgcTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGc                                                                                                   <  2:1264932/150‑1 (MQ=255)
                                                     tGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTg                                                                                                <  1:524761/150‑1 (MQ=255)
                                                     tGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTg                                                                                                <  2:134248/150‑1 (MQ=255)
                                                      gAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGa                                                                                               <  2:424881/150‑1 (MQ=255)
                                                       aaTCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAg                                                                                              >  2:565075/1‑150 (MQ=255)
                                                        aTCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGt                                                                                             >  2:485455/1‑150 (MQ=255)
                                                         tgaGGGTCTGCCACTCGGTTCGCGAACGCACGCCGGTGGCGTAAACCGGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTa                                                                                            <  2:863348/148‑1 (MQ=255)
                                                         tCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTa                                                                                            >  2:691060/1‑150 (MQ=255)
                                                           aGGTTCTGCCACTCGCTTCGCGAACGCACGCAGGTGGCGTAAACCTGGTTTCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATTGAGCTTGAGTAAc                                                                                          <  2:119826/150‑1 (MQ=255)
                                                              gTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCa                                                                                       <  2:92565/150‑1 (MQ=255)
                                                               tCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCaa                                                                                      <  2:1484678/150‑1 (MQ=255)
                                                                 tGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACa                                                                                    <  1:1815053/150‑1 (MQ=255)
                                                                  gCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACAt                                                                                   >  1:1302836/1‑150 (MQ=255)
                                                                  gCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACAt                                                                                   <  2:989251/150‑1 (MQ=255)
                                                                    cACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTa                                                                                 >  1:49753/1‑150 (MQ=255)
                                                                      cTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAg                                                                               <  2:41090/150‑1 (MQ=255)
                                                                       tCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGt                                                                              <  1:350706/150‑1 (MQ=255)
                                                                          cTTCGCGAACGCACGCCGGTGGCGTAAACCCGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCt                                                                           >  2:554945/1‑150 (MQ=255)
                                                                           ttCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCATGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCtt                                                                          <  1:184104/150‑1 (MQ=255)
                                                                            tCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCttt                                                                         <  1:1024560/150‑1 (MQ=255)
                                                                            tCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCttt                                                                         <  1:180294/150‑1 (MQ=255)
                                                                            tCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCttt                                                                         <  1:1937722/150‑1 (MQ=255)
                                                                              gcgAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGa                                                                       >  2:1183894/1‑150 (MQ=255)
                                                                               cgAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGAt                                                                      <  1:169500/150‑1 (MQ=255)
                                                                                gAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATc                                                                     >  1:1906880/1‑150 (MQ=255)
                                                                                 aaCGCACGCCGGTGGCGTAAACCTGGGTGCTGGTGCCGGAGCAGTCTTCCACCAGGCTTTGAACCAGCTGCTGGTTCTCCTTTCGCTTCTCTATGTTTCTGACCAGCCCCGTATGGAGCTTGATTAACTCAACATTAAGTTCTTTGATcc                                                                    <  1:1350990/150‑1 (MQ=255)
                                                                                 aaCGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATcc                                                                    <  2:1106248/150‑1 (MQ=255)
                                                                                    gcacgcCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAAc                                                                 >  1:1789012/1‑150 (MQ=255)
                                                                                        gcCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGt                                                             <  1:444301/150‑1 (MQ=255)
                                                                                        gcCGGTGGCGTAAACCTGGGTGCTGGGCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGt                                                             <  1:1003511/150‑1 (MQ=255)
                                                                                              ggCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTAc                                                       <  1:151648/150‑1 (MQ=255)
                                                                                                    aaCCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGt                                                 >  2:1609664/1‑150 (MQ=255)
                                                                                                      ccTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCa                                               >  2:1244354/1‑150 (MQ=255)
                                                                                                               tGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCt                                      >  1:884067/1‑150 (MQ=255)
                                                                                                               tGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCt                                      >  1:802889/1‑150 (MQ=255)
                                                                                                                ggTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTg                                     >  2:560796/1‑150 (MQ=255)
                                                                                                                ggTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACAAGCGTCAAACCAGCCTg                                     >  2:1468949/1‑150 (MQ=255)
                                                                                                                  tCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGt                                   >  1:606544/1‑150 (MQ=255)
                                                                                                                      ggAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGAc                               >  2:1903411/1‑150 (MQ=255)
                                                                                                                        aGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACgg                             <  1:1810731/150‑1 (MQ=255)
                                                                                                                        aGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACgg                             <  2:1906880/150‑1 (MQ=255)
                                                                                                                         gCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGc                            <  1:691060/150‑1 (MQ=255)
                                                                                                                         gCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGc                            <  1:485455/150‑1 (MQ=255)
                                                                                                                         gCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAACCAGGCTGGtttgacgct                           <  2:1302836/150‑8 (MQ=255)
                                                                                                                            ggCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAc                         <  1:1609664/150‑1 (MQ=255)
                                                                                                                            ggCTTCAACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAc                         >  2:449881/1‑150 (MQ=255)
                                                                                                                             gCTTCCACCAGGCTTTGAACCAGNAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGNGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAcc                        <  2:514454/150‑1 (MQ=255)
                                                                                                                             gCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAcc                        >  1:1933171/1‑150 (MQ=255)
                                                                                                                             gCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAcc                        <  2:507264/150‑1 (MQ=255)
                                                                                                                              cTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAccc                       >  2:1018665/1‑150 (MQ=255)
                                                                                                                              cTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAccc                       <  1:554945/150‑1 (MQ=255)
                                                                                                                                 ccaccaGGCGTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTa                    >  2:1204332/1‑150 (MQ=255)
                                                                                                                                  caccaGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTAc                   >  2:459200/1‑150 (MQ=255)
                                                                                                                                    ccaGGCTTTGAACCATCATCTGGTTCTCCGTTCGCTTCTCAATTTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGTTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTAccc                 <  2:784920/150‑1 (MQ=255)
                                                                                                                                    ccaGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTAccc                 >  2:1634716/1‑150 (MQ=255)
                                                                                                                                    ccaGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTAccc                 >  1:682955/1‑150 (MQ=255)
                                                                                                                                     caGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAc                         <  2:1886034/141‑1 (MQ=255)
                                                                                                                                     caGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTAc                         >  1:1886034/1‑141 (MQ=255)
                                                                                                                                     caGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTAcccc                <  2:830973/150‑1 (MQ=255)
                                                                                                                                        acggagAACCAGCTGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTaa             >  2:1847053/6‑150 (MQ=255)
                                                                                                                                         cTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTaaa            >  2:972974/1‑150 (MQ=255)
                                                                                                                                          tttGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAg           <  1:324859/150‑1 (MQ=255)
                                                                                                                                          tttGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAg           <  1:349524/150‑1 (MQ=255)
                                                                                                                                          ggagAACCAGCTGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAg           >  1:384866/4‑150 (MQ=255)
                                                                                                                                               aCCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTc      >  2:1656209/1‑150 (MQ=255)
                                                                                                                                                ccAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTCa     >  1:942425/1‑150 (MQ=255)
                                                                                                                                                 cagcagCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTCAc    <  1:1661330/150‑1 (MQ=255)
                                                                                                                                                 cagcagCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTCAc    <  1:1244354/150‑1 (MQ=255)
                                                                                                                                                  agcagcTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTCACt   >  1:961752/1‑150 (MQ=255)
                                                                                                                                                   gcagcTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTCACTa  <  1:720186/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTTCGCTTCTCAATGTTTCTGACCAGCCCCGGATGGAGCTTGAGTAACTCAACATTAAGTTCTTTGATCCAACTGGTACTTACCAGCGTCAAACCAGCCTGGTTGACGGCTACCCGTACCCCTAAAGCATTCACTA  >  NC_000913/3401433‑3401729

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