Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,816,199 G→T 5.8% intergenic (+71/+187) 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,1990GT 5.8% 364.0 / 9.2 121intergenic (+71/+187)katE/chbGcatalase HPII, heme d‑containing/chito‑oligosaccharide deacetylase
Reads supporting (aligned to +/- strand):  ref base G (47/67);  new base T (4/3);  total (51/70)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 121 total aligned reads displayed.

ATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCC  >  NC_000913/1816051‑1816346
                                                                                                                                                    |                                                                                                                                                   
aTGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGa                                                                                                                                                    <  2:857967/150‑1 (MQ=255)
    aTGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                                                <  2:1129655/150‑1 (MQ=255)
    aTGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                                                <  1:607563/150‑1 (MQ=255)
     tGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                                               >  2:159288/1‑150 (MQ=255)
     tGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                                               >  1:603889/1‑150 (MQ=255)
     tGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                                               <  1:420130/150‑1 (MQ=255)
       aaCTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                             >  1:1010381/1‑150 (MQ=255)
       aaCTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                             <  1:474496/150‑1 (MQ=255)
           gCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                         >  1:969400/1‑150 (MQ=255)
           gCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                         <  2:73097/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)
                cGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGcat                                                                                                                                    <  2:783405/150‑3 (MQ=255)
                 gcggATGGCAGCACCCCGCGTGGGGGCACGCATTCCTAAGATCGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATa                                                                                                                                   <  1:1190821/147‑1 (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)
                            cacaCCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGAGGCGGCGCAAATGCCCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCCGATGCGGCGTGAACGCCTAAACCGGCCATCAAAACTGTGcc                                                                                                                       >  2:1329556/1‑149 (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)
                                     tgtgGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAat                                                                                                               >  1:213519/1‑150 (MQ=255)
                                     tgtgGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAat                                                                                                               <  1:928215/150‑1 (MQ=255)
                                      gtgtTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAata                                                                                                              <  1:1365431/150‑1 (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)
                                              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)
                                                         aTTGACAAAATTACTGCCTGATGGGAGGAGGCGCAATTGCGCCGCCTCAAATGATTTTCATAGTGCGCTTTGTTTATGCAGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAg                                                                                           <  2:1591430/150‑1 (MQ=255)
                                                            gACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTTTTTATGCCGGATGCGGCGTGAACCCCTTTTCCGGCCTAAAAAAGTGTGCAAATTGAATTTTTTGGAGGAAACACGTAGGcc                                                                                        >  2:1516813/1‑150 (MQ=37)
                                                             aaaaaaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCt                                                                                       <  1:1503467/148‑1 (MQ=255)
                                                              cAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTg                                                                                      <  2:1408394/150‑1 (MQ=255)
                                                                 aaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATa                                                                                   <  2:1653890/150‑1 (MQ=255)
                                                                 aaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATa                                                                                   <  2:1647310/150‑1 (MQ=255)
                                                                 aaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATa                                                                                   <  1:170468/150‑1 (MQ=255)
                                                                   ttCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaag                                                                                 >  1:1165110/1‑150 (MQ=255)
                                                                   ttCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTACCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaag                                                                                 <  1:1635603/150‑1 (MQ=255)
                                                                    tCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaagc                                                                                <  1:1181983/150‑1 (MQ=255)
                                                                     cctgcctgATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaagcg                                                                               <  2:799143/150‑1 (MQ=255)
                                                                     cctgcctgATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaagcg                                                                               <  1:439153/150‑1 (MQ=255)
                                                                        gcctgATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaagcgaag                                                                            <  2:235079/150‑1 (MQ=255)
                                                                              tGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATc                                                                      <  2:1545859/150‑1 (MQ=255)
                                                                               gggAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCa                                                                     <  2:811623/150‑1 (MQ=255)
                                                                                ggAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAg                                                                    >  1:333155/1‑150 (MQ=255)
                                                                                  aGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGc                                                                  <  2:946745/150‑1 (MQ=255)
                                                                                   ggcggcGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCa                                                                 >  2:1519111/1‑150 (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)
                                                                                               gcgcCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTg                                                     <  1:1137028/150‑1 (MQ=255)
                                                                                                cgccgcCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGt                                                    >  1:378105/1‑150 (MQ=255)
                                                                                                  ccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTca                                                  <  1:990732/150‑1 (MQ=255)
                                                                                                  ccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTca                                                  <  1:11314/150‑1 (MQ=255)
                                                                                                  ccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTca                                                  <  2:927885/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:541276/150‑1 (MQ=255)
                                                                                                    gccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagc                                                <  1:300317/150‑1 (MQ=255)
                                                                                                      cTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagcag                                              <  2:257322/150‑1 (MQ=255)
                                                                                                      cTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagcag                                              >  1:1617317/1‑150 (MQ=255)
                                                                                                      cTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagcag                                              <  2:454998/150‑1 (MQ=255)
                                                                                                        cAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGtc                                            >  2:1682546/1‑150 (MQ=255)
                                                                                                        cAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGtc                                            >  2:1418753/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:891543/150‑1 (MQ=255)
                                                                                                                     tAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAg                               >  2:546153/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)
                                                                                                                                     ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTGTGTCAGCATTCTCAAGCGGCGCAGTTACGCCGCCTTTGTa               >  2:1124735/2‑150 (MQ=21)
                                                                                                                                          gataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAAt          >  2:638444/6‑150 (MQ=25)
                                                                                                                                              cggcgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTCTGTAGGAATTAAt      >  2:784368/1‑150 (MQ=255)
                                                                                                                                               ggcgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATc     >  2:1177173/1‑150 (MQ=255)
                                                                                                                                               ggcgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATc     >  1:546663/1‑150 (MQ=255)
                                                                                                                                                gcgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCg    <  2:965949/150‑1 (MQ=255)
                                                                                                                                                  gTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGcc  <  1:209547/150‑1 (MQ=255)
                                                                                                                                                  gTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGcc  <  2:333155/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCC  >  NC_000913/1816051‑1816346

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: