Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,331,859 A→T 6.0% pseudogene (2544/4605 nt) yfaS ← pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,331,8590AT 6.0% 282.7 / 6.6 83pseudogene (2544/4605 nt)yfaSpseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (34/44);  new base T (2/3);  total (36/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GGGTTAACGAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCATTTCATC  >  NC_000913/2331711‑2332005
                                                                                                                                                    |                                                                                                                                                  
gggTTAACGAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCa                                                                                                                                                    <  3:686134/149‑1 (MQ=255)
 ggTTAACGAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAg                                                                                                                                                   <  3:808728/149‑1 (MQ=255)
 ggTTAACGAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAg                                                                                                                                                   <  2:755610/149‑1 (MQ=255)
  gTTAACGAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAgg                                                                                                                                                  >  3:301054/1‑149 (MQ=255)
  gTTAACGAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAgg                                                                                                                                                  >  4:110613/1‑149 (MQ=255)
   tcaaCGAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGgc                                                                                                                                                 <  2:372515/147‑1 (MQ=255)
       cGAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCTGGCGCaa                                                                                                                                             <  4:856024/149‑1 (MQ=255)
        gAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAAc                                                                                                                                            <  1:239796/149‑1 (MQ=255)
          aTCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCg                                                                                                                                          >  2:553862/1‑149 (MQ=255)
           tCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCgg                                                                                                                                         <  1:163454/149‑1 (MQ=255)
           tCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCgg                                                                                                                                         >  2:103578/1‑149 (MQ=255)
            cAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCTGGCGCAActaac                                                                                                                                        <  2:808428/149‑5 (MQ=255)
            cAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGc                                                                                                                                        <  2:838025/149‑1 (MQ=255)
                 aTCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCt                                                                                                                                   <  2:167718/149‑1 (MQ=255)
                     ggAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCTGGCGCaa                                                                                                                                             <  4:160869/135‑1 (MQ=255)
                     ggAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCTGGCGCaa                                                                                                                                             >  3:160869/1‑135 (MQ=255)
                      gAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGgaga                                                                                                                              <  3:529057/149‑1 (MQ=255)
                            gAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCt                                                                                                                        <  2:654094/149‑1 (MQ=255)
                            gAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCt                                                                                                                        >  4:721307/1‑149 (MQ=255)
                            gAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCt                                                                                                                        <  3:528149/149‑1 (MQ=255)
                                  cGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGt                                                                                                                  <  3:505032/149‑1 (MQ=255)
                                   gCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTa                                                                                                                 <  3:721307/149‑1 (MQ=255)
                                    cTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAg                                                                                                                <  4:244576/149‑1 (MQ=255)
                                        gCCTTGTTTATCGGTTGTGAGTGACGGCATCCCTGCCGCGGTATCCACCTCTTCACGCCGTGAACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGa                                                                                                            <  2:355646/149‑1 (MQ=255)
                                        gCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGa                                                                                                            <  3:420069/149‑1 (MQ=255)
                                         ccTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGAt                                                                                                           >  3:26872/1‑149 (MQ=255)
                                             gTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAAc                                                                                                       <  4:102244/149‑1 (MQ=255)
                                             gTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAAc                                                                                                       <  3:293538/149‑1 (MQ=255)
                                             gTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAAc                                                                                                       <  1:870015/149‑1 (MQ=255)
                                               ttATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGTCGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGa                                                                                                     <  2:71058/149‑1 (MQ=255)
                                               ttATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCCCTGGAGTGCGCCTGGTAGTAGCTGATAAACGa                                                                                                     >  2:407335/1‑149 (MQ=255)
                                                 aTCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACa                                                                                                   <  1:53168/149‑1 (MQ=255)
                                                    ggTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGAc                                                                                                <  1:690212/149‑1 (MQ=255)
                                                    ggTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGAc                                                                                                <  3:110613/149‑1 (MQ=255)
                                                     gTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGgaga                                                                                                                              <  4:620002/118‑1 (MQ=255)
                                                     gTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGgaga                                                                                                                              >  3:620002/1‑118 (MQ=255)
                                                        gtgagtgaCGGCATCCATGCCGCGGTATCCACCTCTTCACGCCATGGACGTTCAAGCATTTTTACTCGCCGCTCGCTCCGGTTAGTTGCGCCAAGCGTAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGa                                                                                            <  4:548761/149‑1 (MQ=255)
                                                          gagtgaCGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGc                                                                                          >  1:267024/1‑149 (MQ=255)
                                                          gagtgaCGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGGTGCGCCAGGCTCAACCGGCTCGCTGGTGAGCGCCTGGTCGTAGCTGTTAAACGACAGACTGGAGc                                                                                          >  1:389629/1‑149 (MQ=255)
                                                             tgaCGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTgg                                                                                       <  4:141992/149‑1 (MQ=255)
                                                                 ggCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACg                                                                                   >  3:470082/1‑149 (MQ=255)
                                                                     tCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGcaca                                                                               >  4:656166/1‑149 (MQ=255)
                                                                      ccATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACAt                                                                              >  3:765768/1‑149 (MQ=255)
                                                                       cATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACAtt                                                                             <  1:468830/149‑1 (MQ=255)
                                                                       cATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACAtt                                                                             <  2:65247/149‑1 (MQ=255)
                                                                        aTGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTg                                                                            >  1:441514/1‑149 (MQ=255)
                                                                           ccGCGGTATCCACCTCTTCACGCCGTCGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTAAGTTGAGCCATGCGCAGCCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGGTAAATGACAGTCTCGAGCTTGTACGCACATTGTTa                                                                         <  4:739568/149‑1 (MQ=255)
                                                                                     cACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCg                                                               >  4:316841/1‑149 (MQ=255)
                                                                                      aCCTCTTCATGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCgg                                                              >  4:760356/1‑149 (MQ=255)
                                                                                      aCCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCgg                                                              >  2:785388/1‑149 (MQ=255)
                                                                                      aCCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCgg                                                              >  2:747520/1‑149 (MQ=255)
                                                                                       ccTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTTGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGa                                                             <  4:313310/149‑1 (MQ=255)
                                                                                       ccTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGa                                                             <  4:26872/149‑1 (MQ=255)
                                                                                        ctctTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGAt                                                            >  2:300373/1‑149 (MQ=255)
                                                                                         tctTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGa                                                             <  4:266603/147‑1 (MQ=255)
                                                                                         tctTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGa                                                             >  3:266603/1‑147 (MQ=255)
                                                                                         tctTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGa                                                             >  1:416641/1‑147 (MQ=255)
                                                                                         tctTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGa                                                             <  2:416641/147‑1 (MQ=255)
                                                                                         tctTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATa                                                           <  2:218663/149‑1 (MQ=255)
                                                                                                       ccGTTCAAGCATAATTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGTCTCGCTGGAGAGCGCCTGGTCGTAGCTTATAAACGACAGACTAGAGCTGGTACGCACATTGTTACGCCTCAGCAGATAGAAAAATTTGCCGa                                             >  4:376561/2‑149 (MQ=255)
                                                                                                               gCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGca                                                                                 <  2:117803/105‑1 (MQ=255)
                                                                                                               gCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGca                                                                                 >  1:117803/1‑105 (MQ=255)
                                                                                                               gCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGgc                                     >  1:139845/1‑149 (MQ=255)
                                                                                                               gCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGgc                                     >  2:134675/1‑149 (MQ=255)
                                                                                                               gCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGgc                                     <  4:765768/149‑1 (MQ=255)
                                                                                                                 aTTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGgcgc                                   >  2:627935/1‑149 (MQ=255)
                                                                                                                   ttttACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGa                                 <  4:323898/149‑1 (MQ=255)
                                                                                                                      tACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGAttt                              <  2:696302/149‑1 (MQ=255)
                                                                                                                      tACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGAttt                              <  3:352895/149‑1 (MQ=255)
                                                                                                                         tCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGTTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTg                           >  3:412715/1‑149 (MQ=255)
                                                                                                                         tCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAAtt                                                   >  2:175284/1‑125 (MQ=255)
                                                                                                                         tCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAAtt                                                   <  1:175284/125‑1 (MQ=255)
                                                                                                                           gccgcTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGt                         >  3:68853/1‑149 (MQ=255)
                                                                                                                           gccgcTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGt                         >  3:88555/1‑149 (MQ=255)
                                                                                                                            ccgcTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGtt                        <  2:302503/149‑1 (MQ=255)
                                                                                                                             cgctcgctGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTg                       >  2:882947/1‑149 (MQ=255)
                                                                                                                             cgctcgctGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTCCCGATATTCGGCGCGATTTCTGGTTg                       >  2:532846/1‑149 (MQ=255)
                                                                                                                               ctcgctGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCa                     <  4:470082/149‑1 (MQ=255)
                                                                                                                                        gTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGa            <  2:594470/149‑1 (MQ=255)
                                                                                                                                        gTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGa            >  4:176135/1‑149 (MQ=255)
                                                                                                                                         ttAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGAt           <  3:656166/149‑1 (MQ=255)
                                                                                                                                            gTTGCGCCTGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCAt        >  4:308220/1‑149 (MQ=255)
                                                                                                                                               gcgcCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCATTTc     >  3:104060/1‑149 (MQ=255)
                                                                                                                                                  ccAGGCGCAACCGGCTCGCTGGAGAGCGACTGGTCGTAGCTGATAAACGACAGACTGGGGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGACCATTTCATc  <  4:759421/149‑1 (MQ=255)
                                                                                                                                                  caaGGCGCATCCGGCTCGCTGGAGAGCGCCTGGTCGTTGCTGATAGACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGTAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCATTTCagc  <  4:205105/147‑3 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGGTTAACGAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCATTTCATC  >  NC_000913/2331711‑2332005

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