Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 747,900 G→C 9.0% intergenic (‑131/+21) 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,9000GC 9.0% 359.4 / 18.6 108intergenic (‑131/+21)abrB/ybgOregulator of aidB expression; inner membrane protein/putative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (50/48);  new base C (3/7);  total (53/55)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

CCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCG  >  NC_000913/747752‑748047
                                                                                                                                                     |                                                                                                                                                   
ccACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑Ccggcg                                                                                                                                                     >  4:367955/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)
        aaaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAggg                                                                                                                                             <  4:604855/149‑1 (MQ=255)
         aaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑GCCGCCGGAAAGGGc                                                                                                                                            <  4:166318/149‑1 (MQ=255)
         aaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑GCCGCCGGAAAGGGc                                                                                                                                            <  1:537349/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)
                aTCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATgg                                                                                                                                     >  4:629489/1‑149 (MQ=255)
                aTCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATgg                                                                                                                                     >  3:612802/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)
                               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)
                                              gtgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                       <  1:55449/149‑1 (MQ=255)
                                               tgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                                      <  1:487989/149‑1 (MQ=255)
                                               tgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                                      <  2:289755/149‑1 (MQ=255)
                                                 tATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcacca                                                                                                    >  2:860424/1‑149 (MQ=255)
                                                   tCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACt                                                                                                  >  1:319641/1‑149 (MQ=255)
                                                           gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                          <  1:138605/149‑1 (MQ=255)
                                                           gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                          <  3:342517/149‑1 (MQ=255)
                                                           gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                          <  2:42312/149‑1 (MQ=255)
                                                           gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                          <  2:416080/149‑1 (MQ=255)
                                                             aTAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                                 <  1:461198/124‑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                                                                                                <  2:434291/136‑1 (MQ=255)
                                                                  gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGt                                                                                                >  1:434291/1‑136 (MQ=255)
                                                                  gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCa                                                                                   >  1:706148/1‑149 (MQ=255)
                                                                  gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCa                                                                                   >  3:689727/1‑149 (MQ=255)
                                                                  gaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCa                                                                                   <  2:535563/149‑1 (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                                                                       <  4:3435/149‑1 (MQ=255)
                                                                              aGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                       <  2:559541/149‑1 (MQ=255)
                                                                              aGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                       <  2:581698/149‑1 (MQ=255)
                                                                                gtaaaaaTAACAGAAATAAAGGAAAGGGGGAAAATCAGCAATTTTCTGAAAAAGATCCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAGTCAGTTCTCGCTACTCTTAAAAAGCACCACTGTATCGGAAGTGCCAGCCTCGATGGgtgt                                                                     >  2:62701/3‑149 (MQ=255)
                                                                                gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                     >  2:227459/1‑149 (MQ=255)
                                                                                gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                     >  3:763888/1‑149 (MQ=255)
                                                                                gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                     >  1:376335/1‑149 (MQ=255)
                                                                                gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                     >  2:671343/1‑149 (MQ=255)
                                                                                gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                     >  1:715877/1‑149 (MQ=255)
                                                                                    aaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                                 <  3:749044/149‑1 (MQ=255)
                                                                                    aaaTAACAGCAAGCAGGGAAAAGTGGAAAAGCAGCAACTTTCTGAAAATGCTGCCCTTT‑GCGGCGGCAAAGGGCATCACGGTAAATGAGTTATAGTTAATCTGAAAAAGCACCACTGTATCGATAGTGCCAGCCTCGATTGTGGTGcat                                                                 >  4:655282/1‑147 (MQ=255)
                                                                                     aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                                <  1:496348/149‑1 (MQ=255)
                                                                                     aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                                >  4:415663/1‑149 (MQ=255)
                                                                                     aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                                >  2:849945/1‑149 (MQ=255)
                                                                                     aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                                >  3:868701/1‑149 (MQ=255)
                                                                                     aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                                >  3:532018/1‑149 (MQ=255)
                                                                                        aaCAGCAATCAAGGAAAAGGGGAAAATCAGTAATTTTCTGAAAGATATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAACCAGTTATAGTTAATCTTAAAAAGCACCGCTGTATCGAAGGGGCAAGCTTCGATGGGGGTGCCTGGaa                                                             <  2:629258/149‑1 (MQ=255)
                                                                                         aCAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAAt                                                            <  2:243578/149‑1 (MQ=255)
                                                                                           aGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                          >  4:506103/1‑149 (MQ=255)
                                                                                           aGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                          >  3:213383/1‑149 (MQ=255)
                                                                                            gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                       <  1:393300/103‑1 (MQ=255)
                                                                                            gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                       >  2:393300/1‑103 (MQ=255)
                                                                                            gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                         <  2:551985/149‑1 (MQ=255)
                                                                                            gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                         <  3:271012/149‑1 (MQ=255)
                                                                                             cAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGAc                                                        <  2:706148/149‑1 (MQ=255)
                                                                                             cAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGAc                                                        <  1:382992/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‑CCGGCGGCAAAGGGCATCATGGTAAATCTGTTATAGTTAATCTTAAAGAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGttt                                  <  1:802163/149‑1 (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)
                                                                                                                             cTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTATATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTGCGTAGAc                        <  2:717910/149‑1 (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                                                                                 >  4:4426/1‑79 (MQ=255)
                                                                                                                                           cTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTgg          >  1:41488/1‑149 (MQ=255)
                                                                                                                                                    cggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGTACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCg  <  1:262821/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCG  >  NC_000913/747752‑748047

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