Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,103,527 T→A 5.2% intergenic (+31/+75) 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,5270TA 5.2% 395.3 / 5.3 131intergenic (+31/+75)yiiM/cpxA6‑N‑hydroxylaminopurine resistance protein/sensory histidine kinase in two‑component regulatory system with CpxR
Reads supporting (aligned to +/- strand):  ref base T (53/69);  new base A (3/5);  total (56/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

ATCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACC  >  NC_000913/4103378‑4103674
                                                                                                                                                     |                                                                                                                                                   
aTCAATAGCACCGATTACTCTCCGCTGCCGGGTTATCGAACAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTAt                                                                                                                                                     <  2:1264914/150‑1 (MQ=255)
 tCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGATGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCAATc                                                                                                                                                    >  1:992962/1‑150 (MQ=255)
  cAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCa                                                                                                                                                   >  2:1187653/1‑150 (MQ=255)
    atatCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGc                                                                                                                                                 <  2:417628/150‑1 (MQ=255)
       tCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTc                                                                                                                                              <  1:1267229/150‑1 (MQ=255)
       tCACCGTTTACTCTCCGCTGCCGGGGTATCGAAAAGCTGGACGCGTAAGAGGCAAAAGAGACGACGGAGCGGCTAGATCGAATCTTTTTCAAGCCGGTTTTGGGTAAAATAACCCACGGATTTTAGGCCTGATAATACGCTATCAGCGTc                                                                                                                                              <  1:1453041/150‑1 (MQ=255)
         aCCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGc                                                                                                                                            <  2:534287/150‑1 (MQ=255)
          ccGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCa                                                                                                                                           >  1:551771/1‑150 (MQ=255)
           cgattACTCTCCGCTGCCGGGTTATCGTAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCAt                                                                                                                                          <  2:1481625/147‑1 (MQ=255)
              ttACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAg                                                                                                                                       <  1:734786/150‑1 (MQ=255)
               tACTCTCCGCTGCCGGGGTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATGCGACGCTGATAGCGTCTTATCAgg                                                                                                                                      <  1:1007907/150‑1 (MQ=255)
               tACTCTCCGCGGCCGGGTTAGCGAAAAGCGGGAATCGGAAGATGCAAAAGCGCCGACTGAGCGGCAAGATCTAAGATTTTTCCCGCCGGCTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAgg                                                                                                                                      <  1:174113/150‑1 (MQ=255)
               tACTCTCCGCAGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCTCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATGCGACGCTGATAGCGTCTTATCAgg                                                                                                                                      <  1:575554/150‑1 (MQ=255)
                aCTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAATGCTCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATGCGACGCTGATAGCGTCTTATCAGGc                                                                                                                                     <  2:281672/150‑1 (MQ=255)
                 ttctcCTCTGCAGTGTTTTCTATAAGTTGTAGGCTGACGATGCAAATGCGCCGACTGTGCGGCAAGATCGAAGTTTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATGCGACGCTGATAGCGTCTTATCAGGcc                                                                                                                                    <  1:1807843/149‑2 (MQ=255)
                 ctctCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAATGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCa                                                                                                                                    >  2:1473964/1‑150 (MQ=255)
                  tctcCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCAt                                                                                                                                   >  1:1525963/1‑150 (MQ=255)
                    tcCGCTACCGGGTTATCGAAAAGCTGGACTCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATcc                                                                                                                                 <  2:701363/150‑1 (MQ=255)
                        cTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCt                                                                                                                             >  2:1144690/1‑150 (MQ=255)
                           ccGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCaa                                                                                                                          >  2:1732405/1‑150 (MQ=255)
                            cGGGTTATCGAAAAGCTGGACGAGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTAACAGCGTCGCATCAGGCATCCTGCTCaaa                                                                                                                         >  1:356553/1‑150 (MQ=255)
                              ggTTATGCAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTTAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATg                                                                                                                       <  1:882100/150‑1 (MQ=255)
                                ttATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGcc                                                                                                                     <  1:1550790/150‑1 (MQ=255)
                                 tATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCg                                                                                                                    <  2:427687/150‑1 (MQ=255)
                                 tATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCg                                                                                                                    >  2:38644/1‑150 (MQ=255)
                                    cGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGAt                                                                                                                 >  1:457134/1‑150 (MQ=255)
                                     gAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATg                                                                                                                <  2:188373/150‑1 (MQ=255)
                                     gAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATg                                                                                                                >  2:1275609/1‑150 (MQ=255)
                                          gCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCCTCCTGCTCAAATGCCGGATgcggcg                                                                                                           >  1:104469/1‑150 (MQ=255)
                                          gCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATgcggcg                                                                                                           >  2:1156264/1‑150 (MQ=255)
                                           cTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGt                                                                                                          >  1:530975/1‑150 (MQ=255)
                                            tGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTa                                                                                                         <  2:1697207/150‑1 (MQ=255)
                                               ccgcgGACGATGCAAAAGCGCCGACTGAGCGGCATGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGGGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAAc                                                                                                      <  2:897627/149‑1 (MQ=255)
                                                  cgGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGGATCCTGCTCAAATGCCGGATGCGGCGTAAACGcc                                                                                                   >  2:269790/1‑150 (MQ=255)
                                                  cgGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGcc                                                                                                   >  2:1550272/1‑150 (MQ=255)
                                                   ggACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCt                                                                                                  <  1:553123/150‑1 (MQ=255)
                                                   ggACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCt                                                                                                  <  1:1572302/150‑1 (MQ=255)
                                                    gACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCtt                                                                                                 <  2:1542890/150‑1 (MQ=255)
                                                      cGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCCTCAGGCAGCCTGCTCAAATGCCGGATGCGGCGTAACCGCCGTAt                                                                                               >  1:673639/1‑150 (MQ=255)
                                                        aTGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcc                                                                                             <  2:6878/150‑1 (MQ=255)
                                                          gCAAACGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcctg                                                                                           <  2:1503449/150‑1 (MQ=255)
                                                           cAAAAGCGCCGCCTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATACCCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcctgc                                                                                          >  1:1190303/1‑150 (MQ=255)
                                                           cAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCATAAACGCCTTATcctgc                                                                                          <  1:832728/150‑1 (MQ=255)
                                                            aaaaGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcctgcc                                                                                         >  1:1189070/1‑150 (MQ=255)
                                                             aaaGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATcctgcct                                                                                        <  2:1805424/150‑1 (MQ=255)
                                                                    cGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATg                                                                                 <  1:1628995/150‑1 (MQ=255)
                                                                      aCTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCg                                                                               <  2:1170140/150‑1 (MQ=255)
                                                                       cTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGa                                                                              >  1:746546/1‑150 (MQ=255)
                                                                        tGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGaa                                                                             <  1:426273/150‑1 (MQ=255)
                                                                          aGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGt                                                                           >  1:1031179/1‑150 (MQ=255)
                                                                            cGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGttt                                                                         >  1:593385/1‑150 (MQ=255)
                                                                             ggCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTa                                                                        >  1:1111792/1‑150 (MQ=255)
                                                                             ggCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTa                                                                        <  1:1045101/150‑1 (MQ=255)
                                                                              gCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTaa                                                                       >  1:428226/1‑150 (MQ=255)
                                                                               cAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGACAATACCCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCCAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAAc                                                                      >  1:29183/1‑150 (MQ=255)
                                                                               cAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAAc                                                                      <  1:986117/150‑1 (MQ=255)
                                                                               cAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAAc                                                                      <  1:1810004/150‑1 (MQ=255)
                                                                               cAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAAc                                                                      <  1:1106814/150‑1 (MQ=255)
                                                                                    tCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGc                                                                 <  1:927685/150‑1 (MQ=255)
                                                                                     cGAAGATTTTTCCCGCCGGTTGTGGTGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCt                                                                <  1:850874/150‑1 (MQ=255)
                                                                                     cGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAACCGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCt                                                                >  1:247066/1‑150 (MQ=255)
                                                                                     cGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCt                                                                >  2:892658/1‑150 (MQ=255)
                                                                                      gAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCtt                                                               <  1:1762887/150‑1 (MQ=255)
                                                                                         gATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCGTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTtata                                                            <  1:741291/150‑1 (MQ=255)
                                                                                           tttttCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACa                                                          >  1:961028/1‑150 (MQ=255)
                                                                                               tCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCgg                                                      <  2:120217/150‑1 (MQ=255)
                                                                                               tCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCgg                                                      <  2:595599/150‑1 (MQ=255)
                                                                                               tCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCgg                                                      <  1:648233/150‑1 (MQ=255)
                                                                                                 ccgcctGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTTCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCa                                                    <  1:740746/150‑1 (MQ=255)
                                                                                                          gtgGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCa                                           <  2:330434/150‑1 (MQ=255)
                                                                                                          gtgGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCa                                           <  1:1431524/150‑1 (MQ=255)
                                                                                                            ggggAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCAcc                                         >  1:1426286/1‑150 (MQ=255)
                                                                                                             gggAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCa                                        <  2:1705571/150‑1 (MQ=255)
                                                                                                              ggAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCac                                       >  1:1638395/1‑149 (MQ=255)
                                                                                                              ggAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAg                                       <  1:1001072/150‑1 (MQ=255)
                                                                                                               gAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGc                                      <  1:1245823/150‑1 (MQ=255)
                                                                                                                 aaaTAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCg                                    <  1:587931/150‑1 (MQ=255)
                                                                                                                   aTAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCCCCAGCCGTa                                  >  1:1338342/1‑150 (MQ=255)
                                                                                                                      aCCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAAc                               >  2:1471024/1‑150 (MQ=255)
                                                                                                                            ggAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCa                         <  2:495993/150‑1 (MQ=255)
                                                                                                                               gtgtAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTATACGCCTTATCCTGCCTGCATATGCGAAGGTTTCCTCCGCTTAGTCAGCGGCAACCAGAGCACCAGCCGTAATCCGCCCcgcg                      >  2:906055/1‑150 (MQ=25)
                                                                                                                               caggatgCCTGATGCGACGCTGATAGCGTCTTATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTAAACTCCGCTTATACAGCGGCAACAAAATCACCAGCCGTAAACCGCCCAGCg                      >  2:1000758/7‑150 (MQ=37)
                                                                                                                                  tAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAACTCACCAGCCGTAAACCGCCCAGCGGGc                   <  2:540889/150‑1 (MQ=255)
                                                                                                                                  tAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGc                   <  1:174357/150‑1 (MQ=255)
                                                                                                                                   aGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCt                  >  2:129141/1‑150 (MQ=255)
                                                                                                                                    ggCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTg                 <  1:1156264/150‑1 (MQ=255)
                                                                                                                                    ggCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTg                 <  1:1626975/150‑1 (MQ=255)
                                                                                                                                      ccTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAGCTCCGCTTATACAGCGTCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTTtc               <  2:1031179/150‑1 (MQ=255)
                                                                                                                                       cTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtct              <  2:1785657/150‑1 (MQ=255)
                                                                                                                                        tGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGATGTGTAACTCCGCTTATACAGCGGCAACCAACTCACCAGCCGTAAACCGCCCAGCGGGCTGtctt             >  2:1107554/1‑150 (MQ=255)
                                                                                                                                         gATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAGTGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtcttc            >  1:109230/1‑150 (MQ=255)
                                                                                                                                         gATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGtcttc            <  1:1557501/150‑1 (MQ=255)
                                                                                                                                           tAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTg          <  1:1068409/150‑1 (MQ=255)
                                                                                                                                           aaaGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTg          <  1:309170/149‑1 (MQ=255)
                                                                                                                                             aGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGcc        <  2:1445346/150‑1 (MQ=255)
                                                                                                                                              gACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCt       >  1:396808/1‑150 (MQ=255)
                                                                                                                                                cGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTc     >  1:643797/1‑150 (MQ=255)
                                                                                                                                                 gCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCa    <  1:1421892/150‑1 (MQ=255)
                                                                                                                                                  cTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCAc   <  1:1005212/150‑1 (MQ=255)
                                                                                                                                                   tATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCAcc  >  1:207941/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACC  >  NC_000913/4103378‑4103674

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