Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 747,902 C→G 7.6% intergenic (‑133/+19) abrB ← / ← ybgO regulator of aidB expression; inner membrane protein/putative fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913747,9020CG 7.6% 385.1 / 19.4 127intergenic (‑133/+19)abrB/ybgOregulator of aidB expression; inner membrane protein/putative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base C (63/54);  new base G (4/5);  total (67/60)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-02

Only 100 of 134 total aligned reads displayed.

CACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGT  >  NC_000913/747753‑748051
                                                                                                                                                     |                                                                                                                                                     
   tGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATATCAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCaaa                                                                                                                                                    >  2:1767269/1‑150 (MQ=255)
    gCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAg                                                                                                                                                   <  2:1707151/150‑1 (MQ=255)
    gCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGACCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAg                                                                                                                                                   >  2:1449571/1‑150 (MQ=255)
     cAAAACTTGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAgg                                                                                                                                                  <  2:1650560/150‑1 (MQ=255)
      aaaaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAggg                                                                                                                                                 >  2:661837/1‑150 (MQ=255)
          cTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGcatc                                                                                                                                             <  1:772655/150‑1 (MQ=255)
            ggCATCCTTCGCACCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGcatctc                                                                                                                                           <  1:684316/150‑3 (MQ=255)
             gCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATg                                                                                                                                          <  1:1842028/150‑1 (MQ=255)
              cATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAACTCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATgg                                                                                                                                         <  1:1739922/150‑1 (MQ=255)
              cATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATgg                                                                                                                                         <  2:1419815/150‑1 (MQ=255)
              cATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATgg                                                                                                                                         >  1:1477268/1‑150 (MQ=255)
                tCCTTCGCTCCCTCATTATTTACTCTTTGGTGTATCGATACAGCATAAAGACAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTa                                                                                                                                       <  1:265758/150‑1 (MQ=255)
                 ccTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGcatctctttca                                                                                                                                      <  1:1727703/150‑8 (MQ=255)
                  cTTCTCTACCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTaaa                                                                                                                                     <  2:1919287/150‑1 (MQ=255)
                    tCGCTCCCTCATTATTTTCTCTTTCGAGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATATTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATc                                                                                                                                   <  2:283818/150‑1 (MQ=255)
                       cTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                                                >  1:1280217/1‑150 (MQ=255)
                       cTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                                                <  1:212162/150‑1 (MQ=255)
                        tCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAATTCAGtt                                                                                                                               <  1:1120481/150‑1 (MQ=255)
                         cccTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAAAAACCGCAATCAAGGAAAAAGGGAAAATCAGCAATTTTGAGAAAGAAAATCCCGTTGCGGCGGCAAAGGGCAACATGGTAAATAAGttt                                                                                                                              >  2:1977050/1‑149 (MQ=255)
                          cctacttattTTCTCTTTGGTGTATCGATACAGCATACAGAACCATGGTCCATGAAAAAATAACACAAATCAAGGAAACTGGGAAATTCAGCAATTTTCTGAAATAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTtat                                                                                                                             <  2:964803/145‑1 (MQ=255)
                           cTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                            <  2:1038655/150‑1 (MQ=255)
                              attattTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCTAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGtt                                                                                                                         >  1:1797680/1‑150 (MQ=255)
                                  ttttCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                                     >  1:1009831/1‑150 (MQ=255)
                                        ctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGCTGCACTTTCCGGCGGCTAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                               >  1:129745/1‑150 (MQ=255)
                                        ctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                               >  2:1729026/1‑150 (MQ=255)
                                        ctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                               <  1:1113804/150‑1 (MQ=255)
                                           tGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAACCAGTTATAGTTAATCTTAAAAAGc                                                                                                            >  1:1974773/1‑150 (MQ=255)
                                            ggTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                           >  1:992851/1‑150 (MQ=255)
                                            ggTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                           <  1:1382280/150‑1 (MQ=255)
                                              tgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcacc                                                                                                         <  1:563748/150‑1 (MQ=255)
                                                tATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcaccac                                                                                                       >  2:1269890/1‑150 (MQ=255)
                                                   cGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGt                                                                                                    >  1:1578051/1‑150 (MQ=255)
                                                    gATACAGCATAAAGAAAGCTGGGCCATGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAATGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTa                                                                                                   <  1:1901054/150‑1 (MQ=255)
                                                    gATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTa                                                                                                   >  1:1097893/1‑150 (MQ=255)
                                                     aTACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTAt                                                                                                  >  1:1732722/1‑150 (MQ=255)
                                                      tACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTAta                                                                                                 >  1:823914/1‑149 (MQ=255)
                                                        cAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCATTTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATcaa                                                                                               >  1:1774279/1‑148 (MQ=255)
                                                        cAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGa                                                                                               >  1:318893/1‑150 (MQ=255)
                                                           cATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAgg                                                                                            >  1:493081/1‑150 (MQ=255)
                                                            aTAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGAATCGAAggg                                                                                           >  2:402259/1‑150 (MQ=255)
                                                             tAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGATAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCTAAGGGCATCATGGTAAATCAGTTATAGTTAAACTTAAAAAGCACCACTGTATCGAAgcgg                                                                                          >  1:294343/1‑150 (MQ=255)
                                                              aaagaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCATTTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGc                                                                                         >  1:30950/1‑150 (MQ=255)
                                                                agaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACAGTATCGAAGGGGCCa                                                                                       >  2:558458/1‑150 (MQ=255)
                                                                 gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAg                                                                                      >  2:920214/1‑150 (MQ=255)
                                                                 gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAg                                                                                      >  1:1930372/1‑150 (MQ=255)
                                                                 gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAg                                                                                      >  1:1872354/1‑150 (MQ=255)
                                                                     gCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTc                                                                                  >  1:527279/1‑150 (MQ=255)
                                                                      cTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCg                                                                                 <  1:1511262/150‑1 (MQ=255)
                                                                        gggCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGAt                                                                               >  2:1929627/1‑150 (MQ=255)
                                                                        gggCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGAt                                                                               >  2:1387903/1‑150 (MQ=255)
                                                                         ggCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTAGAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATg                                                                              >  1:139667/1‑150 (MQ=255)
                                                                              ggAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                         >  1:1559630/1‑150 (MQ=255)
                                                                              ggAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                         >  2:1619800/1‑150 (MQ=255)
                                                                                 aaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGcc                                                                      <  2:57608/150‑1 (MQ=255)
                                                                                 aaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGcc                                                                      >  1:1413854/1‑150 (MQ=255)
                                                                                 aaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGcc                                                                      <  1:721642/150‑1 (MQ=255)
                                                                                   aaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                                    <  1:570650/150‑1 (MQ=255)
                                                                                   aaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAAATTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                                    <  1:683996/150‑1 (MQ=255)
                                                                                    aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTgg                                                                   >  1:1686106/1‑150 (MQ=255)
                                                                                        aCAGCAATCCAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTATCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAACAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATg                                                               <  1:1845256/150‑1 (MQ=255)
                                                                                        aCAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATAATGGTAATTCAGTTATAGATAATCTTAAAACGCAACACTGTATCGAAGTGGAAAGCATAGATGGGGGTTCATTGAATg                                                               >  1:87076/1‑150 (MQ=255)
                                                                                           gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCGCGATGGGGGTGCCTGGAATGGAc                                                            <  1:1836458/150‑1 (MQ=255)
                                                                                             aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGa                                                          >  2:1712774/1‑150 (MQ=255)
                                                                                              aTCCAGTAAAAGTGGTATATCAGCAATAATCTGAAAGAGATTCCCTTTCCGGCTGTTAAGGTCTTCATTGTAACTCAGTTATAGTTAATCTTAAAAAGCCCCACTGTATCGAAGGGGCCAGCCGCGATGGGGGGGCCTGGAATGGACGaa                                                         <  2:1269103/150‑1 (MQ=255)
                                                                                                   ggAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTc                                                    <  2:1923628/150‑1 (MQ=255)
                                                                                                   ggAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCAATGAATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTc                                                    >  2:120375/1‑150 (MQ=255)
                                                                                                   ggAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAACCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTc                                                    <  2:1009831/150‑1 (MQ=255)
                                                                                                     aaaaGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAg                                                  <  1:1381998/150‑1 (MQ=255)
                                                                                                        aGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCg                                               <  1:1441083/150‑1 (MQ=255)
                                                                                                        aGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCg                                               >  2:1712681/1‑150 (MQ=255)
                                                                                                         ggggAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGt                                              >  1:912174/1‑150 (MQ=255)
                                                                                                          gggAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCCGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGta                                             >  2:781294/1‑150 (MQ=255)
                                                                                                           ggAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGTCCCGCCTCCATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtat                                            >  2:299481/1‑150 (MQ=255)
                                                                                                           ggAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtat                                            <  1:1793865/150‑1 (MQ=255)
                                                                                                           ggAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtat                                            >  2:1227637/1‑150 (MQ=255)
                                                                                                           ggAAAATCAGCAATTATCTGAAAGAGATGCCCTTTCCGGCGGCATAGGTCATCATGGAAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGTTGCCTGGAATGGACGAAAGTTCAGCCGtat                                            <  2:106424/150‑1 (MQ=255)
                                                                                                              aaaTCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCTAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAgg                                         <  2:1961086/150‑1 (MQ=255)
                                                                                                              aaaTCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAgg                                         <  1:1148651/150‑1 (MQ=255)
                                                                                                              aaaTCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAgg                                         <  2:493081/150‑1 (MQ=255)
                                                                                                                    gCAATTTTCTGAAAGAGATGCCCTGTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAACGCACCCATGTATCGAAGGGGCCAGACTCGCGTGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTAAAGGtttgtt                                   >  1:367408/1‑150 (MQ=255)
                                                                                                                     cAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCAGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTc                                  >  2:1517575/1‑150 (MQ=255)
                                                                                                                     cAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGGGCCTGGAATGAACGAAAGTTCAGCCGTATAGGTTTTTTc                                  >  1:1020591/1‑150 (MQ=255)
                                                                                                                       aTTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTTTAGTTAATCTTAAAAAGCACCACTGTATCGACGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGt                                <  1:1269890/150‑1 (MQ=255)
                                                                                                                       aTTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGt                                <  2:762915/150‑1 (MQ=255)
                                                                                                                        ttttCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTa                               <  2:1732722/150‑1 (MQ=255)
                                                                                                                            cTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGAcc                           >  1:1664988/1‑150 (MQ=255)
                                                                                                                              gAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCaa                         >  2:1025879/1‑150 (MQ=255)
                                                                                                                              gAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCaa                         >  1:635461/1‑150 (MQ=255)
                                                                                                                                aaGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACa                       <  2:477474/150‑1 (MQ=255)
                                                                                                                                 agagATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACAt                      <  2:1991331/150‑1 (MQ=255)
                                                                                                                                  gagaAGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATg                     <  2:1413854/150‑1 (MQ=255)
                                                                                                                                  atgaTGCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATg                     >  2:1022404/3‑150 (MQ=255)
                                                                                                                                    gaTGCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTc                   >  2:711859/1‑150 (MQ=255)
                                                                                                                                     aTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCa                  >  2:764939/1‑150 (MQ=255)
                                                                                                                                      tGCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCAc                 >  1:1340022/1‑150 (MQ=255)
                                                                                                                                            ttCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGa           <  2:1572381/150‑1 (MQ=255)
                                                                                                                                                ggcggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCCGCCTCGAAGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATc       >  2:1269435/1‑150 (MQ=255)
                                                                                                                                                    gcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCg   >  2:1645786/1‑150 (MQ=255)
                                                                                                                                                    gcAAAGGGCATCATGGTAAATCAGATATAGTTAATCTTAAAAAGCACCACTGTCTGGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCg   <  2:1247265/150‑1 (MQ=255)
                                                                                                                                                     cAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGt  >  1:537864/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGT  >  NC_000913/747753‑748051

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