Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,816,207 T→C 6.1% intergenic (+79/+179) katE → / ← chbG catalase HPII, heme d‑containing/chito‑oligosaccharide deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,816,2070TC 6.1% 346.3 / 9.4 116intergenic (+79/+179)katE/chbGcatalase HPII, heme d‑containing/chito‑oligosaccharide deacetylase
Reads supporting (aligned to +/- strand):  ref base T (46/62);  new base C (4/3);  total (51/65)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 124 total aligned reads displayed.

AACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGG  >  NC_000913/1816058‑1816356
                                                                                                                                                     |                                                                                                                                                     
aaCTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                       >  1:1010381/1‑150 (MQ=255)
aaCTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                       <  1:474496/150‑1 (MQ=255)
     cTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                  >  2:508727/1‑150 (MQ=255)
      tAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                                 <  2:1161369/150‑1 (MQ=255)
       aaCGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                                >  2:1137028/1‑150 (MQ=255)
       aaCGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                                >  2:634082/1‑150 (MQ=255)
        gccgcatccggcatgaacaaagcacactttgccaacaatctaaatctgtttaattgccggatgtcagacatccggcaattaaacagatttagattgttggcaaAGTGTGCTTTGTTCATGCCGGATGAGGCGTGAACGCCTTATCCGGcc                                                                                                                                               <  2:1578400‑M1/47‑1 (MQ=255)
         cGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                              <  2:844303/150‑1 (MQ=255)
         cGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                              <  2:1155232/150‑1 (MQ=255)
         cGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGcat                                                                                                                                              <  2:783405/150‑3 (MQ=255)
           cTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCGCAAATGATTTAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                            >  1:1354567/1‑150 (MQ=255)
            catccggcatgaacaaagcacactttgccaacaatctaaatctgtttaattgccggatgtcagacatacggcaattaaacagatttagattgttggcaaAGTGTGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGgccttca                                                                                                                                           <  1:1428778‑M1/51‑4 (MQ=255)
               tGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCCAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                                                                        >  1:865552/1‑150 (MQ=255)
                 gcagcaCACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAGTTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACt                                                                                                                                      >  2:443781/1‑150 (MQ=255)
                  gcatgaacaaagcacactttgccaacaatctaaatctgtttaattgccggatgtcagacatccggcaattaaacagatttagattgttggcaaAGTGTGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCctacaaatcct                                                                                                                                     <  2:910502‑M1/57‑5 (MQ=255)
                  cagcaCACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACtg                                                                                                                                     <  2:1095686/150‑1 (MQ=255)
                  cagcaCACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCCGCGCCATTTCGCCGCCTAAAATTATTTAAAATGTGCGCTTTTTTTATGCCGGATTCGGCGTGAAAGCCATATACCGGCAAAAAAACtg                                                                                                                                     >  1:1550282/1‑150 (MQ=255)
                     cacaCCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGc                                                                                                                                  <  2:138987/150‑1 (MQ=255)
                      actcctcgTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGTGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCTCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCa                                                                                                                                 <  1:1191940/144‑1 (MQ=255)
                         ccGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCAGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAAt                                                                                                                              >  1:11710/1‑150 (MQ=255)
                          cgcgTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCTGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCTAAtt                                                                                                                             >  1:268419/1‑150 (MQ=255)
                          cgcgTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAAtt                                                                                                                             >  2:11314/1‑150 (MQ=255)
                                  gTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATTCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAatatat                                                                                                                     >  1:1443303/1‑150 (MQ=255)
                                    cACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTg                                                                                                                   <  2:970970/150‑1 (MQ=255)
                                    cACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTAAATATAttt                                                                                                                   >  2:1650024/1‑149 (MQ=255)
                                      cGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCa                                                                                                                 >  2:135739/1‑150 (MQ=255)
                                      cGCATTCCTAAGATTGACAAAATTCCTGCCTGAGGGGAGGCGGCGCAAATGCGCCGCCGCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCa                                                                                                                 <  1:1302453/150‑1 (MQ=255)
                                       gCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAg                                                                                                                <  2:1444732/150‑1 (MQ=255)
                                         aTTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGa                                                                                                              <  2:932935/150‑1 (MQ=255)
                                             cTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAaca                                                                                                          >  1:965949/1‑150 (MQ=255)
                                                  aTTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAg                                                                                                     <  2:1692832/150‑1 (MQ=255)
                                                     gtcAAAATTCCTGCCTGATGGGAGTCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGcc                                                                                                  <  1:335041/148‑1 (MQ=255)
                                                     gACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTTTTTATGCCGGATGCGGCGTGAACCCCTTTTCCGGCCTAAAAAAGTGTGCAAATTGAATTTTTTGGAGGAAACACGTAGGcc                                                                                                  >  2:1516813/1‑150 (MQ=37)
                                                      aaaaaaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCt                                                                                                 <  1:1503467/148‑1 (MQ=255)
                                                       cAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTg                                                                                                <  2:1408394/150‑1 (MQ=255)
                                                         aaaTTCCTGCCTGTTGGGTGGCGGCGCAATTGCGCGGAATCAAATGTTTTAGATAGTGCGCTTTGTTTTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAATACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGAt                                                                                              <  2:670715/150‑1 (MQ=21)
                                                          aaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATa                                                                                             <  2:1647310/150‑1 (MQ=255)
                                                          aaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATa                                                                                             <  1:170468/150‑1 (MQ=255)
                                                          aaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATa                                                                                             <  1:228952/150‑1 (MQ=255)
                                                          aaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATa                                                                                             <  2:1653890/150‑1 (MQ=255)
                                                            ttCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaag                                                                                           >  1:1165110/1‑150 (MQ=255)
                                                             tCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaagc                                                                                          <  1:1181983/150‑1 (MQ=255)
                                                              cctgcctgATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaagcg                                                                                         <  2:799143/150‑1 (MQ=255)
                                                              catgcctgATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaagcg                                                                                         >  1:946745/3‑150 (MQ=255)
                                                               ctgcctgATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCATACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaagcga                                                                                        <  2:1579425/150‑1 (MQ=255)
                                                                 gcctgATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaagcgaag                                                                                      <  2:235079/150‑1 (MQ=255)
                                                                       tGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATc                                                                                <  2:1545859/150‑1 (MQ=255)
                                                                        gggAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCa                                                                               <  2:286354/150‑1 (MQ=255)
                                                                        gggAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCa                                                                               <  2:811623/150‑1 (MQ=255)
                                                                          gAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCTCATCAgg                                                                             >  2:1492858/1‑150 (MQ=255)
                                                                           aGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGc                                                                            <  2:946745/150‑1 (MQ=255)
                                                                                gcgcAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGttt                                                                       >  1:950679/1‑150 (MQ=255)
                                                                                 cgcAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGtttt                                                                      >  1:1008118/1‑150 (MQ=255)
                                                                                 cgcAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGtttt                                                                      >  2:661994/1‑150 (MQ=255)
                                                                                       tGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAATTTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGttt                                                                >  1:343973/1‑150 (MQ=255)
                                                                                       tGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCGTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGttt                                                                <  1:820549/150‑1 (MQ=255)
                                                                                         cgccgcCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGt                                                              >  1:378105/1‑150 (MQ=255)
                                                                                           ccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTca                                                            <  1:11314/150‑1 (MQ=255)
                                                                                           ccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTca                                                            <  2:927885/150‑1 (MQ=255)
                                                                                           ccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTca                                                            <  1:990732/150‑1 (MQ=255)
                                                                                             gccTCAAATGATTTACATATTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACTCGTATGCGTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagc                                                          <  2:377690/150‑1 (MQ=255)
                                                                                             gccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCTGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCAACAGGCAGTTTTGCGTTTGTcagc                                                          >  1:205494/1‑150 (MQ=255)
                                                                                             gccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagc                                                          <  1:159288/150‑1 (MQ=255)
                                                                                             gccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagc                                                          <  1:300317/150‑1 (MQ=255)
                                                                                               cTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagcag                                                        <  1:443781/150‑1 (MQ=255)
                                                                                               cTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagcag                                                        >  1:1617317/1‑150 (MQ=255)
                                                                                                 cAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGtc                                                      >  2:1682546/1‑150 (MQ=255)
                                                                                                    aTGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCa                                                   <  2:1450159/150‑1 (MQ=255)
                                                                                                         ttACATAGTGATCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAgcgg                                              <  1:1118548/150‑1 (MQ=255)
                                                                                                              tAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAACCTGTGCAAATTCAATATATTGCAGGAAACACGTCGGCCTGATAAGCGCCGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAg                                         >  1:435309/1‑150 (MQ=255)
                                                                                                              tAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAg                                         >  2:546153/1‑150 (MQ=255)
                                                                                                              tAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAg                                         <  2:891543/150‑1 (MQ=255)
                                                                                                                gTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGTGGCGCAGtt                                       >  1:751980/1‑150 (MQ=255)
                                                                                                                gTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGtt                                       >  1:1151520/1‑150 (MQ=255)
                                                                                                                 tGCGCTTTGTGTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTa                                      <  2:187783/150‑1 (MQ=255)
                                                                                                                        tGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACgccgcc                               >  1:939532/1‑150 (MQ=255)
                                                                                                                          tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCtt                             >  1:153309/1‑150 (MQ=255)
                                                                                                                           ttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGTAACTGcgccgcttga                            <  2:1354567/150‑5 (MQ=18)
                                                                                                                            taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAACACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCACGCGGCGCAGTTACGCCGCCTTTg                           >  1:322399/4‑150 (MQ=21)
                                                                                                                            taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTg                           >  1:327218/4‑150 (MQ=21)
                                                                                                                            tATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGTAACTGcgccgcttgag                           <  2:1010381/150‑6 (MQ=18)
                                                                                                                            tATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTg                           <  2:1550282/150‑1 (MQ=255)
                                                                                                                             aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGt                          >  2:467769/1‑150 (MQ=255)
                                                                                                                                   gataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAAt                    >  2:638444/6‑150 (MQ=25)
                                                                                                                                    aTGCGGCGTGAACGCCTTATCAGGCCTAGAACACAGTGCAAATGCAATATATTGCAGTAAACCCGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGGCTCAAGCGGCGCAGGTACGCCGCCTTTGTAGGAAtt                   <  2:1066151/150‑1 (MQ=21)
                                                                                                                                       cggcgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTCTGTAGGAATTAAt                >  2:784368/1‑150 (MQ=255)
                                                                                                                                        ggcgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATc               >  1:546663/1‑150 (MQ=255)
                                                                                                                                        ggcgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATc               >  2:1177173/1‑150 (MQ=255)
                                                                                                                                        ggcgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCGTGGTAGGATTTaagg               >  1:1655318/1‑148 (MQ=255)
                                                                                                                                         gcgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCg              <  2:965949/150‑1 (MQ=255)
                                                                                                                                           gTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGcc            <  1:209547/150‑1 (MQ=255)
                                                                                                                                           gTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGcc            <  2:333155/150‑1 (MQ=255)
                                                                                                                                                cGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATg       >  2:553176/1‑150 (MQ=255)
                                                                                                                                                  ccTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGACGCCTTTGTAGGAATTAATCGACGGATGCa     >  2:268632/1‑150 (MQ=255)
                                                                                                                                                    ttATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAg   <  1:1014452/150‑1 (MQ=255)
                                                                                                                                                    ttATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAg   <  1:508727/150‑1 (MQ=255)
                                                                                                                                                    ttATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAg   >  1:456048/1‑150 (MQ=255)
                                                                                                                                                    ttATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACCTAGTCCTGATAAGCGAAGAGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAg   <  1:1269042/150‑1 (MQ=255)
                                                                                                                                                     tATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAgg  <  1:1380763/150‑1 (MQ=255)
                                                                                                                                                     tATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAgg  >  2:1674194/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGG  >  NC_000913/1816058‑1816356

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