Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 747,897 G→C 7.1% intergenic (‑128/+24) 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,8970GC 7.1% 388.4 / 18.9 127intergenic (‑128/+24)abrB/ybgOregulator of aidB expression; inner membrane protein/putative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (63/55);  new base C (4/5);  total (67/60)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

Only 100 of 134 total aligned reads displayed.

CCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATC  >  NC_000913/747750‑748046
                                                                                                                                                   |                                                                                                                                                     
ccccACTCCAACCCTCGCATCCAGCGCTCCCACAGAATCTACACATCGGCGGCTCGCTCCAGCAAAACGCAACCTGGCCCAGGAAAAAATAACAGCAATCAAGGAAAAGCGGAAACTCACCAAATTTCTGAAAGAGCTGCCCTCTCcggc                                                                                                                                                     <  2:1576696/150‑1 (MQ=255)
  ccACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCcggcgg                                                                                                                                                   <  2:160862/150‑1 (MQ=255)
      tGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATATCAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCaaa                                                                                                                                               >  2:1767269/1‑150 (MQ=255)
       gCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAg                                                                                                                                              <  2:1707151/150‑1 (MQ=255)
        cAAAACTTGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAgg                                                                                                                                             <  2:1650560/150‑1 (MQ=255)
        cAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                                             >  1:1812673/1‑150 (MQ=255)
         aaaaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAggg                                                                                                                                            >  2:661837/1‑150 (MQ=255)
             cTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGcatc                                                                                                                                        <  1:1577738/150‑1 (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)
                   tCCTTCGCTCCCTCATTATTTACTCTTTGGTGTATCGATACAGCATAAAGACAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTa                                                                                                                                  <  1:265758/150‑1 (MQ=255)
                     cTTCTCTACCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTaaa                                                                                                                                <  2:1919287/150‑1 (MQ=255)
                     cTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTaaa                                                                                                                                >  2:1325666/1‑150 (MQ=255)
                       tCGCTCCCTCATTATTTTCTCTTTCGAGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATATTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATc                                                                                                                              <  2:283818/150‑1 (MQ=255)
                          cTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                                           >  1:1280217/1‑150 (MQ=255)
                          cTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                                           <  1:212162/150‑1 (MQ=255)
                          cTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                                           >  1:1961086/1‑150 (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                                                                                                          <  1:1113804/150‑1 (MQ=255)
                                           ctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                          >  2:1729026/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)
                                                       gATACAGCATAAAGAAAGCTGGGCCATGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAATGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTa                                                                                              <  1:1901054/150‑1 (MQ=255)
                                                       gATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTa                                                                                              >  1:1097893/1‑150 (MQ=255)
                                                        aTACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTAt                                                                                             >  1:1732722/1‑150 (MQ=255)
                                                           cAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGa                                                                                          >  1:318893/1‑150 (MQ=255)
                                                              cATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAgg                                                                                       >  1:493081/1‑150 (MQ=255)
                                                               aTAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGAATCGAAggg                                                                                      >  2:402259/1‑150 (MQ=255)
                                                                 aaagaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCATTTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGc                                                                                    >  1:30950/1‑150 (MQ=255)
                                                                   agaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACAGTATCGAAGGGGCCa                                                                                  >  2:558458/1‑150 (MQ=255)
                                                                    taaagCTGGGCCAGCAAAAAATAACAGCAATCAAGGAACAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATAAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAg                                                                                 <  1:1744706/149‑1 (MQ=255)
                                                                    gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAg                                                                                 >  2:920214/1‑150 (MQ=255)
                                                                    gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAg                                                                                 >  1:1930372/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)
                                                                            ggCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTAGAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATg                                                                         >  1:139667/1‑150 (MQ=255)
                                                                             gCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgg                                                                        >  1:716184/1‑150 (MQ=255)
                                                                              ccAGGAAAAAATAACAGCAATCAAGGACAAGGGGAAACTCAGCAATTTTCTGAAAGATATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                                       <  1:586266/150‑1 (MQ=255)
                                                                               cAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                      <  2:1263071/150‑1 (MQ=255)
                                                                                 ggAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                    >  2:1619800/1‑150 (MQ=255)
                                                                                 ggAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                    >  1:1559630/1‑150 (MQ=255)
                                                                                    aaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGcc                                                                 <  2:57608/150‑1 (MQ=255)
                                                                                    aaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGcc                                                                 >  1:1413854/1‑150 (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)
                                                                                                aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGa                                                     >  2:1712774/1‑150 (MQ=255)
                                                                                                   cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAg                                                  >  1:1868178/1‑150 (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)
                                                                                                         aaaGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTTATCTTAAAACGCACCACTGTATCGAAGGGGCCCGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGc                                            >  2:1023192/1‑150 (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                                       >  2:1227637/1‑150 (MQ=255)
                                                                                                              ggAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtat                                       <  1:1793865/150‑1 (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)
                                                                                                                    tCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAATGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGGTCAGCCGGATAGCttt                                 >  1:814390/1‑150 (MQ=255)
                                                                                                                    tCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGttt                                 <  2:1130254/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)
                                                                                                                          aTTTTCTGAAAGAGAAGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGt                           <  2:142848/150‑1 (MQ=255)
                                                                                                                           ttttCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTa                          <  2:1732722/150‑1 (MQ=255)
                                                                                                                               cTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGAcc                      >  1:1664988/1‑150 (MQ=255)
                                                                                                                                tGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCa                     >  1:331423/1‑150 (MQ=255)
                                                                                                                                   aaGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACa                  <  2:477474/150‑1 (MQ=255)
                                                                                                                                   aaGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAACCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACa                  <  1:976677/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)
                                                                                                                                          gCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGacacg           >  1:1713002/1‑149 (MQ=255)
                                                                                                                                               ttCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGa      <  2:1572381/150‑1 (MQ=255)
                                                                                                                                                   ggcggcAAAGGGCATCATGTTAAATCAGTTATAGTTAATCTTAAAAAGCCCCACTGTATCGAAGGGGCCAGCCTAGATGGGGGTGCATGGTATGGACGAAAGTACAGCCGTATTGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATc  >  1:270299/1‑150 (MQ=255)
                                                                                                                                                   ggcggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCCGCCTCGAAGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATc  >  2:1269435/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATC  >  NC_000913/747750‑748046

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: