Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,814,742 C→G 5.4% H292Q (CAC→CAG)  katE → catalase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,814,7420CG 5.4% 312.0 / 7.0 93H292Q (CAC→CAG) katEcatalase II
Reads supporting (aligned to +/- strand):  ref base C (48/39);  new base G (2/3);  total (51/42)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

TTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGAT  >  NC_000913/1814594‑1814888
                                                                                                                                                    |                                                                                                                                                  
tttCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCCTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACt                                                                                                                                                   >  1:698757/1‑150 (MQ=255)
 ttCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTg                                                                                                                                                  <  1:1377075/150‑1 (MQ=255)
  tctcTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCAGTgg                                                                                                                                                 >  2:716876/1‑150 (MQ=255)
    tctGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGTAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGaa                                                                                                                                               <  2:664913/150‑1 (MQ=255)
    tctGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGaa                                                                                                                                               >  2:290720/1‑150 (MQ=255)
      tGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAAc                                                                                                                                             >  1:786749/1‑150 (MQ=255)
       gCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAAcc                                                                                                                                            <  1:1377937/150‑1 (MQ=255)
        cAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCAGTGGAAAcgt                                                                                                                                           <  2:151041/150‑3 (MQ=255)
         aaCCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCAc                                                                                                                                          >  2:436051/1‑150 (MQ=255)
          aCCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACt                                                                                                                                         <  1:939275/150‑1 (MQ=255)
          aCCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACt                                                                                                                                         <  1:945222/150‑1 (MQ=255)
           ccTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCAGTGGAAAcgtaca                                                                                                                                        <  1:1320954/150‑6 (MQ=255)
           acttaaaCTCTTCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCAGCAGTTACCGCACCATGGAAGGCTTCTGTATTCACAACTTCAGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACAACTg                                                                                                                                        <  2:679862/146‑1 (MQ=255)
            cTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCAGTGGAAAcgtacaa                                                                                                                                       <  1:685468/150‑7 (MQ=255)
                    ctGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTaaa                                                                                                                               >  1:422437/1‑150 (MQ=255)
                     tGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAg                                                                                                                              <  1:1542899/150‑1 (MQ=255)
                     tGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCCGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAg                                                                                                                              >  2:893090/1‑150 (MQ=255)
                      gCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGGTTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGc                                                                                                                             >  1:765079/1‑150 (MQ=255)
                       cacaACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGcc                                                                                                                            >  2:1327378/1‑150 (MQ=255)
                          aaCGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCa                                                                                                                         >  2:1606422/1‑150 (MQ=255)
                           aCGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCAc                                                                                                                        <  1:1606422/150‑1 (MQ=255)
                             gTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTc                                                                                                                      >  1:11796/1‑150 (MQ=255)
                                aTGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCGTGATTAATGCCGAAGGGAAGGCAAGGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCATTcagc                                                                                                                   >  1:643914/1‑147 (MQ=255)
                                  gtgGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACAGGCAGGTAAAGCCTCACTCGTTTg                                                                                                                 >  2:585024/1‑150 (MQ=255)
                                    gggCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTc                                                                                                                          <  1:437873/139‑1 (MQ=255)
                                    gggCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTc                                                                                                                          >  2:437873/1‑139 (MQ=255)
                                    gggCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTggg                                                                                                               >  1:285912/1‑150 (MQ=255)
                                    gggCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTggg                                                                                                               <  1:436051/150‑1 (MQ=255)
                                    gggCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTggg                                                                                                               >  1:921387/1‑150 (MQ=255)
                                    gggCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTggg                                                                                                               <  1:1255234/150‑1 (MQ=255)
                                     ggCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGa                                                                                                              <  2:1172373/150‑1 (MQ=255)
                                             cGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAGGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGcac                                                                                                      <  2:1287552/150‑1 (MQ=255)
                                                 tCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACaaaa                                                                                                  <  1:290720/150‑1 (MQ=255)
                                                   gcgGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAAc                                                                                                >  1:1446709/1‑150 (MQ=255)
                                                   gcgGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAAc                                                                                                <  2:833412/150‑1 (MQ=255)
                                                        aTCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCAcc                                                                                           <  1:200460/150‑1 (MQ=255)
                                                          cccccGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCgg                                                                                         <  2:285912/150‑1 (MQ=255)
                                                          cccccGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCgg                                                                                         >  1:1222590/1‑150 (MQ=255)
                                                           ccccGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGa                                                                                        >  1:1375199/1‑150 (MQ=255)
                                                             ccGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCCGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACg                                                                                      <  1:893090/150‑1 (MQ=255)
                                                                  gTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGAcc                                                                                 >  2:1358085/1‑150 (MQ=255)
                                                                  gTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGAcc                                                                                 >  1:626510/1‑150 (MQ=255)
                                                                    tACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCg                                                                               >  1:469510/1‑150 (MQ=255)
                                                                    tACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCg                                                                               <  2:1037533/150‑1 (MQ=255)
                                                                        gCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGa                                                                                        >  2:215194/1‑137 (MQ=255)
                                                                        gCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGa                                                                                        <  1:215194/137‑1 (MQ=255)
                                                                          aCCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTc                                                                         <  2:1504992/150‑1 (MQ=255)
                                                                            cATGGAAGGCTTCGGTCTAGACTGCTGCCGCCTGATTAATGCCGATGGGAAGGCTACGTTTGTACGTTTCCACTCGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCa                                                                       <  2:1393951/150‑1 (MQ=255)
                                                                             aTGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCAc                                                                      <  1:1358085/150‑1 (MQ=255)
                                                                              tGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTTGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCAcc                                                                     <  2:387919/150‑1 (MQ=255)
                                                                              tGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCAcc                                                                     <  1:965874/150‑1 (MQ=255)
                                                                                    gCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCcgcg                                                               <  1:828459/150‑1 (MQ=255)
                                                                                     cTTCGGTAATCACACCTTCCGCCTGATTAATGCCGAAGGGAACGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCCCCGGACGTGCCCCGGACTTCCACCGCCGCGa                                                              >  2:1149531/1‑150 (MQ=255)
                                                                                     cTTCGGAATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGa                                                              >  2:42489/1‑150 (MQ=255)
                                                                                      ttCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAg                                                             <  2:970849/150‑1 (MQ=255)
                                                                                      ttCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGAAGTGACCCGGACTTCCACCGCCGCGAg                                                             <  2:921387/150‑1 (MQ=255)
                                                                                       tCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGt                                                            >  1:149471/1‑150 (MQ=255)
                                                                                        cGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGtt                                                           <  2:11796/150‑1 (MQ=255)
                                                                                         ggTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTtg                                                          >  1:906325/1‑150 (MQ=255)
                                                                                             ttCACACCTTCATCCTTATTAATGACGAAGGTAAGGCAACGTTTTTACGTTTCCACTGGAAACAACTGGAAGGTAAAGACTCACTCGTTTGGGATGAATCACAAAAAGTCACCGGAATTGACCCGGACTTCGACCGCCGCGAGTTGTggg                                                      <  2:324657/150‑1 (MQ=255)
                                                                                              tCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGATTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACCTGACCCGGACTTCCACCGCCGCGAGTTGTGGGa                                                     >  1:1409298/1‑150 (MQ=255)
                                                                                               cacaCCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTTTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAATAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGaa                                                    >  2:1091016/1‑150 (MQ=255)
                                                                                               cacaCCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGaa                                                    >  1:566855/1‑150 (MQ=255)
                                                                                                       ccGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGa                                            >  1:276371/1‑150 (MQ=255)
                                                                                                       ccGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGa                                            <  1:1327378/150‑1 (MQ=255)
                                                                                                         gCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAg                                          >  2:924498/1‑150 (MQ=255)
                                                                                                         gCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAg                                          >  2:309203/1‑150 (MQ=255)
                                                                                                           cTGGTTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCa                                        <  2:765079/150‑1 (MQ=255)
                                                                                                            tGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAg                                       >  2:1277243/1‑150 (MQ=255)
                                                                                                                  aTGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGAtt                                 >  2:939237/1‑150 (MQ=255)
                                                                                                                   tGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGAttt                                >  1:1649961/1‑150 (MQ=255)
                                                                                                                   tGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGAttt                                >  2:1433858/1‑150 (MQ=255)
                                                                                                                   tGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGAttt                                >  2:896242/1‑150 (MQ=255)
                                                                                                                    gCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGAtttt                               >  2:1526485/1‑150 (MQ=255)
                                                                                                                    gCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGAtttt                               <  1:309203/150‑1 (MQ=255)
                                                                                                                    gCCGAAGAGAAGGCAACGATTGTACGTTTACAATGCAAACCACTGGAAGCTAAAGCCTCAATCGTTTGGGATGAAGCACAAAAACTAACCGGACGTGAACCGGACTTACACCCCCGCGATTTGTGGCAAGACAAAGAAGAAGGCGAtttt                               >  2:888395/1‑150 (MQ=255)
                                                                                                                        aagggaaggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCgg                           <  2:267396/150‑1 (MQ=255)
                                                                                                                            gaaggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATa                       >  1:1261423/1‑150 (MQ=255)
                                                                                                                            gaaggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATa                       >  2:1522889/1‑150 (MQ=255)
                                                                                                                              aggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACg                     >  1:268553/1‑150 (MQ=255)
                                                                                                                              aggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACg                     <  2:1375199/150‑1 (MQ=255)
                                                                                                                                gCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGaa                   >  1:428758/1‑150 (MQ=255)
                                                                                                                                    cGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCTCTCGATTGGGATGAAGCAGTAAATGTCACCGGAGGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGATTACGAACTgg               <  1:1457823/150‑1 (MQ=255)
                                                                                                                                       ccagtggTTTCCAGTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCt            >  1:1616784/7‑150 (MQ=255)
                                                                                                                                        tGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCtt           >  2:67730/1‑150 (MQ=255)
                                                                                                                                          tACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCAGTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTcc         <  2:910831/150‑1 (MQ=255)
                                                                                                                                           aCGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCCCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCa        <  1:541510/150‑1 (MQ=255)
                                                                                                                                            cGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTAGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAg       >  2:256432/1‑150 (MQ=255)
                                                                                                                                             gTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGt      >  1:587935/1‑150 (MQ=255)
                                                                                                                                             gTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGt      >  1:1251484/1‑150 (MQ=255)
                                                                                                                                               ttCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTg    <  2:566855/150‑1 (MQ=255)
                                                                                                                                                tCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGa   >  1:20932/1‑150 (MQ=255)
                                                                                                                                                 ccACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGAt  >  1:615910/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTTCTCTGCAACCTGAAACTCTGCACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGAT  >  NC_000913/1814594‑1814888

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