Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,816,199 G→T 7.2% intergenic (+71/+187) katE → / ← chbG catalase II/chitin disaccharide deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,816,1990GT 7.2% 365.3 / 10.2 113intergenic (+71/+187)katE/chbGcatalase II/chitin disaccharide deacetylase
Reads supporting (aligned to +/- strand):  ref base G (45/57);  new base T (2/7);  total (48/65)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

TATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGG  >  NC_000913/1816050‑1816348
                                                                                                                                                     |                                                                                                                                                     
tATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTg                                                                                                                                                       <  2:367256/150‑1 (MQ=255)
tATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTg                                                                                                                                                       <  1:117803/150‑1 (MQ=255)
tATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTg                                                                                                                                                       <  2:1076093/150‑1 (MQ=255)
    gATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACg                                                                                                                                                   <  1:1625934/150‑1 (MQ=255)
    gATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACg                                                                                                                                                   <  1:1517181/150‑1 (MQ=255)
       gAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                <  1:310260/150‑1 (MQ=255)
        aaCTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                               >  2:1490239/1‑150 (MQ=255)
         aCTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                              >  2:1085080/1‑150 (MQ=255)
                aCGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                       <  1:1049358/150‑1 (MQ=255)
                aCGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                       >  1:1252047/1‑150 (MQ=255)
                 cGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                      <  2:529962/150‑1 (MQ=255)
                 cGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                         >  1:175952/1‑147 (MQ=255)
                 cGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                         <  2:175952/147‑1 (MQ=255)
                 cGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTCACGCCGCATCCGGCATAAAca                                                                                                                                                   >  1:1132874/1‑136 (MQ=255)
                 cGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTCACGCCGCATCCGGCATAAAca                                                                                                                                                   <  2:1132874/137‑2 (MQ=255)
                  gCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATa                                                                                                                                     <  1:207091/150‑1 (MQ=255)
                   cggATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                    <  2:475745/148‑1 (MQ=255)
                   cTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATaa                                                                                                                                    <  2:1336973/150‑2 (MQ=255)
                   cTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATaa                                                                                                                                    <  2:1342765/150‑2 (MQ=255)
                    tGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa                                                                                                                                   >  2:93980/1‑150 (MQ=255)
                        ggCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAc                                                                                                                               <  2:1346069/150‑1 (MQ=255)
                         gcagcaCACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACt                                                                                                                              >  1:1052705/1‑150 (MQ=255)
                          cagcaCACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACtg                                                                                                                             <  1:1233847/150‑1 (MQ=255)
                          cagcaCACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACtg                                                                                                                             >  2:63887/1‑150 (MQ=255)
                             cacaCCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGc                                                                                                                          <  1:1102752/150‑1 (MQ=255)
                               cacCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCaa                                                                                                                        >  1:315856/1‑150 (MQ=255)
                               cacCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCaa                                                                                                                        >  2:222799/1‑150 (MQ=255)
                               cacCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCaa                                                                                                                        >  2:71798/1‑150 (MQ=255)
                                acCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCaaa                                                                                                                       >  2:1074905/1‑150 (MQ=255)
                                   gcgTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTc                                                                                                                    >  2:1173400/1‑150 (MQ=255)
                                    cgTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCa                                                                                                                   <  1:635855/150‑1 (MQ=255)
                                    cgTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCa                                                                                                                   <  1:73900/150‑1 (MQ=255)
                                     gtgtGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCaa                                                                                                                  <  1:661069/150‑1 (MQ=255)
                                      tgtgGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAat                                                                                                                 <  2:1043531/150‑1 (MQ=255)
                                      tgtgGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAat                                                                                                                 <  2:1016551/150‑1 (MQ=255)
                                         ggTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAatata                                                                                                              >  1:1305566/1‑150 (MQ=255)
                                          gTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa                                                                                                                                   <  2:1136772/128‑1 (MQ=255)
                                           tCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTACAAATTCAATATAtt                                                                                                            <  2:1497871/150‑1 (MQ=255)
                                            cACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTg                                                                                                           >  2:1181341/1‑150 (MQ=255)
                                             aCGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGc                                                                                                          <  2:661427/150‑1 (MQ=255)
                                             aCGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGc                                                                                                          <  2:1305566/150‑1 (MQ=255)
                                             aCGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGc                                                                                                          <  2:50322/150‑1 (MQ=255)
                                                     cTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCGTGCAAATTCAATATATTGCAGGAAaca                                                                                                  >  2:46242/1‑150 (MQ=255)
                                                         gATTGACAAAATTCCTGCCTGATGGGAGGCGTCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGCTGTGCAAATTCAATATATTGCAGGAAACACGTa                                                                                              <  2:526140/150‑1 (MQ=255)
                                                         gATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTa                                                                                              <  2:933098/150‑1 (MQ=255)
                                                          aTTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAg                                                                                             <  2:335135/150‑1 (MQ=255)
                                                           ttGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAgg                                                                                            <  2:1565674/150‑1 (MQ=255)
                                                            tGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGc                                                                                           <  1:447838/150‑1 (MQ=255)
                                                               cAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTg                                                                                        <  1:1327627/150‑1 (MQ=255)
                                                                aaaaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGa                                                                                       <  2:1252047/150‑1 (MQ=255)
                                                                 aaaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCTTTAATTATTTATCCGGCATAAAAAAGTTTTAAAATTCTATATATTGTAGGAAAAAAGTAGGctttat                                                                                      <  2:1250289/150‑6 (MQ=255)
                                                                 aaaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGAt                                                                                      >  2:374372/1‑150 (MQ=255)
                                                                 aaaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGAt                                                                                      <  2:1343406/150‑1 (MQ=255)
                                                                     tCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaagc                                                                                  >  1:131241/1‑150 (MQ=255)
                                                                      cctgcctgATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGCTGCGGCGTGAACGCCTTATCCGGCCTACAAATCTGTGCAAATTCAATATATTGCAGGAAACACGTGGGCCTGATaagcg                                                                                 >  2:843230/1‑150 (MQ=255)
                                                                      cctgcctgATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaagcg                                                                                 >  2:1644178/1‑150 (MQ=255)
                                                                       ctgcctgATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaagcga                                                                                <  2:119866/150‑1 (MQ=255)
                                                                             gATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCa                                                                          >  2:1242800/1‑150 (MQ=255)
                                                                                gggAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCa                                                                       <  2:1378625/150‑1 (MQ=255)
                                                                                gggAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCa                                                                       <  1:1085080/150‑1 (MQ=255)
                                                                                  gAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAgg                                                                     <  1:1216679/150‑1 (MQ=255)
                                                                                  gAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAgg                                                                     >  2:667604/1‑150 (MQ=255)
                                                                                     gcggcgCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAg                                                                  <  2:471451/150‑1 (MQ=255)
                                                                                         cgcAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAc                                                                                                                               >  2:766333/1‑85 (MQ=255)
                                                                                         cgcAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAc                                                                                                                               <  1:766333/85‑1 (MQ=255)
                                                                                           cAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGc                                                            >  2:1477502/1‑150 (MQ=255)
                                                                                               tGCGCCGCCTCAAATGATTTACATAGTGCGATTTGTTTATGCCGGATGCGGCGTTAACGCCTTTTCCGGCGTACTAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGttt                                                        <  1:763517/150‑1 (MQ=21)
                                                                                                gcgcCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTg                                                       <  2:1385877/150‑1 (MQ=255)
                                                                                                gcgcCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTg                                                       >  2:1169982/1‑150 (MQ=255)
                                                                                                  gccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTc                                                     <  1:336035/150‑1 (MQ=255)
                                                                                                  gccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTc                                                     <  2:1052705/150‑1 (MQ=255)
                                                                                                  gccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTc                                                     >  2:507884/1‑150 (MQ=255)
                                                                                                    cgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcag                                                   <  2:1437944/150‑1 (MQ=255)
                                                                                                      ccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagca                                                 >  2:819190/1‑150 (MQ=255)
                                                                                                        tCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGt                                               >  1:577960/1‑150 (MQ=255)
                                                                                                         cAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGtc                                              >  2:681956/1‑150 (MQ=255)
                                                                                                          aaaTGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAATCGCATCGGGCAGTTTTGCGTTTGTCAGCAGtct                                             <  2:749868/150‑1 (MQ=255)
                                                                                                          aaaTGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGtct                                             <  2:1228925/150‑1 (MQ=255)
                                                                                                                    cATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTTAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGgcgc                                   >  1:814183/1‑150 (MQ=255)
                                                                                                                      tAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAg                                 <  2:425396/150‑1 (MQ=255)
                                                                                                                      tAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAg                                 >  1:715017/1‑150 (MQ=255)
                                                                                                                       aGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGt                                >  1:1286894/1‑150 (MQ=255)
                                                                                                                        gTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGtt                               >  2:1167657/1‑150 (MQ=255)
                                                                                                                        gTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGtt                               >  2:1546289/1‑150 (MQ=255)
                                                                                                                              tttGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTAcgccg                         >  2:691682/1‑150 (MQ=255)
                                                                                                                                  tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCtt                     >  1:335907/1‑150 (MQ=255)
                                                                                                                                   ttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGTAACTGcgccgcttga                    <  1:40255/150‑5 (MQ=18)
                                                                                                                                     aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGt                  >  2:663244/3‑150 (MQ=21)
                                                                                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAg                <  1:93980/150‑1 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTGATCGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAgg               >  2:269695/1‑150 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTGAACGCCTTATCCGGGCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAgg               <  1:1162291/150‑1 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAgg               >  2:889747/1‑150 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAgg               <  1:1219860/150‑1 (MQ=255)
                                                                                                                                          gtaTGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGaa             <  1:1137672/148‑1 (MQ=255)
                                                                                                                                            aTGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAAtt           <  1:71798/150‑1 (MQ=255)
                                                                                                                                            aTGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAAtt           <  1:906921/150‑1 (MQ=255)
                                                                                                                                                  cgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGc     <  1:1101731/150‑1 (MQ=255)
                                                                                                                                                     gAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGc     >  1:581090/1‑147 (MQ=255)
                                                                                                                                                     gAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGc     <  2:581090/147‑1 (MQ=255)
                                                                                                                                                     gAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCgg  >  1:1502748/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGG  >  NC_000913/1816050‑1816348

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