Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 747,902 C→G 7.5% 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.5% 359.1 / 15.0 108intergenic (‑133/+19)abrB/ybgOregulator of aidB expression; inner membrane protein/putative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base C (51/48);  new base G (2/6);  total (53/55)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 110 total aligned reads displayed.

CAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGC  >  NC_000913/747758‑748049
                                                                                                                                                 |                                                                                                                                                   
cAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAg                                                                                                                                                 <  3:379647/149‑1 (MQ=255)
 aaaaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                                                <  1:114956/149‑1 (MQ=255)
  aaaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAggg                                                                                                                                               <  4:604855/149‑1 (MQ=255)
   aaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGc                                                                                                                                              <  4:166318/149‑1 (MQ=255)
   aaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGc                                                                                                                                              <  1:537349/149‑1 (MQ=255)
       ggCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGcatct                                                                                                                                          <  3:89230/149‑2 (MQ=255)
       ggCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcatca                                                                                                                                          <  3:381058/149‑1 (MQ=255)
       ggCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcatca                                                                                                                                          <  1:18424/149‑1 (MQ=255)
        gCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcatcat                                                                                                                                         >  4:271012/1‑149 (MQ=255)
          aTCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATgg                                                                                                                                       >  4:629489/1‑149 (MQ=255)
          aTCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATgg                                                                                                                                       >  3:612802/1‑149 (MQ=255)
           tCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGt                                                                                                                                      >  3:870543/1‑149 (MQ=255)
            ccTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGcatctctttc                                                                                                                                     <  4:449152/149‑7 (MQ=255)
              ttCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTaaa                                                                                                                                   >  2:677135/1‑149 (MQ=255)
                   tCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGt                                                                                                                                      <  3:148672/141‑1 (MQ=255)
                   tCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGt                                                                                                                                      >  4:148672/1‑141 (MQ=255)
                      cTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTtat                                                                                                                           <  1:677135/149‑1 (MQ=255)
                      cTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTtat                                                                                                                           >  4:70212/1‑149 (MQ=255)
                       tCATTATTTTCTCTTTGGTTTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                          <  2:641116/149‑1 (MQ=255)
                       tCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                          >  4:598572/1‑149 (MQ=255)
                         attattTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGt                                                                                                                        >  3:345369/1‑149 (MQ=255)
                            attTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGAGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGttta                                                                                                                     >  2:474590/1‑147 (MQ=255)
                             ttttCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAAt                                                                                                                    <  4:119289/149‑1 (MQ=255)
                                tctcTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCtt                                                                                                                 >  1:527051/1‑149 (MQ=255)
                                    tttGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                             >  2:311763/1‑149 (MQ=255)
                                        gtgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                         <  1:55449/149‑1 (MQ=255)
                                         tgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAGGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                                        <  2:92003/149‑1 (MQ=255)
                                         tgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                                        <  1:487989/149‑1 (MQ=255)
                                         tgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                                        <  2:289755/149‑1 (MQ=255)
                                           tATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcacca                                                                                                      >  2:860424/1‑149 (MQ=255)
                                             tCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACt                                                                                                    >  1:319641/1‑149 (MQ=255)
                                                     gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                            <  2:42312/149‑1 (MQ=255)
                                                     gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                            <  2:416080/149‑1 (MQ=255)
                                                     gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                            <  1:138605/149‑1 (MQ=255)
                                                       aTAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                                   >  2:461198/1‑124 (MQ=255)
                                                       aTAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                                   <  1:461198/124‑1 (MQ=255)
                                                        tAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                                   <  3:775659/123‑1 (MQ=255)
                                                        tAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                                   >  4:775659/1‑123 (MQ=255)
                                                          aagaaagCTGGGCCAGGAAAAAATTACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGc                                                                                       >  4:640821/1‑149 (MQ=255)
                                                          aagaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATAATGGTAAATAAGATATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGc                                                                                       >  2:744866/1‑149 (MQ=255)
                                                            gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGt                                                                                                  <  2:434291/136‑1 (MQ=255)
                                                            gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGt                                                                                                  >  1:434291/1‑136 (MQ=255)
                                                            gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCa                                                                                     <  2:535563/149‑1 (MQ=255)
                                                            gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCa                                                                                     >  1:706148/1‑149 (MQ=255)
                                                             aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAg                                                                                    >  1:771175/1‑149 (MQ=255)
                                                              aagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGc                                                                                   <  1:53814/149‑1 (MQ=255)
                                                                gCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATg                                                                                                                                        <  3:723617/94‑1 (MQ=255)
                                                                gCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATg                                                                                                                                        >  4:723617/1‑94 (MQ=255)
                                                                  tGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCg                                                                               >  2:591271/1‑149 (MQ=255)
                                                                        aGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                         <  2:559541/149‑1 (MQ=255)
                                                                        aGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                         <  4:3435/149‑1 (MQ=255)
                                                                        aGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                         <  2:581698/149‑1 (MQ=255)
                                                                          gtaaaaaTAACAGAAATAAAGGAAAGGGGGAAAATCAGCAATTTTCTGAAAAAGATCCCCTTTCCGGCGGCAAAGGGCATCATGGTAAGTCAGTTCTCGCTACTCTTAAAAAGCACCACTGTATCGGAAGTGCCAGCCTCGATGGgtgt                                                                       >  2:62701/3‑149 (MQ=255)
                                                                          gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                       >  2:671343/1‑149 (MQ=255)
                                                                          gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                       >  1:715877/1‑149 (MQ=255)
                                                                          gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                       >  3:763888/1‑149 (MQ=255)
                                                                          gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                       >  2:227459/1‑149 (MQ=255)
                                                                          gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                       >  1:376335/1‑149 (MQ=255)
                                                                              aaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                                   <  3:749044/149‑1 (MQ=255)
                                                                              aaaTAACAGCAAGCAGGGAAAAGTGGAAAAGCAGCAACTTTCTGAAAATGCTGCCCTTTGCGGCGGCAAAGGGCATCACGGTAAATGAGTTATAGTTAATCTGAAAAAGCACCACTGTATCGATAGTGCCAGCCTCGATTGTGGTGcat                                                                   >  4:655282/1‑147 (MQ=255)
                                                                               aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                                  >  4:415663/1‑149 (MQ=255)
                                                                               aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                                  >  2:849945/1‑149 (MQ=255)
                                                                               aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                                  >  3:532018/1‑149 (MQ=255)
                                                                               aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                                  <  1:496348/149‑1 (MQ=255)
                                                                               aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                                  >  3:868701/1‑149 (MQ=255)
                                                                                  aaCAGCAATCAAGGAAAAGGGGAAAATCAGTAATTTTCTGAAAGATATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAACCAGTTATAGTTAATCTTAAAAAGCACCGCTGTATCGAAGGGGCAAGCTTCGATGGGGGTGCCTGGaa                                                               <  2:629258/149‑1 (MQ=255)
                                                                                  aaCAGCAATCAAGGAAAAGGAGAAAATCAGCAATTACCTGAAAGAGATGCCCTTTCCGGCCGCCAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTCCCTGGaa                                                               <  2:105197/149‑1 (MQ=255)
                                                                                     aGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                            >  4:506103/1‑149 (MQ=255)
                                                                                     aGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                            >  3:213383/1‑149 (MQ=255)
                                                                                      gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                         <  1:393300/103‑1 (MQ=255)
                                                                                      gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                         >  2:393300/1‑103 (MQ=255)
                                                                                      gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                           <  2:551985/149‑1 (MQ=255)
                                                                                      gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                           <  3:271012/149‑1 (MQ=255)
                                                                                      gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAATTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                           <  2:63616/149‑1 (MQ=255)
                                                                                       cAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGAc                                                          <  1:382992/149‑1 (MQ=255)
                                                                                           cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                                      <  1:576982/149‑1 (MQ=255)
                                                                                           cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                                      <  3:254915/149‑1 (MQ=255)
                                                                                           cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCCTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                                      <  3:79667/149‑1 (MQ=255)
                                                                                              ggAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGtt                                                   >  4:553850/1‑149 (MQ=255)
                                                                                                             cAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCTGTTATAGTTAATCTTAAAGAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGttt                                    <  1:802163/149‑1 (MQ=255)
                                                                                                             cAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGttt                                    <  2:895637/149‑1 (MQ=255)
                                                                                                               gCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGttttt                                  >  1:811307/1‑149 (MQ=255)
                                                                                                                    tttCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAGTCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTa                             <  3:566950/149‑1 (MQ=255)
                                                                                                                    tttCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTa                             <  3:640821/149‑1 (MQ=255)
                                                                                                                     ttCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAg                            >  1:484976/1‑149 (MQ=255)
                                                                                                                       cTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTATATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTGCGTAGAc                          <  2:717910/149‑1 (MQ=255)
                                                                                                                          aaaGAGATGCCCTTTCCGGCCGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCaa                       >  1:714578/1‑149 (MQ=255)
                                                                                                                           aaGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACGAAc                      >  3:501661/1‑149 (MQ=255)
                                                                                                                           aaGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAAc                      >  3:439766/1‑149 (MQ=255)
                                                                                                                                 tGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCa                <  4:400508/149‑1 (MQ=255)
                                                                                                                                  gCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtata                                         >  2:272294/1‑123 (MQ=255)
                                                                                                                                  gCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtata                                         <  1:272294/123‑1 (MQ=255)
                                                                                                                                  gCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCAc               >  4:696918/1‑149 (MQ=255)
                                                                                                                                   cccTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                                      >  3:850287/1‑109 (MQ=255)
                                                                                                                                   cccTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                                      <  4:850287/109‑1 (MQ=255)
                                                                                                                                    ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGc                                                                                   >  4:4426/1‑79 (MQ=255)
                                                                                                                                    ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGc                                                                                   <  3:4426/79‑1 (MQ=255)
                                                                                                                                     cTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTgg            >  1:41488/1‑149 (MQ=255)
                                                                                                                                              cggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGTACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCg    <  1:262821/149‑1 (MQ=255)
                                                                                                                                                gcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGc  <  2:200203/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGC  >  NC_000913/747758‑748049

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