Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,491,906 A→T 5.0% intergenic (+56/+16) ydcA → / ← hokB putative periplasmic protein/toxic polypeptide, small

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,491,9060AT 5.0% 410.8 / 7.3 138intergenic (+56/+16)ydcA/hokBputative periplasmic protein/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base A (74/57);  new base T (4/3);  total (78/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 143 total aligned reads displayed.

GTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACACACCAC  >  NC_000913/1491756‑1492053
                                                                                                                                                      |                                                                                                                                                   
gTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTGTCTGTAATGATGGTTCTATTAGTGCATCGCAAAAAACATGCACTAAATTAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGaaaa                                                                                                                                                     >  1:945735/1‑150 (MQ=255)
 ttCGAAACTGGTTATCTCACACAGTACACCAGGAGGCAAATTTGTATGTAATTATGGTTCTATTAGTGCATAGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGaaaa                                                                                                                                                     >  2:668510/1‑150 (MQ=255)
      aaGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATcccc                                                                                                                                                <  1:723614/150‑1 (MQ=255)
      aaGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATcccc                                                                                                                                                >  2:1793868/1‑150 (MQ=255)
      aaGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATcccc                                                                                                                                                <  2:1894245/150‑1 (MQ=255)
       tgggtggATCTCACACTGTACGTCATGCGGAAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACAATCACTAACTGAAGTGTAAAATGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCcg                                                                                                                                               <  1:355903/149‑1 (MQ=255)
         gggggATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCcgcg                                                                                                                                             >  2:510617/1‑150 (MQ=255)
          ggggATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCgcgc                                                                                                                                            >  2:1061240/1‑150 (MQ=255)
           gggATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAATTGTAAAAGGGGTGCAATGAGAAAATGAATCTTATAGTCAAATTCGCGGGGAAAATCACAGCGCt                                                                                                                                           >  2:535722/1‑150 (MQ=255)
           gggATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCt                                                                                                                                           <  1:1089731/150‑1 (MQ=255)
           gggATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCt                                                                                                                                           >  1:1306555/1‑150 (MQ=255)
            ggATCTCACACTGTACGTCAGGCGGAAAATTTGTCTGTAATGATGGTTCTATTAGTGCATAGAAAAAAACATGCACTAACTGAAGTGTAAATTGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGATTTTCCCCgcgcat                                                                                                                                          <  2:1523102/150‑3 (MQ=255)
             gATCTCACACTGTAAGGCAGGCGTTTAATTTGTCTTTTATGATGGTGGTATTAGTGGATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGTGATTTTCCCCgcgcatt                                                                                                                                         <  2:128161/150‑4 (MQ=255)
               tctcACACTGTACGGAAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGcc                                                                                                                                       <  1:468001/150‑1 (MQ=255)
                ctcACACTGTACGTCAGGCGGCAAATTTGTCTGTAATGATGGGTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCATATGCGCGGGGATTTTCCCCgcgcatttga                                                                                                                                      <  1:964040/150‑7 (MQ=255)
                ctcACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGccc                                                                                                                                      >  2:1932807/1‑150 (MQ=255)
                ctcACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGccc                                                                                                                                      >  1:1682216/1‑150 (MQ=255)
                ctcACACTGTACGGCAGGCGGCAAATATGTCTGAAATCATGGTTCTATTAGAGCATCGAAAAAAACATGCACTAACTGAAGCGAAAAAGGGGTGCCATGAGAAAATGAAACTGCTAGACAAATGCGCGGGGAAAATCCCCGCGCTTGccc                                                                                                                                      <  1:259763/150‑1 (MQ=255)
                ctcACACTGTACGGAAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGccc                                                                                                                                      <  1:265157/150‑1 (MQ=255)
                        gTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTgg                                                                                                                              >  2:898769/1‑150 (MQ=255)
                          aCGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGAc                                                                                                                            >  2:485745/1‑150 (MQ=255)
                          aCGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGAc                                                                                                                            >  2:1000466/1‑150 (MQ=255)
                             gCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGTGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTg                                                                                                                         >  2:1123302/1‑150 (MQ=255)
                              cAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGc                                                                                                                        >  2:758830/1‑150 (MQ=255)
                                ggcggcTAATTTGTCTGTAATGATGCTTCTCTTAGTGCATCGAAAAAAACATCCATTAACTTAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAATATCCCCGCGCTTGCCCTTACCTGGACGTGCAg                                                                                                                      <  1:302008/150‑1 (MQ=255)
                                    gcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAACCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCa                                                                                                                  >  1:395113/1‑150 (MQ=255)
                                     cAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCAt                                                                                                                 >  2:961703/1‑150 (MQ=255)
                                     cAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCAt                                                                                                                 >  2:1860962/1‑150 (MQ=255)
                                        aTTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAg                                                                                                              >  2:1043228/1‑150 (MQ=255)
                                           tGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGACAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAgcgc                                                                                                           <  1:1804612/150‑1 (MQ=255)
                                            gTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgca                                                                                                          <  2:456016/150‑1 (MQ=255)
                                             tCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgcag                                                                                                         <  1:2015179/150‑1 (MQ=255)
                                              cTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgcagc                                                                                                        >  1:89949/1‑150 (MQ=255)
                                                    tgatgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAcctc                                                                                                  >  2:1578582/1‑150 (MQ=255)
                                                     gatgGTTCTTTTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAcctcc                                                                                                 <  1:1650066/150‑1 (MQ=255)
                                                      atgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCACcctcct                                                                                                >  1:1718864/1‑150 (MQ=255)
                                                       tgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCtt                                                                                               >  2:1598903/1‑150 (MQ=255)
                                                        ggTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGTTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCACCGCGCTTGCCCTTACCTGTACGTGCAGGCCATGAGCGCAGCAACATCCTTa                                                                                              >  2:204382/1‑150 (MQ=255)
                                                        ggTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTa                                                                                              <  1:802511/150‑1 (MQ=255)
                                                           tCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCa                                                                                           <  1:604859/150‑1 (MQ=255)
                                                           tCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCa                                                                                           <  1:837496/150‑1 (MQ=255)
                                                                 aGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCACCCTCCTTATCACCGTcc                                                                                     >  2:903242/1‑150 (MQ=255)
                                                                    gCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCgg                                                                                  >  1:1900102/1‑150 (MQ=255)
                                                                      aTCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAAGGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGCGCGCAGCAACCTCCTGATCACCGTCCCGGaa                                                                                >  1:1623269/1‑150 (MQ=255)
                                                                       tCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAAc                                                                               >  2:513926/1‑150 (MQ=255)
                                                                       tCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAAc                                                                               <  1:1300145/150‑1 (MQ=255)
                                                                        cGAAAAAAACATGCACTAACTGAAGTGTAAAAGGTGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAAcc                                                                              <  1:1128233/150‑1 (MQ=255)
                                                                        cGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCATTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAAcc                                                                              >  2:453369/1‑150 (MQ=255)
                                                                         gAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCg                                                                             >  2:558835/1‑150 (MQ=255)
                                                                         gAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCg                                                                             >  2:1090336/1‑150 (MQ=255)
                                                                         gAAAAAAACAGGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCACCCGCCTTATCACCGTCCCGGAAccc                                                                             >  1:1839785/1‑149 (MQ=255)
                                                                            aaaaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAg                                                                          <  1:552210/150‑1 (MQ=255)
                                                                             aaaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGCCGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCAGGAACCGCAGt                                                                         >  2:1758095/1‑150 (MQ=255)
                                                                             aaaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCAATGAGCGCAGCAACCTCCTTATCACAGTCCCGGACACGCAGt                                                                         >  1:1487157/1‑150 (MQ=255)
                                                                               aaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCCCCGTCCCGGAACCGCAGTTc                                                                       >  2:765778/1‑150 (MQ=255)
                                                                                 cATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGt                                                                     <  1:268643/150‑1 (MQ=255)
                                                                                  aTGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGGGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTa                                                                    <  2:799867/150‑1 (MQ=255)
                                                                                  aTGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTa                                                                    <  2:1370104/150‑1 (MQ=255)
                                                                                   tGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTag                                                                   >  2:972212/1‑150 (MQ=255)
                                                                                   tGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTag                                                                   <  1:1154277/150‑1 (MQ=255)
                                                                                       cTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCg                                                               <  2:508944/150‑1 (MQ=255)
                                                                                           cTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTg                                                           >  1:1470827/1‑150 (MQ=255)
                                                                                           cTGAAGTGTAAAAGGGGTGACATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTg                                                           <  2:1439521/150‑1 (MQ=255)
                                                                                               aGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCgg                                                       <  2:1516853/150‑1 (MQ=255)
                                                                                                   tAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCa                                                   <  1:1564476/150‑1 (MQ=255)
                                                                                                     aaaGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAgg                                                 >  2:908858/1‑150 (MQ=255)
                                                                                                     aaaGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAgg                                                 >  2:851876/1‑150 (MQ=255)
                                                                                                     aaaGGGGGGCCATGAGATAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGGAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAgg                                                 <  1:2010326/150‑1 (MQ=255)
                                                                                                       aGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAg                                               <  2:341071/150‑1 (MQ=255)
                                                                                                             gCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGaa                                         <  1:1061240/150‑1 (MQ=255)
                                                                                                               cATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATg                                       <  2:1722695/150‑1 (MQ=255)
                                                                                                                aTGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGt                                      >  2:1142551/1‑150 (MQ=255)
                                                                                                                 tGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTc                                     <  1:758830/150‑1 (MQ=255)
                                                                                                                   agaAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGGCGGGTCAGGAGTGTGAATGTCAg                                   <  2:1306555/150‑1 (MQ=255)
                                                                                                                     aaaaTGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGaa                                 <  1:513926/150‑1 (MQ=255)
                                                                                                                       aaTGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTCCCTGGACGTGAACGCCAAGAGCCCACAAACCTACTTATCACAGTCCCGGAACCTCAGTTCGTAGAGCGGTTGTCGGGTCAGGTGTGTGAATGTCAGAATc                               >  2:1018818/1‑150 (MQ=255)
                                                                                                                         tGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGt                             >  2:1438141/1‑150 (MQ=255)
                                                                                                                         tGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGGAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCATGATCGt                             >  1:1010015/1‑150 (MQ=255)
                                                                                                                           aaTCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTaa                           >  1:1273103/1‑150 (MQ=255)
                                                                                                                             tCTGCTAGTCTAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATg                         <  1:545776/150‑1 (MQ=255)
                                                                                                                             tCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATg                         >  1:1192730/1‑150 (MQ=255)
                                                                                                                                gCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAg                      <  2:1359151/150‑1 (MQ=255)
                                                                                                                                 cTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGa                     >  1:888240/1‑150 (MQ=255)
                                                                                                                                    gTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGGTTGTCGGGTCAGGAGTGTGAATGTCAGAAGCGTAATGCAGATaa                  >  1:282882/1‑150 (MQ=255)
                                                                                                                                      cAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATg                <  2:542268/150‑1 (MQ=255)
                                                                                                                                       gcaagcgcGGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGa               >  1:772356/5‑150 (MQ=255)
                                                                                                                                        caagcgcGGGGATATTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCCGCATCCTCCTATTCACCGTCCCGCTCCCGCTGTTCGTTGCGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGACTCCTAAGGCAGTTTATGAg              >  1:677989/4‑150 (MQ=255)
                                                                                                                                        aaTGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTACTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCTGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAg              <  1:1043228/150‑1 (MQ=255)
                                                                                                                                          tGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCa            <  2:1945065/150‑1 (MQ=255)
                                                                                                                                          agcgcGGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGCGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCa            >  1:1702517/2‑150 (MQ=255)
                                                                                                                                          agcgcGGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCa            >  2:1021074/2‑150 (MQ=255)
                                                                                                                                           gcgcGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGACCCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAg           >  1:50546/1‑150 (MQ=255)
                                                                                                                                           gcgcGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAg           <  2:1446943/150‑1 (MQ=255)
                                                                                                                                            cgcgGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGa          >  1:1389151/1‑150 (MQ=255)
                                                                                                                                                gggAAAATGCCCGCGGGTGCCCTTACCCGGCCCTGGATGCCATGAGGGCAGCAACCTCCGTATCCCCCGCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGacaca      <  2:517603/150‑1 (MQ=255)
                                                                                                                                                gggAAAATCCACGCGCTTGCCCTTACCTGGACGTGCAGGCCTTGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCTGAATCGTAATGCAGATAATGAGCAGacaca      <  1:1741313/150‑1 (MQ=255)
                                                                                                                                                 ggAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCACCCTCCTTATCACCGTCCCGGAGCCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGacccac     >  1:1399872/1‑150 (MQ=255)
                                                                                                                                                  tataaTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCAATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCACTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACAcacc    <  2:351964/147‑1 (MQ=255)
                                                                                                                                                  gAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCACCCTCCTTATCCCCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACAcacc    >  1:1740252/1‑150 (MQ=255)
                                                                                                                                                  gAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACAcacc    >  1:1622922/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACACACCAC  >  NC_000913/1491756‑1492053

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: