Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,103,526 A→G 5.0% intergenic (+30/+76) yiiM → / ← cpxA 6‑N‑hydroxylaminopurine resistance protein/sensory histidine kinase in two‑component regulatory system with CpxR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,103,5260AG 5.0% 358.3 / 6.8 119intergenic (+30/+76)yiiM/cpxA6‑N‑hydroxylaminopurine resistance protein/sensory histidine kinase in two‑component regulatory system with CpxR
Reads supporting (aligned to +/- strand):  ref base A (61/52);  new base G (2/4);  total (63/56)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 121 total aligned reads displayed.

GATCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCC  >  NC_000913/4103377‑4103675
                                                                                                                                                     |                                                                                                                                                     
gATCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGATAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTa                                                                                                                                                       >  1:456617/1‑150 (MQ=255)
 aTCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTAt                                                                                                                                                      >  1:603003/1‑150 (MQ=255)
 aTCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTAt                                                                                                                                                      >  2:1691307/1‑150 (MQ=255)
      tatCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCg                                                                                                                                                 >  1:1598729/1‑150 (MQ=255)
           ccGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTGATAGCGTctta                                                                                                                                            <  1:1603534/150‑4 (MQ=255)
             gTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTGATAGCGTcttatc                                                                                                                                          <  2:1330223/150‑6 (MQ=255)
               ttACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAg                                                                                                                                        <  2:603003/150‑1 (MQ=255)
                tACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAgg                                                                                                                                       <  2:1592532/150‑1 (MQ=255)
                tACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAgg                                                                                                                                       <  2:742476/150‑1 (MQ=255)
                 aCTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGc                                                                                                                                      >  1:1057632/1‑150 (MQ=255)
                  ctctCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCa                                                                                                                                     >  2:827767/1‑150 (MQ=255)
                    gtcCGCTGCCGGGTTATCGAAATGGTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATc                                                                                                                                   <  2:620215/149‑1 (MQ=255)
                      ccGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCt                                                                                                                                 <  1:298543/150‑1 (MQ=255)
                      ccGCTGCCGGGTTAGCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCt                                                                                                                                 <  1:536699/150‑1 (MQ=255)
                         cTGCCGGGGAATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCt                                                                                                                              <  2:224374/150‑1 (MQ=255)
                             cGGGTTATCGATAAGCTGTACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCaaa                                                                                                                          <  1:1694365/150‑1 (MQ=255)
                                gTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGc                                                                                                                       <  1:285874/150‑1 (MQ=255)
                                  tATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCg                                                                                                                     <  1:1211340/150‑1 (MQ=255)
                                  tATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCg                                                                                                                     <  2:1676621/150‑1 (MQ=255)
                                  tATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCg                                                                                                                     <  1:1226272/150‑1 (MQ=255)
                                     cGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGAt                                                                                                                  >  2:154289/1‑150 (MQ=255)
                                      gAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATg                                                                                                                 <  2:1441449/150‑1 (MQ=255)
                                      gAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATg                                                                                                                 >  1:875069/1‑150 (MQ=255)
                                          aGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGAAAGATCGACGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATgcggc                                                                                                             <  2:1171827/150‑1 (MQ=255)
                                           gCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATgcggcg                                                                                                            <  1:432134/150‑1 (MQ=255)
                                            cTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGt                                                                                                           >  1:1469678/1‑150 (MQ=255)
                                              ggACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTaa                                                                                                         >  1:1710120/1‑150 (MQ=255)
                                                   cgGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGcc                                                                                                    <  2:579678/150‑1 (MQ=255)
                                                      aCGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTa                                                                                                 <  1:154289/150‑1 (MQ=255)
                                                       cGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTAt                                                                                                <  2:1320026/150‑1 (MQ=255)
                                                       cGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTCTGGGGAAAATAACCCCCGGAGTGTAGGCCAGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTAt                                                                                                >  2:101345/1‑150 (MQ=255)
                                                        gATGAAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATATGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATc                                                                                               <  1:1076418/150‑1 (MQ=255)
                                                         aTGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcc                                                                                              <  2:480138/150‑1 (MQ=255)
                                                         aTGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcc                                                                                              >  1:1121442/1‑150 (MQ=255)
                                                          tGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAACTGCCGGATGCGGCGTAAACGCCTTATcct                                                                                             >  1:1484941/1‑150 (MQ=255)
                                                            cAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCCTCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcctgc                                                                                           >  1:684870/1‑150 (MQ=255)
                                                             aaaaGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGTGGGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcctgcc                                                                                          <  2:970393/150‑1 (MQ=255)
                                                             aaaaGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcctgcc                                                                                          <  2:875069/150‑1 (MQ=255)
                                                                aGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCctgcctgc                                                                                       <  1:827767/150‑1 (MQ=255)
                                                                 gcgcCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCAGGATAAGGCGTTTACGCCGCATccggcatttg                                                                                      <  2:673382/150‑6 (MQ=38)
                                                                   gcCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCaaa                                                                                    >  1:1583887/1‑150 (MQ=255)
                                                                         tGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGaa                                                                              >  1:911848/1‑150 (MQ=255)
                                                                              ggCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTa                                                                         <  2:1583887/150‑1 (MQ=255)
                                                                                cAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAAc                                                                       <  2:1150882/150‑1 (MQ=255)
                                                                                  aGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGaa                                                                              >  1:510627/1‑141 (MQ=255)
                                                                                  aGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGaa                                                                              <  2:510627/141‑1 (MQ=255)
                                                                                  aGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTc                                                                     >  1:1594102/1‑150 (MQ=255)
                                                                                    aTCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCg                                                                   >  1:1351075/1‑150 (MQ=255)
                                                                                     tCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCTGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGc                                                                  <  1:252451/150‑1 (MQ=255)
                                                                                      cGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTTAACGCCTTATCCTGCCTGCAAATGCGATGTTTAACTCCGCt                                                                 >  1:1630476/1‑150 (MQ=255)
                                                                                      cGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCt                                                                 >  1:1608197/1‑150 (MQ=255)
                                                                                      cGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCt                                                                 >  2:530912/1‑150 (MQ=255)
                                                                                       gAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAACCGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCtt                                                                >  1:1039738/1‑150 (MQ=255)
                                                                                       gAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCtt                                                                >  2:414505/1‑150 (MQ=255)
                                                                                        aaGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAACCGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTta                                                               >  1:655515/1‑150 (MQ=255)
                                                                                          gATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTtata                                                             >  1:1687845/1‑150 (MQ=255)
                                                                                           aTTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATAc                                                            >  2:273773/1‑150 (MQ=255)
                                                                                           aTTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATAc                                                            <  1:476882/150‑1 (MQ=255)
                                                                                            tttttCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACa                                                           <  1:602656/150‑1 (MQ=255)
                                                                                               ttCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCg                                                        >  2:183541/1‑150 (MQ=255)
                                                                                                tCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCgg                                                       <  2:666288/150‑1 (MQ=255)
                                                                                                tCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCgg                                                       >  1:367754/1‑150 (MQ=255)
                                                                                                 cccGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGTATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGc                                                      >  2:1184587/1‑150 (MQ=255)
                                                                                                  ccgccgGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAATTGCGAAGTGTAACTCCGCTTATACAGCGGCa                                                     >  2:261862/1‑150 (MQ=255)
                                                                                                  ccgccgGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCa                                                     >  1:1409822/1‑150 (MQ=255)
                                                                                                  ccgccgGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCa                                                     >  1:1416718/1‑150 (MQ=255)
                                                                                                   cgccgGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCaa                                                    >  1:175407/1‑150 (MQ=255)
                                                                                                      cgGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCa                                                 <  2:520991/150‑1 (MQ=255)
                                                                                                       ggTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCaa                                                <  2:367754/150‑1 (MQ=255)
                                                                                                        gTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCaaa                                               >  1:525334/1‑150 (MQ=255)
                                                                                                          ggtgGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTATCTGCGCTTATACAGCGGCAACCAAATc                                             <  2:1563100/149‑1 (MQ=255)
                                                                                                           gtgGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATcc                                            >  2:1041179/1‑149 (MQ=255)
                                                                                                           gtgGGGAAAAAACCCCCCGGAGTGTCGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTAACCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATCCGTCGGCAACCAAATCa                                            >  2:1041183/1‑150 (MQ=37)
                                                                                                            tgGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCAc                                           >  2:440272/1‑150 (MQ=255)
                                                                                                            tgGGGAAAATAACACCCGGAGTGTAGGCATGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCTGCAACCAAATCAc                                           >  2:15917/1‑150 (MQ=255)
                                                                                                             ggggAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCAcc                                          <  2:162439/150‑1 (MQ=255)
                                                                                                               ggAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAg                                        >  1:1019608/1‑150 (MQ=255)
                                                                                                               ggAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAg                                        >  1:1690975/1‑150 (MQ=255)
                                                                                                                  aaaTAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCg                                     >  2:339052/1‑150 (MQ=255)
                                                                                                                        cccccGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAAcc                               >  1:257134/1‑150 (MQ=255)
                                                                                                                         ccccGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCg                              >  2:1451005/1‑150 (MQ=255)
                                                                                                                           ccGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGAGTAAACGGCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGcc                            <  2:252706/150‑1 (MQ=255)
                                                                                                                             ggAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCa                          >  1:999160/1‑150 (MQ=255)
                                                                                                                                 tgtAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGTgg                      >  1:386642/1‑150 (MQ=255)
                                                                                                                                 tgtAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCgg                      >  2:985350/1‑150 (MQ=255)
                                                                                                                                 tgtAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCgg                      >  2:489154/1‑150 (MQ=255)
                                                                                                                                 tgtAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCgg                      <  1:414505/150‑1 (MQ=255)
                                                                                                                                  ggatgCCTGATGCGACGCTGATAGCGTCTTATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCggg                     <  2:704331/146‑1 (MQ=38)
                                                                                                                                  ggatgCCTGATGCGACGCTGATAGCGTCTTATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCggg                     >  1:704331/5‑150 (MQ=38)
                                                                                                                                     ggCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCACTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTg                  >  2:1634254/1‑150 (MQ=255)
                                                                                                                                       ccTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtc                >  1:886892/1‑150 (MQ=255)
                                                                                                                                        cTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtct               <  2:1594102/150‑1 (MQ=255)
                                                                                                                                        cTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtct               >  2:865815/1‑150 (MQ=255)
                                                                                                                                          gATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtcttc             <  2:1416718/150‑1 (MQ=255)
                                                                                                                                          gATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtcttc             >  1:289456/1‑150 (MQ=255)
                                                                                                                                          gATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtcttc             <  1:1451005/150‑1 (MQ=255)
                                                                                                                                            tAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTg           <  2:1598729/150‑1 (MQ=255)
                                                                                                                                               gACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCt        >  2:923155/1‑150 (MQ=255)
                                                                                                                                               gACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCt        <  2:971775/150‑1 (MQ=255)
                                                                                                                                                     aTCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCAccc  >  2:1024963/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCC  >  NC_000913/4103377‑4103675

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