Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,902,495 A→C 5.2% G200G (GGA→GGC)  ygcE → putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,902,4950AC 5.2% 257.8 / 3.0 90G200G (GGA→GGC) ygcEputative kinase
Reads supporting (aligned to +/- strand):  ref base A (40/45);  new base C (2/3);  total (42/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

ATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACT  >  NC_000913/2902350‑2902643
                                                                                                                                                 |                                                                                                                                                    
ataAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGc                                                                                                                                                  >  2:57450/1‑150 (MQ=255)
 taAGAGCTGGGCATGGAGTGAAGTTGCTGCGGTAATGGATATGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACAGTCCTCGGACATATGACACCACTAGCCGCACTGGCGACACATTTCCCGGCAGGACTGcc                                                                                                                                                 <  2:172979/150‑1 (MQ=255)
 taAGAGCTGGGCATGGAGCGAAGAGGCGGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGcc                                                                                                                                                 <  2:205733/150‑1 (MQ=255)
   agagCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCgg                                                                                                                                               <  2:466674/150‑1 (MQ=255)
   agagCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAACTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTACCCGGCAGGACTGCCgg                                                                                                                                               >  2:242964/1‑150 (MQ=255)
    gagCTGGGCAAGGAGCGAAGATGCTGCGGAAATGGATAAGTGTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggg                                                                                                                                              <  2:102202/150‑2 (MQ=255)
     agCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATTAGTTTAATATCCCCCCTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACTCATTTCCCGGCAGGACTGCCGgtt                                                                                                                                             <  1:160453/150‑1 (MQ=255)
        tGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaat                                                                                                                                          <  2:1109669/150‑6 (MQ=255)
         gggCATGGAGCGAAGATGCTGCGGTAATGTATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaatg                                                                                                                                         <  2:523118/150‑7 (MQ=255)
         gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACAGTCCTCGGACATATCACACAACAAGCCGCACTGGCGCCACATTTCCCGGCAGGCCTGCCGGttgttt                                                                                                                                         >  1:713096/1‑150 (MQ=255)
          ggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGttgtttg                                                                                                                                        <  2:975394/150‑1 (MQ=255)
            cATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTa                                                                                                                                      <  2:803361/150‑1 (MQ=255)
             aTGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTac                                                                                                                                     <  1:410164/150‑1 (MQ=255)
              tGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTacc                                                                                                                                    >  1:1081838/1‑150 (MQ=255)
               ggAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTacca                                                                                                                                   >  1:15660/1‑150 (MQ=255)
                gAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTaccac                                                                                                                                  <  1:404649/150‑1 (MQ=255)
                   cGAAGATGCTGCGGTAATGGATAAGTTTAATATCACCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAg                                                                                                                               >  2:864961/1‑150 (MQ=255)
                      aGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGa                                                                                                                            >  2:717775/1‑150 (MQ=255)
                      aGATGCGGCGGTCATGGATAAGTTTAATATCCCCCGTCATAAGCTGTTTGATGTTCAAACGCCTGGCACCGGCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTACCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGa                                                                                                                            <  2:59620/150‑1 (MQ=255)
                       gATGCTTCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCTGCAGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                                                                           <  1:967843/150‑1 (MQ=255)
                       gATGCTGCGGTAATTGAAAAGTTTAATATCCCCCGTCATAAGCAGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGTCATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                                                                           <  2:235366/150‑1 (MQ=255)
                       gATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                                                                           >  2:1095387/1‑150 (MQ=255)
                       gATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                                                                           <  2:369925/150‑1 (MQ=255)
                       gATCCTGCGGTAATGGATAAGGTTTATATCCCCCCGCATCTGCGGTTTGATGTGCAAATGACTGGCACCGTCCTAGGACATATCACACCCCAAGGCGCACTGGCGACACATTTCCCGGCCGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                                                                           <  2:233673/150‑1 (MQ=255)
                        aTGCAGCGGTAATGGATAAGTTTAAAATCCACCGTCATAAGCGGTTTGATGTGCAATTGCCTTGCACCGTCCTCGGACATATCACTCCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATa                                                                                                                          <  2:155778/150‑1 (MQ=255)
                           ctgcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAAc                                                                                                                       >  1:794392/1‑150 (MQ=255)
                               ggTAATGGATAAGTTTAATATCCCCCGTCATATGCTGGTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGt                                                                                                                   <  1:849065/150‑1 (MQ=255)
                               ggTAAAGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCAGCGTCCTCTGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCTGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGt                                                                                                                   <  2:1236614/150‑1 (MQ=255)
                                    tGGATAAGTTTAATATCCCCCGTCATATTCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAg                                                                                                              <  2:888740/150‑1 (MQ=255)
                                          aGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTgg                                                                                                        >  1:473557/1‑150 (MQ=255)
                                           gTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggg                                                                                                       <  1:712359/150‑1 (MQ=255)
                                            tttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTgggg                                                                                                      >  1:483845/1‑150 (MQ=255)
                                            tttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTgggg                                                                                                      <  2:724122/150‑1 (MQ=255)
                                             ttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggggg                                                                                                     <  2:320995/150‑1 (MQ=255)
                                              tAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCCCACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAATCCGGTAGAAGATCTGGGGGc                                                                                                    >  2:946167/1‑150 (MQ=255)
                                                atatCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCg                                                                                                  <  1:1313589/150‑1 (MQ=255)
                                                  atCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGa                                                                                                <  1:953169/150‑1 (MQ=255)
                                                   tCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGAt                                                                                               >  1:738094/1‑150 (MQ=255)
                                                    cccccGTCATATGCTGTTTGCTGGGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACACGCCGCACTGGCGACACATTTCGCGGCAGGACTGCCGGTTGTTTGTCCCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGAtt                                                                                              <  2:639540/150‑1 (MQ=255)
                                                    cccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGAtt                                                                                              >  2:994551/1‑150 (MQ=255)
                                                     ccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTa                                                                                             >  1:1294806/1‑150 (MQ=255)
                                                       ccGTCCTATGCTGATTGATGGGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACt                                                                                           <  2:105779/150‑1 (MQ=255)
                                                        cGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTg                                                                                          <  1:331157/150‑1 (MQ=255)
                                                        cGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTg                                                                                          <  1:293781/150‑1 (MQ=255)
                                                            atatGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACGTGTGATATACCGGTAGAAGCTCTGGGGGCGGGATTACTGgatg                                                                                      <  2:539359/150‑1 (MQ=255)
                                                            atatGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatg                                                                                      <  2:1183661/150‑1 (MQ=255)
                                                               tGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCAGTGGCGCCACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGaggatg                                                                                   >  1:221512/1‑150 (MQ=255)
                                                               tGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGatgatg                                                                                   <  2:1231895/150‑1 (MQ=255)
                                                                 cTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCAAACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCATTGATAAACCGTTAGAAGCTCTGGGGGCCCTATTACTGCATGAGGaa                                                                                 >  2:1262824/1‑150 (MQ=255)
                                                                     ttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcgg                                                                             <  1:701647/150‑1 (MQ=255)
                                                                        ttgtgCACATGCCCGGCACCGTCCTCGGACATATCACAGCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGGTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcggcgg                                                                          <  1:230892/149‑1 (MQ=255)
                                                                          gtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCTGTTGTTTGTCCCACCAGTGATAACCCGGTAGCAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTa                                                                        >  1:219481/1‑150 (MQ=255)
                                                                          gtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTa                                                                        <  2:9588/150‑1 (MQ=255)
                                                                            gCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGCAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAt                                                                      >  1:224230/1‑150 (MQ=255)
                                                                               aaTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTtc                                                                   <  2:34738/150‑1 (MQ=255)
                                                                                aTGCCTGGCACCGTCCTCGGACATATCCCACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTtct                                                                  >  2:958854/1‑150 (MQ=255)
                                                                                  gCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTctct                                                                <  1:560179/150‑1 (MQ=255)
                                                                                   ccTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                               >  2:275682/1‑150 (MQ=255)
                                                                                    cTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTg                                                              <  1:1064588/150‑1 (MQ=255)
                                                                                    cTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTg                                                              <  1:692224/150‑1 (MQ=255)
                                                                                     tGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTgg                                                             <  1:616646/150‑1 (MQ=255)
                                                                                      ggCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGc                                                            >  2:825245/1‑150 (MQ=255)
                                                                                        cACCGTCCTCGGACATATCACACCACAAGCCGCACTGGAGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCAc                                                          <  2:1258622/150‑1 (MQ=255)
                                                                                             tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACa                                                     >  2:111575/1‑150 (MQ=255)
                                                                                              ccTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAt                                                    <  2:1081838/150‑1 (MQ=255)
                                                                                                tCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCg                                                  <  2:15660/150‑1 (MQ=255)
                                                                                                 cGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGc                                                 >  1:104521/1‑150 (MQ=255)
                                                                                                        atCACACCCCAATCCGCACTGGCGACACATTTCCCGGCATGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatg                                          <  2:803872/150‑1 (MQ=255)
                                                                                                        atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatg                                          >  1:268978/1‑150 (MQ=255)
                                                                                                          cacaCCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCCCCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgat                                        >  2:1214262/1‑150 (MQ=255)
                                                                                                            caccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgatga                                      <  1:395230/150‑1 (MQ=255)
                                                                                                            caccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTAACACCAGTGATAAACCGATAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgatga                                      >  1:1137421/1‑150 (MQ=255)
                                                                                                             accacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGaa                                     >  1:1049761/1‑150 (MQ=255)
                                                                                                                     ccGCCCTGGCGCCACATTTCCCGGCAGGACTGCCGGTTGTGTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGCTTACTGGATGATGCAAGGGCGGTAATTTCTCTTGGCACTTACATGGCAGTGATGATGAACGGCAAAg                             >  2:919032/1‑150 (MQ=255)
                                                                                                                           tGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGc                       <  2:794392/150‑1 (MQ=255)
                                                                                                                            ggCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAGATCTCTGGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGcc                      >  1:667410/1‑150 (MQ=255)
                                                                                                                              cGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGCTTACTGGATGATGAAACGGCGGTATTTTCTCTTGGCACTTACATCTCATTGATGATGAACGGCAAATCAATGACGa                    >  2:655434/1‑150 (MQ=255)
                                                                                                                                   caTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAATCCGGTAGAAGCTCTGGGGGCAGGATTACTGGATGATGAAACGGCGGTAATTTCTCTAGGGCCTTACATCGCATTGATGATTAACGGCAAAGCACTGCCGAAAGAt               >  2:212678/1‑150 (MQ=255)
                                                                                                                                   caTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGAt               >  2:636736/1‑150 (MQ=255)
                                                                                                                                   caTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGAt               <  2:1032572/150‑1 (MQ=255)
                                                                                                                                    aTTTCCCGGCAGGACTGCCGGTTGTTTGTACCCCCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTATCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATc              <  2:6377/150‑1 (MQ=255)
                                                                                                                                     tttCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTTATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATcc             >  1:1103283/1‑150 (MQ=255)
                                                                                                                                      ttCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTTATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCa            >  1:111560/1‑150 (MQ=255)
                                                                                                                                      ttCACGGAAGGACTGACGGTTGTTTGTACCACCAGTGATAAACAGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTATCTTTTCAATTACATCACATTGATGATGAACGGCAAAGCACTGCCGAAAGCTCCa            >  2:45976/1‑150 (MQ=255)
                                                                                                                                       tCCCGGCAGGACTGCCGGTTGTATGTACCACCAGTGATAAAGCGGTAGAAGCTGTGGGGGCCGGATTAGTGGATGATGAAACGGCGGTTATTTCTCTTGGGACTTACATCGCTATGTTGATGAACCGCAAAGCCCTGCCGAAAGATCCAg           >  1:795005/1‑150 (MQ=255)
                                                                                                                                       aaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAACCCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAGCGCATTGAGGATGAACTGCAAAGCACTGCCGAAAGATCCAg           >  2:818686/3‑150 (MQ=255)
                                                                                                                                          cGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAg        >  1:570828/1‑150 (MQ=255)
                                                                                                                                             cAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATTAACGGCAAAGCACTGCCTAAAGATCCAGTAGCGt     >  2:666243/1‑150 (MQ=255)
                                                                                                                                             cAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCAGTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGt     >  2:540696/1‑150 (MQ=255)
                                                                                                                                                gACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACt  >  1:489267/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACT  >  NC_000913/2902350‑2902643

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