Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,902,494 G→T 11.4% G200V (GGA→GTA) ‡ ygcE → putative sugar kinase YgcE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,902,4940GT11.4% 296.9 / 18.3 96G200V (GGA→GTA) ‡ygcEputative sugar kinase YgcE
Reads supporting (aligned to +/- strand):  ref base G (48/36);  new base T (4/7);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GGATTATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTAC  >  NC_000913/2902345‑2902642
                                                                                                                                                     |                                                                                                                                                    
ggATTATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAgg                                                                                                                                                      <  2:85143/150‑1 (MQ=255)
  aTTATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGAc                                                                                                                                                    >  1:1562696/1‑150 (MQ=255)
    tataACAGCTGGCAATGGAGCGAATATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCAACACATTTCCCGGCAGGACTg                                                                                                                                                  >  2:1427151/1‑150 (MQ=255)
        agagCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgg                                                                                                                                              >  1:256754/1‑150 (MQ=255)
         gagCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggg                                                                                                                                             <  2:1624102/150‑2 (MQ=255)
          agCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgtt                                                                                                                                            >  1:891356/1‑150 (MQ=255)
          agCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgtt                                                                                                                                            >  1:822636/1‑150 (MQ=255)
            cTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGTCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaa                                                                                                                                          <  2:788880/150‑5 (MQ=255)
            cTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaa                                                                                                                                          <  2:191483/150‑5 (MQ=255)
             tGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                         <  1:1100900/150‑1 (MQ=255)
             tGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                         <  1:852687/150‑1 (MQ=255)
              gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaatg                                                                                                                                        <  2:482720/150‑7 (MQ=255)
              gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaatg                                                                                                                                        <  2:256754/150‑7 (MQ=255)
              gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaatg                                                                                                                                        <  1:4371/150‑7 (MQ=255)
              gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaat                                                                                                                                         >  2:996891/1‑144 (MQ=255)
              gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaat                                                                                                                                         <  1:996891/149‑6 (MQ=255)
              gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGttgttt                                                                                                                                        >  2:105213/1‑150 (MQ=255)
               ggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGttgtttg                                                                                                                                       <  1:765419/150‑1 (MQ=255)
               ggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGttgtttg                                                                                                                                       <  1:1179432/150‑1 (MQ=255)
                    ggAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTacca                                                                                                                                  >  1:1560921/1‑150 (MQ=255)
                     gAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATCTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTaccac                                                                                                                                 >  2:527324/1‑150 (MQ=255)
                      aGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTaccacc                                                                                                                                <  1:1452572/150‑1 (MQ=255)
                        cGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAg                                                                                                                              >  2:384027/1‑150 (MQ=255)
                          aaGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGGTCTTTGATGTCCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAATCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTg                                                                                                                            >  1:301145/1‑150 (MQ=255)
                           aGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGa                                                                                                                           <  2:891356/150‑1 (MQ=255)
                           aGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTAATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCATTGa                                                                                                                           <  2:321577/150‑1 (MQ=255)
                            gATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTg                                                                                                                            >  1:155470/1‑148 (MQ=255)
                            gATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTg                                                                                                                            <  2:155470/148‑1 (MQ=255)
                                ctgcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAAc                                                                                                                      <  1:1473450/150‑1 (MQ=255)
                                 tgcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAAcc                                                                                                                     >  1:789232/1‑150 (MQ=255)
                                 tgcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAAcc                                                                                                                     <  2:1036431/150‑1 (MQ=255)
                                 tgcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTGGATGTGCAAATGCCTGGCACCGCCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAAcc                                                                                                                     <  2:14239/150‑1 (MQ=255)
                                    ggTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGt                                                                                                                  >  2:1009151/1‑150 (MQ=255)
                                             tAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGctct                                                                                                         >  1:1343691/1‑150 (MQ=255)
                                               aGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTgg                                                                                                       >  1:1323571/1‑150 (MQ=255)
                                                gTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggg                                                                                                      <  1:327614/150‑1 (MQ=255)
                                                gTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggg                                                                                                      >  1:1393772/1‑150 (MQ=255)
                                                gTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggg                                                                                                      >  2:1455935/1‑150 (MQ=255)
                                                    aaTATCCCCCGTCATATGCTGTTTGATCTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGcc                                                                                                  >  1:188441/1‑150 (MQ=255)
                                                     atatCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCg                                                                                                 >  1:465613/1‑150 (MQ=255)
                                                      tatCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCgg                                                                                                >  2:324989/1‑150 (MQ=255)
                                                      tatCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCgg                                                                                                <  2:1031312/150‑1 (MQ=255)
                                                          ccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTa                                                                                            >  2:1525656/1‑150 (MQ=255)
                                                           cccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCCGGGGGCCGGATTAc                                                                                           <  2:900419/150‑1 (MQ=255)
                                                            ccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACt                                                                                          <  2:1101627/150‑1 (MQ=255)
                                                              gTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTgg                                                                                        <  2:1393772/150‑1 (MQ=255)
                                                                  tatGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatga                                                                                    <  1:1455935/150‑1 (MQ=255)
                                                                    tGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGatgatg                                                                                  <  1:1038273/150‑1 (MQ=255)
                                                                     gCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGAtgatga                                                                                 >  1:103179/1‑150 (MQ=255)
                                                                       tGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGaaa                                                                               >  1:104829/1‑150 (MQ=255)
                                                                          ttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcgg                                                                            <  2:1562696/150‑1 (MQ=255)
                                                                              tgtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGt                                                                        >  1:1370847/1‑150 (MQ=255)
                                                                                tgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTaa                                                                      >  2:1311925/1‑150 (MQ=255)
                                                                                 gCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAt                                                                     >  1:1381721/1‑150 (MQ=255)
                                                                                 gCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAt                                                                     <  1:527324/150‑1 (MQ=255)
                                                                                  cAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAtt                                                                    >  2:633781/1‑150 (MQ=255)
                                                                                     aTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTtct                                                                 <  2:675618/150‑1 (MQ=255)
                                                                                       gCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTctct                                                               >  2:541498/1‑150 (MQ=255)
                                                                                        ccTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                              >  1:640942/1‑150 (MQ=255)
                                                                                        ccTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                              >  1:191264/1‑150 (MQ=255)
                                                                                         cTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTg                                                             >  1:66830/1‑150 (MQ=255)
                                                                                           ggCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGc                                                           <  2:994846/150‑1 (MQ=255)
                                                                                             cACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCAc                                                         >  1:740170/1‑150 (MQ=255)
                                                                                             cACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCAc                                                         <  1:105213/150‑1 (MQ=255)
                                                                                                     tCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCg                                                 >  2:1231428/1‑150 (MQ=255)
                                                                                                        gACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCAt                                              <  1:384290/150‑1 (MQ=255)
                                                                                                        gACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCAt                                              >  1:1267630/1‑150 (MQ=255)
                                                                                                         aCATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCAtt                                             <  1:324989/150‑1 (MQ=255)
                                                                                                          cATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtg                                            <  2:730958/150‑1 (MQ=255)
                                                                                                          cATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtg                                            <  2:1381721/150‑1 (MQ=255)
                                                                                                            tatCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgat                                          >  1:497712/1‑150 (MQ=255)
                                                                                                            tatCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgat                                          >  1:281855/1‑150 (MQ=255)
                                                                                                            aatCCCACAACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGACGGTTGTTTGTACCCCCAGTGATAAAAAGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAAGGCGGTAATTTCTCTTGGCAGTTACATCGCATtgat                                          >  2:200448/2‑150 (MQ=255)
                                                                                                              tCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatga                                        <  2:740170/150‑1 (MQ=255)
                                                                                                               cacaCCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgat                                       >  2:1238330/1‑150 (MQ=255)
                                                                                                                   ccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAAc                                   >  2:117253/1‑150 (MQ=255)
                                                                                                                    gacaAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACg                                  >  1:1607197/2‑150 (MQ=255)
                                                                                                                    cacaAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACg                                  >  2:882488/1‑150 (MQ=255)
                                                                                                                      caAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGc                                <  1:1009151/150‑1 (MQ=255)
                                                                                                                        aGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCaa                              <  1:1525656/150‑1 (MQ=255)
                                                                                                                           cGCACTGCCGACAAATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTCATAAACCGGTAGAAGCTCTGGGGGCCGGATTACAGGATGATGAAACGGCGGTAATTTATCTTGGAACTTACATCGCATTGATGATGAACGGCAAAGc                           >  2:885038/1‑150 (MQ=255)
                                                                                                                            gCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCa                          >  1:1628788/1‑150 (MQ=255)
                                                                                                                              aCTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACt                        <  2:1142796/150‑1 (MQ=255)
                                                                                                                                 ggCGACACATTTCCCGGCAGCACTGACGGTTGTTTGTACCACAAGTGATAAACCGCTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATAGCATTGATGATGAACGGCAAAGCACTGcc                     >  2:579716/1‑150 (MQ=255)
                                                                                                                                    gACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                  <  2:640942/150‑1 (MQ=255)
                                                                                                                                    gACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                  >  1:1132300/1‑150 (MQ=255)
                                                                                                                                      cacaTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAg                >  2:1551269/1‑150 (MQ=255)
                                                                                                                                      cacaTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAg                >  2:1580321/1‑150 (MQ=255)
                                                                                                                                       caaacaaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGa               >  1:152660/8‑150 (MQ=255)
                                                                                                                                        caTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGAt              >  1:1552359/1‑150 (MQ=255)
                                                                                                                                         aTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATc             <  1:882488/150‑1 (MQ=255)
                                                                                                                                          tttCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATcc            <  2:1267630/150‑1 (MQ=255)
                                                                                                                                           ttCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCa           <  1:332757/150‑1 (MQ=255)
                                                                                                                                            tCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAg          >  1:10470/1‑150 (MQ=255)
                                                                                                                                            aaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAg          >  1:1069163/3‑150 (MQ=255)
                                                                                                                                              ccGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTa        >  2:1451356/1‑150 (MQ=255)
                                                                                                                                               cGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAg       <  2:465613/150‑1 (MQ=255)
                                                                                                                                                    ggACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTAc  >  1:830345/1‑150 (MQ=255)
                                                                                                                                                    ggACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTAc  <  2:443156/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGATTATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTAC  >  NC_000913/2902345‑2902642

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