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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,902,4950AC 7.8% 223.3 / 11.1 77G200G (GGA→GGC) ‡ygcEputative sugar kinase YgcE
Reads supporting (aligned to +/- strand):  ref base A (38/33);  new base C (2/4);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

ATTATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTG  >  NC_000913/2902347‑2902644
                                                                                                                                                    |                                                                                                                                                     
aTTATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGAc                                                                                                                                                      >  1:1144698/1‑150 (MQ=255)
  tataAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTg                                                                                                                                                    >  1:1537988/1‑150 (MQ=255)
    taAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGcc                                                                                                                                                  <  1:556856/150‑1 (MQ=255)
     aaGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCg                                                                                                                                                 <  1:817304/150‑1 (MQ=255)
       gagCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggg                                                                                                                                               <  2:754598/150‑2 (MQ=255)
       gagCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgt                                                                                                                                               <  1:1286268/150‑1 (MQ=255)
        agCTGGGCATGGAGCGAAGATGCTGCGGTATTGGATAAGTTTAATATGCGCCGTCTTATGCTGTTTGATGTGCAAATGCCTGGGACCGTCGTCGGACTTATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCggga                                                                                                                                              <  2:1102662/150‑3 (MQ=255)
         gCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttg                                                                                                                                             <  1:600095/150‑1 (MQ=255)
         gCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttg                                                                                                                                             <  1:1308067/150‑1 (MQ=255)
         gCTGGGCAGGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTAATATTCTGTTTGATGTGCAAATGCATGTCACCGTACTCTGAAATATCACACCACAAGCCCCCCTGGCGACACATTTCCCGGCAGTCCTTCCgggaa                                                                                                                                             <  2:1523390/150‑4 (MQ=255)
           tGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaat                                                                                                                                           <  1:311823/150‑6 (MQ=255)
           tGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                           <  1:969892/150‑1 (MQ=255)
            gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGttgttt                                                                                                                                          >  2:853664/1‑150 (MQ=255)
               cATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTGTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTa                                                                                                                                       <  2:1299206/150‑1 (MQ=255)
                   gAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTaccac                                                                                                                                   >  2:974341/1‑150 (MQ=255)
                   gAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTaccac                                                                                                                                   >  1:1147128/1‑150 (MQ=255)
                    aGCGAAGATGCTGCGGTAATGGATAAGGTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCCACGTCCTCGGCCATATCACACCACAAGCCGCACTGCCGCCACATTTCCCGGCAGGACTGCCGGTTGTTTGTaccacc                                                                                                                                  >  2:516946/1‑150 (MQ=255)
                      cGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAg                                                                                                                                <  2:1064852/150‑1 (MQ=255)
                        aaGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACAGATTTCCGGGCAGGACTGCCGGTTGTTTGTACCACCATTg                                                                                                                              <  2:603741/150‑1 (MQ=255)
                        aaGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTg                                                                                                                              <  1:1104393/150‑1 (MQ=255)
                          gATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                                                                            <  1:117529/150‑1 (MQ=255)
                            tgctgcGGTAATGGATAAGTTTAATATCCAACGTCATATCCTCTTTGATGTGCAAATGCCTGGCACAGTCCTCGGAAATATCACACAACAAGCCGCACTCGCGACACAATTCCAGGCAGGACTGCACGTTGTTTGTAACAACAGTGATaa                                                                                                                          >  2:1135136/1‑150 (MQ=255)
                              ctgcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAAc                                                                                                                        >  1:1255689/1‑150 (MQ=255)
                              ctgcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATAACACACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAAcc                                                                                                                       >  2:901944/1‑150 (MQ=255)
                                gcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCg                                                                                                                      >  2:1295174/1‑150 (MQ=255)
                                  ggTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGt                                                                                                                    >  2:1154159/1‑150 (MQ=255)
                                        ggATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGc                                                                                                              >  2:116648/1‑150 (MQ=255)
                                         gATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGct                                                                                                             >  1:1457319/1‑150 (MQ=255)
                                         gATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCCCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGct                                                                                                             >  2:594147/1‑150 (MQ=255)
                                           tAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGctct                                                                                                           <  2:542442/150‑1 (MQ=255)
                                            aaGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTg                                                                                                          <  2:1255689/150‑1 (MQ=255)
                                                ttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggggg                                                                                                      <  1:1270962/150‑1 (MQ=255)
                                                       cccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCgg                                                                                                  >  2:1027496/1‑147 (MQ=255)
                                                       cccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCgg                                                                                                  <  1:1027496/147‑1 (MQ=255)
                                                           cGTCATATGCTGTTTGAAGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTg                                                                                           <  2:1084462/150‑1 (MQ=255)
                                                            gTCATATGCTGTTTGATGTGCGAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTgg                                                                                          >  1:5896/1‑150 (MQ=255)
                                                            gTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTgg                                                                                          <  2:711770/150‑1 (MQ=255)
                                                              cATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgat                                                                                        <  1:925072/150‑1 (MQ=255)
                                                                 atGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatgat                                                                                     >  1:459802/1‑150 (MQ=255)
                                                                 atGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatgat                                                                                     >  2:1128430/1‑150 (MQ=255)
                                                                 atGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatgat                                                                                     <  1:77909/150‑1 (MQ=255)
                                                                         tGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGATTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcggc                                                                             >  1:1508132/1‑150 (MQ=255)
                                                                          gATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcggcg                                                                            <  1:974341/150‑1 (MQ=255)
                                                                           aTGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcggcgg                                                                           <  1:831191/150‑1 (MQ=255)
                                                                           aTGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTG‑‑GACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTa                                                                         <  1:949757/150‑1 (MQ=255)
                                                                            tgtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGt                                                                          >  2:1149579/1‑150 (MQ=255)
                                                                             gtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTa                                                                         >  2:574584/1‑150 (MQ=255)
                                                                                      ccTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCAACAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                >  1:83345/1‑150 (MQ=255)
                                                                                         ggCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGc                                                             <  2:1392613/150‑1 (MQ=255)
                                                                                           cACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCAc                                                           >  2:606495/1‑150 (MQ=255)
                                                                                              cGTCCTCGGACATATCACACCACTAGCCGCACTGGCGACACATTTCCCGTCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCTGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTa                                                        <  2:1412471/150‑1 (MQ=255)
                                                                                              cGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTa                                                        <  2:1144698/150‑1 (MQ=255)
                                                                                                    cGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGc                                                  >  1:1364731/1‑150 (MQ=255)
                                                                                                    cGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGc                                                  >  1:415463/1‑150 (MQ=255)
                                                                                                     ggACATATCACACCACCAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCa                                                 <  1:1506974/150‑1 (MQ=255)
                                                                                                           atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatg                                           <  1:573200/150‑1 (MQ=255)
                                                                                                             caccCCACAAGCCGCACTGGCGCCACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTACAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGGACTTACAACGCATtgatgat                                         >  1:1110327/1‑150 (MQ=255)
                                                                                                             cacaCCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgat                                         >  2:451427/1‑150 (MQ=255)
                                                                                                             cacaCCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgat                                         >  1:811192/1‑150 (MQ=255)
                                                                                                             cacaCCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgat                                         <  1:1154159/150‑1 (MQ=255)
                                                                                                              acacCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgatg                                        <  1:173618/150‑1 (MQ=255)
                                                                                                               caccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgatga                                       >  2:1173342/1‑150 (MQ=255)
                                                                                                                  cacaAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACg                                    >  1:909597/1‑150 (MQ=255)
                                                                                                                  cacaAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACg                                    <  2:1147128/150‑1 (MQ=255)
                                                                                                                       gCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCaaa                               <  2:415463/150‑1 (MQ=255)
                                                                                                                       gCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCaaa                               >  1:961089/1‑150 (MQ=255)
                                                                                                                                 cGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGa                     >  2:1454244/1‑150 (MQ=255)
                                                                                                                                   acacATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaaa                   >  2:434808/1‑150 (MQ=255)
                                                                                                                                    cacaTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAg                  >  2:276516/1‑150 (MQ=255)
                                                                                                                                    cacaTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAg                  >  2:954685/1‑150 (MQ=255)
                                                                                                                                      caTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCCAAGCACTGCCGAAAGAt                >  2:357587/1‑150 (MQ=255)
                                                                                                                                      aaacaaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGAt                >  1:1089594/7‑150 (MQ=255)
                                                                                                                                       aacaaccGGCAGTCCTGCAGGTTGTTTGTACCACCACTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTAATGGATGATGAAACGGCGGTAATTTATCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGAta               >  2:645922/6‑149 (MQ=255)
                                                                                                                                       aTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATc               <  1:901944/150‑1 (MQ=255)
                                                                                                                                       aTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATc               >  1:391284/1‑150 (MQ=255)
                                                                                                                                            ccGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGAACCAGTa          >  2:1358274/1‑150 (MQ=255)
                                                                                                                                                  ggACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGTGCCTGATTACTGTATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCTCATTGATGATTAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTAc    <  1:1135136/150‑1 (MQ=255)
                                                                                                                                                    aCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTg  <  1:1454244/150‑1 (MQ=255)
                                                                                                                                                    aCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGACGAAAGATCCAGTAGCGTACTg  >  1:487211/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTG  >  NC_000913/2902347‑2902644

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