Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,795,229 A→T 6.1% intergenic (‑77/+24) btuC ← / ← ihfA vitamin B12 ABC transporter permease/integration host factor (IHF), DNA‑binding protein, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,795,2290AT 6.1% 366.3 / 9.9 109intergenic (‑77/+24)btuC/ihfAvitamin B12 ABC transporter permease/integration host factor (IHF), DNA‑binding protein, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (54/48);  new base T (4/3);  total (58/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 110 total aligned reads displayed.

GCATCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTT  >  NC_000913/1795083‑1795377
                                                                                                                                                  |                                                                                                                                                    
gcaTCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGc                                                                                                                                                    >  1:376607/1‑149 (MQ=255)
 caTCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCg                                                                                                                                                   <  1:63681/149‑1 (MQ=255)
   tCAAACCTGACAGGCATAATCACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCACGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGCAAAAAGGCCGCAGAGCGGc                                                                                                                                                 >  4:688392/1‑149 (MQ=255)
   tCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGc                                                                                                                                                 >  4:869716/1‑149 (MQ=255)
   tCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTAATGGCGTTGATTTTACGGAGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGc                                                                                                                                                 <  2:570516/149‑1 (MQ=255)
        aCTGACAGGCATAATAACCCGCGAATATTTTGTCGCTGGTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCtttt                                                                                                                                            <  1:747306/149‑1 (MQ=255)
         cTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCttttt                                                                                                                                           <  4:178833/149‑1 (MQ=255)
         cTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCttttt                                                                                                                                           <  4:143696/149‑1 (MQ=255)
          tGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTtta                                                                                                                                          <  3:491443/149‑1 (MQ=255)
          tGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCAGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTtta                                                                                                                                          <  4:373154/149‑1 (MQ=255)
          tGACAGGCATAATAACCAGCAAATATTTTGTCGCGGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGCTTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTtta                                                                                                                                          <  2:586569/149‑1 (MQ=255)
           gACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttag                                                                                                                                         <  3:293500/149‑1 (MQ=255)
            aCAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagt                                                                                                                                        <  4:671889/149‑1 (MQ=255)
             cagtcgtaataaCCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCAAGGTATCCGTACTGCTGAAGTGTCAAGGCGTTGATTTTACGGTGACTCTTCGACAGTGATAAGAAAAAAGGCCGCAGTGCGGGCTTTttagtt                                                                                                                                       <  4:679303/143‑1 (MQ=255)
             cAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagtt                                                                                                                                       <  2:512550/149‑1 (MQ=255)
              aGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAACGAAAAAAGGCCGCAGAGCGGCCTTTttagtta                                                                                                                                      <  4:129039/149‑1 (MQ=255)
              aGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttttttcttt                                                                                                                                      <  3:497377/149‑6 (MQ=255)
              aGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagtta                                                                                                                                      >  4:287745/1‑149 (MQ=255)
                gCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGa                                                                                                                                    >  3:843547/1‑149 (MQ=255)
                gCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGa                                                                                                                                    >  3:60298/1‑149 (MQ=255)
                 cATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGAt                                                                                                                                   <  3:985750/149‑1 (MQ=255)
                     ataaCCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGa                                                                                                                               <  3:287745/149‑1 (MQ=255)
                      taaCCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGAt                                                                                                                              >  1:359845/1‑149 (MQ=255)
                         ccAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTAc                                                                                                                           <  2:544987/149‑1 (MQ=255)
                         ccAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTAc                                                                                                                           <  3:869716/149‑1 (MQ=255)
                          cAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagt                                                                                                                                        >  3:509390/1‑135 (MQ=255)
                          cAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagt                                                                                                                                        <  4:509390/135‑1 (MQ=255)
                          cAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACt                                                                                                                          >  2:9959/1‑149 (MQ=255)
                          cAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACt                                                                                                                          >  3:22823/1‑149 (MQ=255)
                           aGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTc                                                                                                                         <  1:9959/149‑1 (MQ=255)
                           aGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTc                                                                                                                         <  1:817007/149‑1 (MQ=255)
                            gCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCg                                                                                                                        >  4:375335/1‑149 (MQ=255)
                            gCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCg                                                                                                                        >  2:49739/1‑149 (MQ=255)
                            gCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCg                                                                                                                        >  2:821225/1‑149 (MQ=255)
                              gAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTtta                                                                                                                                          <  1:543603/129‑1 (MQ=255)
                              gAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTtta                                                                                                                                          >  2:543603/1‑129 (MQ=255)
                                 tatTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCttt                                                                                                                   <  2:9128/149‑1 (MQ=255)
                                  atTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTg                                                                                                                  >  1:890806/1‑149 (MQ=255)
                                       gTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGa                                                                                                             >  1:602111/1‑149 (MQ=255)
                                        tCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGaa                                                                                                            <  2:739554/149‑1 (MQ=255)
                                         cGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAg                                                                                                           <  3:423289/149‑1 (MQ=255)
                                             gttgttGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGtt                                                                                                       <  4:337766/149‑1 (MQ=255)
                                              ttgttgGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGttt                                                                                                      <  4:737503/149‑1 (MQ=255)
                                                gttgGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTc                                                                                                    >  1:511165/1‑149 (MQ=255)
                                                   ggCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGAc                                                                                                 >  4:427973/1‑149 (MQ=255)
                                                      gggCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCg                                                                                              >  4:786810/1‑149 (MQ=255)
                                                       ggCAAGTGTCAGCATGGTTTCCGTTCTGCAGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAATTAGATCAGATTACTCGTCTTTGAGAGAAGCGTTTTCGGCCCgg                                                                                             >  2:826097/1‑149 (MQ=255)
                                                                 aGCATGGTATCCGTTCTGCTGAAGAGTCATGGCGTTGATTTTACGGTGACTCTTCTACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACt                                                                                   >  2:566250/1‑149 (MQ=255)
                                                                  gCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACtt                                                                                  <  3:625845/149‑1 (MQ=255)
                                                                        tATCCGTTCTGCTGAAGTGTCATGGCCTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGccc                                                                            <  3:541325/149‑1 (MQ=255)
                                                                         aTCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCg                                                                           >  2:941711/1‑149 (MQ=255)
                                                                         aTCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCg                                                                           >  1:685951/1‑149 (MQ=255)
                                                                          tCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGt                                                                                                                       >  1:312610/1‑104 (MQ=255)
                                                                          tCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCgg                                                                          >  2:840980/1‑149 (MQ=255)
                                                                             gTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTc                                                                       >  3:751947/1‑149 (MQ=255)
                                                                                 tgctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGAt                                                                                                                              <  3:436497/90‑1 (MQ=255)
                                                                                  gctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTc                                                                                                                      >  2:402828/1‑97 (MQ=255)
                                                                                  gctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTc                                                                                                                      <  1:402828/97‑1 (MQ=255)
                                                                                  gctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAg                                                                  <  1:1038882/149‑1 (MQ=255)
                                                                                    tgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGaa                                                                   >  4:40381/1‑146 (MQ=255)
                                                                                    tgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGt                                                                <  3:257155/149‑1 (MQ=255)
                                                                                     gAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGaa                                                                   <  3:297763/145‑1 (MQ=255)
                                                                                     gAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGaa                                                                   >  4:297763/1‑145 (MQ=255)
                                                                                      aaGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTca                                                              <  4:60298/149‑1 (MQ=255)
                                                                                      aaGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTca                                                              <  2:511165/149‑1 (MQ=255)
                                                                                                 agTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCTGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGt                                                   <  2:40909/148‑1 (MQ=255)
                                                                                                  gTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTg                                                  >  4:1031329/1‑149 (MQ=255)
                                                                                                   ttGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGc                                                 >  3:480459/1‑149 (MQ=255)
                                                                                                   ttGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGc                                                 >  1:355044/1‑149 (MQ=255)
                                                                                                           aCGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATgg                                         <  2:291693/149‑1 (MQ=255)
                                                                                                           aCGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATgg                                         <  3:499329/149‑1 (MQ=255)
                                                                                                             cgTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGa                                       <  1:941711/148‑1 (MQ=255)
                                                                                                                   tctTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATcc                                 >  2:693092/1‑149 (MQ=255)
                                                                                                                   tctTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATcc                                 >  1:546543/1‑149 (MQ=255)
                                                                                                                   tctTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATcc                                 >  1:30475/1‑149 (MQ=255)
                                                                                                                   tctTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAACGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTaa                                            <  3:248961/138‑1 (MQ=255)
                                                                                                                   tctTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAACGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTaa                                            >  4:248961/1‑138 (MQ=255)
                                                                                                                        gACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGcc                            >  2:978106/1‑149 (MQ=255)
                                                                                                                         aCAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGccc                           >  3:676306/1‑149 (MQ=255)
                                                                                                                            gTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGtt                        >  1:603194/1‑149 (MQ=255)
                                                                                                                            gTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGtt                        >  1:89261/1‑149 (MQ=255)
                                                                                                                              tctaactaaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGtttt                      >  4:196626/8‑149 (MQ=255)
                                                                                                                               aaaagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGc                                                 >  1:961636/1‑121 (MQ=255)
                                                                                                                               aaaagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTc                     >  1:1019168/1‑149 (MQ=255)
                                                                                                                                aaagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCg                    <  4:284885/149‑1 (MQ=255)
                                                                                                                                aaagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCg                    >  2:814378/1‑149 (MQ=255)
                                                                                                                                 aagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCgg                   <  3:169304/149‑1 (MQ=255)
                                                                                                                                   gtaaaaaGGCCGCAGAGCAGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCTAAGCGTTTTCGTACCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTCATGGGAATATCTTCGCCCGGTGTCGGGt                 >  2:742913/3‑149 (MQ=255)
                                                                                                                                   gaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGt                 <  4:751947/149‑1 (MQ=255)
                                                                                                                                   gaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGATTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGt                 >  3:374535/1‑149 (MQ=255)
                                                                                                                                    taaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGtt                >  1:852889/2‑149 (MQ=255)
                                                                                                                                    taaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGtt                >  2:312536/2‑149 (MQ=255)
                                                                                                                                     aaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTa               >  2:861754/1‑149 (MQ=255)
                                                                                                                                     aaaaaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTa               >  1:192765/1‑149 (MQ=255)
                                                                                                                                     aaaaaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTa               >  1:499294/1‑149 (MQ=255)
                                                                                                                                        aaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGt            >  1:486622/1‑149 (MQ=255)
                                                                                                                                           gCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTccc         >  4:590380/1‑149 (MQ=255)
                                                                                                                                           gCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTccc         >  2:192726/1‑149 (MQ=255)
                                                                                                                                             cGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCgg       <  1:287889/149‑1 (MQ=255)
                                                                                                                                                  agCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGtt  <  4:26449/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCATCAAAACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTT  >  NC_000913/1795083‑1795377

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