Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,103,528 C→T 5.1% intergenic (+32/+74) 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,5280CT 5.1% 386.5 / 4.9 117intergenic (+32/+74)yiiM/cpxA6‑N‑hydroxylaminopurine resistance protein/sensory histidine kinase in two‑component regulatory system with CpxR
Reads supporting (aligned to +/- strand):  ref base C (60/51);  new base T (2/4);  total (62/55)
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.62e-01

Only 100 of 119 total aligned reads displayed.

TCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAG  >  NC_000913/4103379‑4103677
                                                                                                                                                     |                                                                                                                                                     
tCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACACTATc                                                                                                                                                       <  2:445492/150‑1 (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:742476/150‑1 (MQ=255)
              tACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAgg                                                                                                                                         <  2:1592532/150‑1 (MQ=255)
              tACTCTCCGCGGCCGGGGTATCGAAAAGCTGGACGCGGACGAGGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATGCGACGCTGATAGCGTCTTATCAgg                                                                                                                                         <  1:692076/150‑1 (MQ=255)
               aCTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGc                                                                                                                                        >  1:1057632/1‑150 (MQ=255)
               aCTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGACGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCTCATCAGGc                                                                                                                                        >  2:1557940/1‑150 (MQ=255)
                ctctCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCa                                                                                                                                       >  2:827767/1‑150 (MQ=255)
                  gtcCGCTGCCGGGTTATCGAAATGGTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATc                                                                                                                                     <  2:620215/149‑1 (MQ=255)
                   tcCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGGAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATcc                                                                                                                                    <  2:558888/150‑1 (MQ=255)
                   tcCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATcc                                                                                                                                    >  1:174422/1‑150 (MQ=255)
                    ccGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCt                                                                                                                                   <  1:298543/150‑1 (MQ=255)
                    ccGCTGCCGGGTTAGCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCt                                                                                                                                   <  1:536699/150‑1 (MQ=255)
                       cTGCCGGGGAATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCt                                                                                                                                <  2:224374/150‑1 (MQ=255)
                         gCCGGGTTAGCGACCTGCTGGACGCGGAGGATGCACAAGCGCCGCCTGAGCGGCAAGTTGGATGATTTTTCCCGCCGGGGCTGGGGATAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCa                                                                                                                              <  2:302485/150‑1 (MQ=255)
                           cGGGTTATCGATAAGCTGTACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCaaa                                                                                                                            <  1:1694365/150‑1 (MQ=255)
                              gTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGc                                                                                                                         <  1:285874/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                                                                                                                   >  1:875069/1‑150 (MQ=255)
                                    gAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATg                                                                                                                   <  2:1441449/150‑1 (MQ=255)
                                       aaGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATgcgg                                                                                                                <  2:1407758/150‑1 (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)
                                                  ggACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCt                                                                                                     <  1:1188275/150‑1 (MQ=255)
                                                    aCGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTa                                                                                                   <  1:154289/150‑1 (MQ=255)
                                                     cGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCCTCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTAt                                                                                                  >  2:1100241/1‑150 (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)
                                                            aaaGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcctgcct                                                                                           >  1:45184/1‑150 (MQ=255)
                                                              aGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCctgcctgc                                                                                         <  1:827767/150‑1 (MQ=255)
                                                                 gcCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCaaa                                                                                      >  1:1583887/1‑150 (MQ=255)
                                                                       tGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGaa                                                                                >  1:911848/1‑150 (MQ=255)
                                                                                aGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGaa                                                                                <  2:510627/141‑1 (MQ=255)
                                                                                aGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGaa                                                                                >  1:510627/1‑141 (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)
                                                                                     gAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAACCGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCtt                                                                  >  1:1039738/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)
                                                                                               cccGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGc                                                        >  1:1432661/1‑150 (MQ=255)
                                                                                                ccgccgGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAATTGCGAAGTGTAACTCCGCTTATACAGCGGCa                                                       >  2:261862/1‑150 (MQ=255)
                                                                                                ccgccgGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCa                                                       >  1:1416718/1‑150 (MQ=255)
                                                                                                ccgccgGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCa                                                       <  1:567971/150‑1 (MQ=255)
                                                                                                ccgccgGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCa                                                       >  1:1409822/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)
                                                                                                       ttGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAAt                                                <  1:1641920/150‑1 (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)
                                                                                                            gggAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCa                                           >  2:1596091/1‑150 (MQ=255)
                                                                                                             ggAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAg                                          >  1:1019608/1‑150 (MQ=255)
                                                                                                              gAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGc                                         <  2:1710120/150‑1 (MQ=255)
                                                                                                                aaaTAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCg                                       >  2:339052/1‑150 (MQ=255)
                                                                                                                      cccccGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAAcc                                 >  1:965332/1‑150 (MQ=255)
                                                                                                                       ccccGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCg                                >  2:1451005/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)
                                                                                                                                 tAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGc                      <  1:1590599/150‑1 (MQ=255)
                                                                                                                                   ggCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCACTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTg                    >  2:1634254/1‑150 (MQ=255)
                                                                                                                                     ccTGATGCGACGCTGATAGCGTCTTATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtc                  >  2:907039/1‑150 (MQ=38)
                                                                                                                                     ccTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtc                  >  1:886892/1‑150 (MQ=255)
                                                                                                                                      cTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtct                 >  2:865815/1‑150 (MQ=255)
                                                                                                                                      cTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtct                 <  2:1594102/150‑1 (MQ=255)
                                                                                                                                        gATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtcttc               <  2:676851/150‑1 (MQ=255)
                                                                                                                                        gATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtcttc               >  1:289456/1‑150 (MQ=255)
                                                                                                                                         aTAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtcttct              >  1:29086/1‑150 (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)
                                                                                                                                                     cAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAg  >  2:1458014/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAG  >  NC_000913/4103379‑4103677

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: