Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,816,192 G→A 5.3% intergenic (+64/+194) 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,1920GA 5.3% 361.6 / 5.8 113intergenic (+64/+194)katE/chbGcatalase HPII, heme d‑containing/chito‑oligosaccharide deacetylase
Reads supporting (aligned to +/- strand):  ref base G (41/65);  new base A (3/3);  total (45/68)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

Only 100 of 117 total aligned reads displayed.

GTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATT  >  NC_000913/1816046‑1816339
                                                                                                                                                  |                                                                                                                                                   
tttttttGGATGAACTGCTAACGCTGATGGCAGCACACCGCGGGTGGTCACGAATTCCTAAGATTGACAAAATTCCTGCCTTATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                                                                  <  1:1578418/149‑1 (MQ=255)
tttttATGGATTAACTGCTAACGTTGATGTCAGCACACTGCGTGTGGTCACGCATTCCTAAGATTGTCAAAATTCCTGCCTGTTGGGAGGCGGCGCATTTGCGCGGCCTCAAATTATTTACATAGTGCGCTTTGTTTATGCCGGATgcgt                                                                                                                                                  <  1:1308071/149‑2 (MQ=255)
gTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTCATGCCGGATgcgg                                                                                                                                                  <  1:143775/150‑1 (MQ=255)
  tttATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                                                                <  2:423668/150‑1 (MQ=255)
   ttATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTTGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                                                               <  1:756664/150‑1 (MQ=255)
   ttATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                                                               <  2:577080/150‑1 (MQ=255)
     aTGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGa                                                                                                                                             <  2:857967/150‑1 (MQ=255)
         aTGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                                         <  1:607563/150‑1 (MQ=255)
         aTGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                                         <  2:1129655/150‑1 (MQ=255)
          tGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                                        <  1:420130/150‑1 (MQ=255)
          tGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                                        >  1:603889/1‑150 (MQ=255)
          tGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                                        >  2:159288/1‑150 (MQ=255)
            aaCTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                      >  1:1010381/1‑150 (MQ=255)
                gCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                  >  1:969400/1‑150 (MQ=255)
                gCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                  <  2:73097/150‑1 (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:1155232/150‑1 (MQ=255)
                     cGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                             <  2:844303/150‑1 (MQ=255)
                      gcggATGGCAGCACCCCGCGTGGGGGCACGCATTCCTAAGATCGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATa                                                                                                                            <  1:1190821/147‑1 (MQ=255)
                       cTGATGGCAGCACACCGCGTGTGGTCACACATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATaa                                                                                                                           <  2:823037/150‑2 (MQ=255)
                        catccggcatgaacaaagcacactttgccaacaatctaaatctgtttaattgccggatgtcagacatacggcaattaaacagatttagattgttggcaaAGTGTGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGgccttca                                                                                                                          <  1:1428778‑M1/51‑4 (MQ=255)
                             gcagcaCACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAGTTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACt                                                                                                                     >  2:443781/1‑150 (MQ=255)
                              cagcaCACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACtg                                                                                                                    <  2:1095686/150‑1 (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)
                                          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)
                                                  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)
                                                              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)
                                                                     aaaTTCCTGCCTGTTGGGTGGCGGCGCAATTGCGCGGAATCAAATGTTTTAGATAGTGCGCTTTGTTTTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAATACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGAt                                                                             <  2:670715/150‑1 (MQ=21)
                                                                      aaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATa                                                                            <  1:228952/150‑1 (MQ=255)
                                                                      aaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATa                                                                            <  1:170468/150‑1 (MQ=255)
                                                                      aaTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATa                                                                            <  2:1653890/150‑1 (MQ=255)
                                                                        ttCGTGGCTGATGGGAGGCTGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAACTGTGCTAATTCAATATATTGCAGGTAACACGTAGGCCTGATaag                                                                          <  2:932871/150‑1 (MQ=25)
                                                                        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)
                                                                          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:811623/150‑1 (MQ=255)
                                                                                    gggAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCa                                                              <  2:286354/150‑1 (MQ=255)
                                                                                     ggAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAg                                                             >  1:333155/1‑150 (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)
                                                                                                    gcgcCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTg                                              <  1:1137028/150‑1 (MQ=255)
                                                                                                     cgccgcCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGt                                             >  1:378105/1‑150 (MQ=255)
                                                                                                       ccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTca                                           <  2:927885/150‑1 (MQ=255)
                                                                                                       ccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTca                                           <  1:11314/150‑1 (MQ=255)
                                                                                                         gccTCAAATGATTTACATATTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACTCGTATGCGTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagc                                         <  2:377690/150‑1 (MQ=255)
                                                                                                         gccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCTGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCAACAGGCAGTTTTGCGTTTGTcagc                                         >  1:205494/1‑150 (MQ=255)
                                                                                                         gccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagc                                         <  1:541276/150‑1 (MQ=255)
                                                                                                         gccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagc                                         <  1:159288/150‑1 (MQ=255)
                                                                                                         gccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagc                                         <  1:300317/150‑1 (MQ=255)
                                                                                                           cTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagcag                                       <  1:443781/150‑1 (MQ=255)
                                                                                                           cTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagcag                                       <  2:454998/150‑1 (MQ=255)
                                                                                                           cTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagcag                                       >  1:1617317/1‑150 (MQ=255)
                                                                                                             cAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGtc                                     >  2:1418753/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:891543/150‑1 (MQ=255)
                                                                                                                          tAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAg                        >  2:546153/1‑150 (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)
                                                                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTGTGTCAGCATTCTCAAGCGGCGCAGTTACGCCGCCTTTGTa        >  2:1124735/2‑150 (MQ=21)
                                                                                                                                                aTGCGGCGTGAACGCCTTATCAGGCCTAGAACACAGTGCAAATGCAATATATTGCAGTAAACCCGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGGCTCAAGCGGCGCAGGTACGCCGCCTTTGTAGGAAtt  <  2:1066151/150‑1 (MQ=21)
                                                                                                                                                  |                                                                                                                                                   
GTTTTATGGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATT  >  NC_000913/1816046‑1816339

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: