Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,984,241 T→G 5.3% intergenic (+48/+111) aslB → / ← aslA putative AslA‑specific sulfatase‑maturating enzyme/putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,984,2410TG 5.3% 408.8 / 12.5 114intergenic (+48/+111)aslB/aslAputative AslA‑specific sulfatase‑maturating enzyme/putative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base T (54/54);  new base G (3/3);  total (57/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 120 total aligned reads displayed.

GCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTT  >  NC_000913/3984093‑3984389
                                                                                                                                                    |                                                                                                                                                    
gCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCt                                                                                                                                                      <  1:1763089/149‑1 (MQ=255)
   aTCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                   >  1:529770/1‑149 (MQ=255)
     cTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                                 <  4:345924/149‑1 (MQ=255)
       aCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                               <  1:731899/149‑1 (MQ=255)
       aCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGAGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                               >  2:1772046/1‑149 (MQ=255)
        ccGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                              >  1:1230288/1‑149 (MQ=255)
           ccATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                           >  2:1272636/1‑149 (MQ=255)
            cATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                                          <  2:529770/149‑1 (MQ=255)
            cATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                                          <  4:230236/149‑1 (MQ=255)
            cATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                                          <  1:1081734/149‑1 (MQ=255)
              tatCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAt                                                                                                                                        <  4:1745298/149‑1 (MQ=255)
                  ttAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTATGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                              >  1:305014/1‑139 (MQ=255)
                    aaaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGa                                                                                                                                      >  1:509501/1‑145 (MQ=255)
                    aaaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGa                                                                                                                                      <  2:509501/145‑1 (MQ=255)
                    aaaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGa                                                                                                                                      <  4:598318/145‑1 (MQ=255)
                    aaaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGa                                                                                                                                      >  3:598318/1‑145 (MQ=255)
                    aaaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCAGCCTACagaacgct                                                                                                                                  >  1:322426/1‑142 (MQ=255)
                     aaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTg                                                                                                                                 >  1:812352/1‑149 (MQ=255)
                             gCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCtt                                                                                                                                                     <  1:256906/121‑1 (MQ=255)
                             gCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCtt                                                                                                                                                     >  2:256906/1‑121 (MQ=255)
                             gCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCa                                                                                                                         >  1:656601/1‑149 (MQ=255)
                              cTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCaa                                                                                                                        >  4:394/1‑149 (MQ=255)
                              cTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCaa                                                                                                                        <  3:394/149‑1 (MQ=255)
                               tGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCaaa                                                                                                                             <  1:1181994/143‑1 (MQ=255)
                               tGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCaaa                                                                                                                             >  2:1181994/1‑143 (MQ=255)
                                 aTTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCaaa                                                                                                                             <  4:1834479/141‑1 (MQ=255)
                                 aTTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCCGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGACCGTGCaaa                                                                                                                             >  3:1834479/1‑141 (MQ=255)
                                       tGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCaaa                                                                                                                             <  2:948799/135‑1 (MQ=255)
                                       tGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCa                                                                                                               <  4:464622/149‑1 (MQ=255)
                                         gcgcACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCaa                                                                                                                        >  1:1402827/1‑138 (MQ=255)
                                         gcgcACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCaa                                                                                                                        <  2:1402827/138‑1 (MQ=255)
                                         gcgcACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGc                                                                                                             <  4:1123484/149‑1 (MQ=255)
                                                 tCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCaa                                                                                                                        <  1:1361099/130‑1 (MQ=255)
                                                 tCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCaa                                                                                                                        >  2:1361099/1‑130 (MQ=255)
                                                 tCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAAtt                                                                                                                  <  3:1479255/136‑1 (MQ=255)
                                                 tCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAAtt                                                                                                                  >  4:1479255/1‑136 (MQ=255)
                                                 tCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCa                                                                                                               <  1:1145681/139‑1 (MQ=255)
                                                 tCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCa                                                                                                               >  2:1145681/1‑139 (MQ=255)
                                                           aGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGc                                                                                                                                >  1:413316/1‑112 (MQ=255)
                                                           aGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGc                                                                                                                                <  2:413316/112‑1 (MQ=255)
                                                           aGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATaa                                                                                           >  1:1162236/1‑149 (MQ=255)
                                                             cGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGa                                                                                         <  3:970223/149‑1 (MQ=255)
                                                             cGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCCGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGa                                                                                         <  2:812352/149‑1 (MQ=255)
                                                                  ttATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATg                                                                                       <  3:454677/146‑1 (MQ=255)
                                                                  ttATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGt                                                                                    <  2:1230288/149‑1 (MQ=255)
                                                                      gcatgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgc                                                                                <  1:709645/149‑1 (MQ=255)
                                                                        atgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCTGGCCTACATGATCGTGCAAATTCAATaaa                                                                                                                    >  4:379931/1‑111 (MQ=255)
                                                                        atgcGCAGTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCTGGCCTACATGATCGTGCAAATTCAATaaa                                                                                                                    <  3:379931/111‑1 (MQ=255)
                                                                           cgcATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcgc                                                                           <  1:1272636/149‑1 (MQ=255)
                                                                              aTTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATc                                                                                                                                    >  3:863747/1‑89 (MQ=255)
                                                                              aTTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATc                                                                                                                                    <  4:863747/89‑1 (MQ=255)
                                                                               tttGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAgg                                                                                                   <  1:1633140/121‑1 (MQ=255)
                                                                                    tggtggTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTc                                                                                                                          <  2:1548746/93‑1 (MQ=255)
                                                                                    tggtggTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTc                                                                                                                          >  1:1548746/1‑93 (MQ=255)
                                                                                     ggtggtGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCaaa                                                                 <  1:1772046/149‑1 (MQ=255)
                                                                                      gtggtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCt                                                                                                 <  2:1049258/116‑1 (MQ=255)
                                                                                      gtggtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCt                                                                                                 >  1:1049258/1‑116 (MQ=255)
                                                                                         gtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                             <  4:1037866/69‑1 (MQ=255)
                                                                                         gtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                             >  3:1037866/1‑69 (MQ=255)
                                                                                          tgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATAAGGCGTTTACGCCGCATCCGGca                                                                                                                                             <  4:58395/68‑2 (MQ=17)
                                                                                              taagtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcat                                                                              >  4:251896/1‑127 (MQ=255)
                                                                                              taagtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcat                                                                              <  3:251896/127‑1 (MQ=255)
                                                                                                 gtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGc                                                                                                             <  2:751795/93‑1 (MQ=255)
                                                                                                 gtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGc                                                                                                             >  1:751795/1‑93 (MQ=255)
                                                                                                          ggcggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTgg                                                                                               >  3:1711602/1‑98 (MQ=255)
                                                                                                             ggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCgg                                         >  3:408340/1‑149 (MQ=255)
                                                                                                                   cTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGt                                   <  2:656601/149‑1 (MQ=255)
                                                                                                                        gTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGAt                              <  4:408340/149‑1 (MQ=255)
                                                                                                                         ttGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatc                                                                             <  2:1765088/101‑1 (MQ=255)
                                                                                                                             ttGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCg                                          <  1:1380684/132‑1 (MQ=255)
                                                                                                                             ttGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCg                                          >  2:1380684/1‑132 (MQ=255)
                                                                                                                               gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTACAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGa                                                                                              <  1:178280/78‑1 (MQ=255)
                                                                                                                               gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTACAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGa                                                                                              >  2:178280/1‑78 (MQ=255)
                                                                                                                               gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGAt                                                          >  2:1349277/1‑114 (MQ=37)
                                                                                                                               gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGAt                                                          >  4:1889741/1‑114 (MQ=37)
                                                                                                                               gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGAt                                                          <  3:1889741/114‑1 (MQ=37)
                                                                                                                               gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGAt                                                          <  1:1349277/114‑1 (MQ=37)
                                                                                                                                   gATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGt                                                                                                                                  >  2:400048/1‑38 (MQ=255)
                                                                                                                                        ggcgTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCgtg              >  4:217999/1‑149 (MQ=255)
                                                                                                                                         gcgTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCgtgg             >  1:1434713/1‑149 (MQ=255)
                                                                                                                                           gTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCgtggtg           >  2:1028375/1‑149 (MQ=255)
                                                                                                                                            tAAACGCCTTATCCGGCCTACATGATCGTGCAAAtt                                                                                                                           <  3:160379/36‑1 (MQ=255)
                                                                                                                                            tAAACGCCTTATCCGGCCTACATGATCGTGCAAAtt                                                                                                                           >  4:160379/1‑36 (MQ=255)
                                                                                                                                            tAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCa                                                             <  2:711555/98‑1 (MQ=255)
                                                                                                                                            tAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCa                                                             >  1:711555/1‑98 (MQ=255)
                                                                                                                                               aCGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGata       >  3:584549/1‑149 (MQ=255)
                                                                                                                                                 gCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCa                                                             >  2:1130331/1‑93 (MQ=255)
                                                                                                                                                 gCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCa                                                             <  1:1130331/93‑1 (MQ=255)
                                                                                                                                                 gCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGt                                   <  2:165117/119‑1 (MQ=255)
                                                                                                                                                 gCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGt                                   >  1:165117/1‑119 (MQ=255)
                                                                                                                                                 gCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATAtt     >  4:1709745/1‑149 (MQ=255)
                                                                                                                                                 gCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATAtt     >  2:1209326/1‑149 (MQ=255)
                                                                                                                                                 gCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATAtt     >  1:489540/1‑149 (MQ=255)
                                                                                                                                                   cTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCaaa                                                                 >  2:1620995/1‑87 (MQ=255)
                                                                                                                                                   cTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCaaa                                                                 <  1:1620995/87‑1 (MQ=255)
                                                                                                                                                   cTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTaa                          <  1:1007384/126‑1 (MQ=255)
                                                                                                                                                   cTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTaa                          >  2:1007384/1‑126 (MQ=255)
                                                                                                                                                    ttATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCg                                                                                     >  3:946184/1‑66 (MQ=255)
                                                                                                                                                    ttATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCg                                                                                     <  4:946184/66‑1 (MQ=255)
                                                                                                                                                    ttATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATAttttt  >  1:1409458/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTT  >  NC_000913/3984093‑3984389

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: