Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 747,895 C→G 8.0% intergenic (‑126/+26) 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,8950CG 8.0% 365.9 / 15.5 111intergenic (‑126/+26)abrB/ybgOregulator of aidB expression; inner membrane protein/putative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base C (52/50);  new base G (3/6);  total (55/56)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 112 total aligned reads displayed.

ATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGG  >  NC_000913/747747‑748039
                                                                                                                                                     |                                                                                                                                                
aTTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑c                                                                                                                                                  <  4:139651/149‑1 (MQ=255)
   ccccACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑Ccgg                                                                                                                                               >  1:816883/1‑149 (MQ=255)
     ccACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑Ccggcg                                                                                                                                             <  1:656015/149‑1 (MQ=255)
           cAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑GCCGCCGGAAAg                                                                                                                                       <  4:133725/149‑1 (MQ=255)
           cAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAg                                                                                                                                       <  3:379647/149‑1 (MQ=255)
            aaaaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAgg                                                                                                                                      <  1:114956/149‑1 (MQ=255)
              aaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑GCCGCCGGAAAGGGc                                                                                                                                    <  1:537349/149‑1 (MQ=255)
              aaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑GCCGCCGGAAAGGGc                                                                                                                                    <  4:166318/149‑1 (MQ=255)
                  ggCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGG‑AAAGGGcatct                                                                                                                                <  3:89230/149‑2 (MQ=255)
                  ggCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGcatca                                                                                                                                <  1:18424/149‑1 (MQ=255)
                  ggCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGcatca                                                                                                                                <  3:381058/149‑1 (MQ=255)
                   gCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGcatcat                                                                                                                               >  4:271012/1‑149 (MQ=255)
                     aTCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATgg                                                                                                                             >  4:629489/1‑149 (MQ=255)
                      tCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGt                                                                                                                            >  3:870543/1‑149 (MQ=255)
                       ccTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑GCCGCCGGAAAGGGcatctctttc                                                                                                                           <  4:449152/149‑7 (MQ=255)
                        cTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGGCCTTT‑CCGGCGGCAAAGGGCATCATGGTaa                                                                                                                          >  1:581698/1‑149 (MQ=255)
                         ttCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTaaa                                                                                                                         >  2:677135/1‑149 (MQ=255)
                              tCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGt                                                                                                                            >  4:148672/1‑141 (MQ=255)
                              tCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGt                                                                                                                            <  3:148672/141‑1 (MQ=255)
                                 cTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtat                                                                                                                 >  4:70212/1‑149 (MQ=255)
                                 cTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtat                                                                                                                 <  1:677135/149‑1 (MQ=255)
                                  tCATTATTTTCTCTTTGGTTTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                <  2:641116/149‑1 (MQ=255)
                                  tCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                >  2:538529/1‑149 (MQ=255)
                                  tCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                >  4:598572/1‑149 (MQ=255)
                                    attattTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGt                                                                                                              >  3:345369/1‑149 (MQ=255)
                                       attTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGAGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGttta                                                                                                           >  2:474590/1‑147 (MQ=255)
                                        ttttCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAAt                                                                                                          <  4:119289/149‑1 (MQ=255)
                                           tctcTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCtt                                                                                                       >  1:527051/1‑149 (MQ=255)
                                               tttGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                   >  2:311763/1‑149 (MQ=255)
                                                    tgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAGGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                              <  2:92003/149‑1 (MQ=255)
                                                    tgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                              <  2:289755/149‑1 (MQ=255)
                                                    tgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                              <  1:487989/149‑1 (MQ=255)
                                                      tATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcacca                                                                                            >  2:860424/1‑149 (MQ=255)
                                                        tCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACt                                                                                          >  1:319641/1‑149 (MQ=255)
                                                                gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                  <  2:42312/149‑1 (MQ=255)
                                                                gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                  <  1:138605/149‑1 (MQ=255)
                                                                  aTAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                         <  1:461198/124‑1 (MQ=255)
                                                                  aTAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                         >  2:461198/1‑124 (MQ=255)
                                                                   tAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                         <  3:775659/123‑1 (MQ=255)
                                                                   tAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                         >  4:775659/1‑123 (MQ=255)
                                                                     aagaaagCTGGGCCAGGAAAAAATTACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGc                                                                             >  4:640821/1‑149 (MQ=255)
                                                                     aagaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATAATGGTAAATAAGATATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGc                                                                             >  2:744866/1‑149 (MQ=255)
                                                                       gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGt                                                                                        >  1:434291/1‑136 (MQ=255)
                                                                       gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGt                                                                                        <  2:434291/136‑1 (MQ=255)
                                                                       gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCa                                                                           >  3:689727/1‑149 (MQ=255)
                                                                       gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCa                                                                           <  2:535563/149‑1 (MQ=255)
                                                                       gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCa                                                                           >  1:706148/1‑149 (MQ=255)
                                                                        aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAg                                                                          >  1:771175/1‑149 (MQ=255)
                                                                         aagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGc                                                                         <  1:53814/149‑1 (MQ=255)
                                                                           gCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATg                                                                                                                              >  4:723617/1‑94 (MQ=255)
                                                                           gCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATg                                                                                                                              <  3:723617/94‑1 (MQ=255)
                                                                           gCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCt                                                                       >  4:566950/1‑149 (MQ=255)
                                                                             tGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCg                                                                     >  2:591271/1‑149 (MQ=255)
                                                                                   aGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                               <  2:559541/149‑1 (MQ=255)
                                                                                   aGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                               <  4:3435/149‑1 (MQ=255)
                                                                                     gtaaaaaTAACAGAAATAAAGGAAAGGGGGAAAATCAGCAATTTTCTGAAAAAGATCCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAGTCAGTTCTCGCTACTCTTAAAAAGCACCACTGTATCGGAAGTGCCAGCCTCGATGGgtgt                                                             >  2:62701/3‑149 (MQ=255)
                                                                                     gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                             >  1:715877/1‑149 (MQ=255)
                                                                                     gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                             >  1:376335/1‑149 (MQ=255)
                                                                                     gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                             >  3:763888/1‑149 (MQ=255)
                                                                                     gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                             >  2:227459/1‑149 (MQ=255)
                                                                                     gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                             >  2:671343/1‑149 (MQ=255)
                                                                                         aaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                         <  3:749044/149‑1 (MQ=255)
                                                                                         aaaTAACAGCAAGCAGGGAAAAGTGGAAAAGCAGCAACTTTCTGAAAATGCTGCCCTTT‑GCGGCGGCAAAGGGCATCACGGTAAATGAGTTATAGTTAATCTGAAAAAGCACCACTGTATCGATAGTGCCAGCCTCGATTGTGGTGcat                                                         >  4:655282/1‑147 (MQ=255)
                                                                                          aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                        >  4:415663/1‑149 (MQ=255)
                                                                                          aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                        >  3:532018/1‑149 (MQ=255)
                                                                                          aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                        <  1:496348/149‑1 (MQ=255)
                                                                                          aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                        >  2:849945/1‑149 (MQ=255)
                                                                                             aaCAGCAATCAAGGAAAAGGGGAAAATCAGTAATTTTCTGAAAGATATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAACCAGTTATAGTTAATCTTAAAAAGCACCGCTGTATCGAAGGGGCAAGCTTCGATGGGGGTGCCTGGaa                                                     <  2:629258/149‑1 (MQ=255)
                                                                                             aaCAGCAATCAAGGAAAAGGAGAAAATCAGCAATTACCTGAAAGAGATGCCCTTT‑CCGGCCGCCAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTCCCTGGaa                                                     <  2:105197/149‑1 (MQ=255)
                                                                                              aCAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAAt                                                    <  2:243578/149‑1 (MQ=255)
                                                                                                aGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                  >  3:213383/1‑149 (MQ=255)
                                                                                                aGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                  >  4:506103/1‑149 (MQ=255)
                                                                                                 gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                               <  1:393300/103‑1 (MQ=255)
                                                                                                 gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                               >  2:393300/1‑103 (MQ=255)
                                                                                                 gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                 <  3:271012/149‑1 (MQ=255)
                                                                                                 gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                 <  2:551985/149‑1 (MQ=255)
                                                                                                 gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAATTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                 <  2:63616/149‑1 (MQ=255)
                                                                                                  cAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGAc                                                <  2:706148/149‑1 (MQ=255)
                                                                                                  cAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGAc                                                <  1:382992/149‑1 (MQ=255)
                                                                                                      cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                            <  1:576982/149‑1 (MQ=255)
                                                                                                      cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                            <  3:254915/149‑1 (MQ=255)
                                                                                                      cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCCTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                            <  3:79667/149‑1 (MQ=255)
                                                                                                         ggAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGtt                                         >  4:553850/1‑149 (MQ=255)
                                                                                                                        cAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGttt                          <  2:895637/149‑1 (MQ=255)
                                                                                                                          gCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGttttt                        >  1:811307/1‑149 (MQ=255)
                                                                                                                               tttCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAGTCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTa                   <  3:566950/149‑1 (MQ=255)
                                                                                                                               tttCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTa                   <  3:640821/149‑1 (MQ=255)
                                                                                                                                ttCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAg                  >  1:484976/1‑149 (MQ=255)
                                                                                                                                     aaaGAGATGCCCTTT‑CCGGCCGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCaa             >  1:714578/1‑149 (MQ=255)
                                                                                                                                      aaGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACGAAc            >  3:501661/1‑149 (MQ=255)
                                                                                                                                      aaGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAAc            >  3:439766/1‑149 (MQ=255)
                                                                                                                                            tGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCa      <  4:400508/149‑1 (MQ=255)
                                                                                                                                             gCCCTTT‑GCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtata                               >  2:272294/1‑123 (MQ=255)
                                                                                                                                             gCCCTTT‑GCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtata                               <  1:272294/123‑1 (MQ=255)
                                                                                                                                             gCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCAc     >  4:696918/1‑149 (MQ=255)
                                                                                                                                              cccTTT‑GCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                            >  3:850287/1‑109 (MQ=255)
                                                                                                                                              cccTTT‑GCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                            <  4:850287/109‑1 (MQ=255)
                                                                                                                                               ccTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGc                                                                         <  3:4426/79‑1 (MQ=255)
                                                                                                                                               ccTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGc                                                                         >  4:4426/1‑79 (MQ=255)
                                                                                                                                                cTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTgg  >  1:41488/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
ATTCCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGG  >  NC_000913/747747‑748039

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: ‑