Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,857,696 G→C 7.7% P237A (CCA→GCA)  yidJ ← sulfatase/phosphatase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,857,6960GC 7.7% 403.1 / 21.2 117P237A (CCA→GCA) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (50/58);  new base C (4/5);  total (54/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-01

Only 100 of 117 total aligned reads displayed.

CGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAA  >  NC_000913/3857549‑3857841
                                                                                                                                                   |                                                                                                                                                 
cGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGc                                                                                                                                                  >  1:513053/1‑149 (MQ=255)
 gTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCa                                                                                                                                                 <  4:284504/149‑1 (MQ=255)
   gATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATc                                                                                                                                               <  4:720137/149‑1 (MQ=255)
   gATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTCGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGCCATc                                                                                                                                               <  3:239828/149‑1 (MQ=255)
       ggAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCt                                                                                                                                           >  4:559438/1‑149 (MQ=255)
          ggTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGCCATCGccagtc                                                                                                                                        <  2:472102/149‑5 (MQ=255)
            tataAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGCCATCGccagtc                                                                                                                                        <  1:130539/147‑5 (MQ=255)
            tataAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGCCATCGccagtc                                                                                                                                        >  2:130539/1‑143 (MQ=255)
             ataaataaCCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGccc                                                                                                                                     <  2:300511/149‑1 (MQ=255)
             ataaataaCCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGccc                                                                                                                                     <  2:360048/149‑1 (MQ=255)
              taaataaCCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCa                                                                                                                                    <  3:728766/149‑1 (MQ=255)
               aaataaCCCACGTATTTTCACGTTGCTCTGGCGTTAAGGAATTGATGACCCGTCCGATTTGGTCATCAACAAAGTAATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATTGCATCGCCGGCGCCCAt                                                                                                                                   <  2:119139/149‑1 (MQ=255)
                  taaCCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAg                                                                                                                                <  3:386529/149‑1 (MQ=255)
                   aaCCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGc                                                                                                                               >  2:136829/1‑149 (MQ=255)
                        aCGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGa                                                                                                                          <  1:548877/149‑1 (MQ=255)
                        aCGTATTTGCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGa                                                                                                                          >  3:699883/1‑149 (MQ=255)
                         cGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGAt                                                                                                                         <  2:46677/149‑1 (MQ=255)
                          gTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATg                                                                                                                        >  3:523359/1‑149 (MQ=255)
                                   cGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGttt                                                                                                               <  3:673112/149‑1 (MQ=255)
                                    gTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTg                                                                                                              >  3:5355/1‑149 (MQ=255)
                                    gTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTg                                                                                                              <  2:767546/149‑1 (MQ=255)
                                        ctctGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAg                                                                                                          <  4:151008/149‑1 (MQ=255)
                                         tctGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGc                                                                                                         >  1:292102/1‑149 (MQ=255)
                                         tctGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGCCATCGccagtc                                                                                                                                        <  1:276300/118‑5 (MQ=255)
                                         tctGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGCCATCGccagtc                                                                                                                                        >  2:276300/1‑114 (MQ=255)
                                           tGGCGTTAAGGCATTGATGCCCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCa                                                                                                       <  4:517277/149‑1 (MQ=255)
                                           tGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCa                                                                                                       <  1:685307/149‑1 (MQ=255)
                                            ggCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAg                                                                                                      >  4:267863/1‑149 (MQ=255)
                                               gTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTc                                                                                                   >  2:392440/1‑149 (MQ=255)
                                               gTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATAGCCTGCGCCAATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTc                                                                                                   <  2:396668/149‑1 (MQ=255)
                                                  aaGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATc                                                                                                >  2:361852/1‑149 (MQ=255)
                                                   aGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATcc                                                                                               <  3:240492/149‑1 (MQ=255)
                                                     gCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCtg                                                                                             >  4:380429/1‑149 (MQ=255)
                                                      cATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCtgt                                                                                            <  3:559438/149‑1 (MQ=255)
                                                        ttGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGc                                                                                          <  3:520965/149‑1 (MQ=255)
                                                          gatgaCCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCtt                                                                                        >  1:643486/1‑149 (MQ=255)
                                                           atgaCCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAg                                                                                                      >  3:98167/1‑134 (MQ=255)
                                                              aCCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCATCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTtctc                                                                                    >  2:596586/1‑149 (MQ=255)
                                                                ccGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGc                                                                                  >  3:322836/1‑149 (MQ=255)
                                                                 cGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGcc                                                                                 >  3:237589/1‑149 (MQ=255)
                                                                   tCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCa                                                                               >  2:424832/1‑149 (MQ=255)
                                                                     cGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGc                                                                             <  4:421341/149‑1 (MQ=255)
                                                                       aTTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCtgtg                                                                                           >  3:359804/1‑133 (MQ=255)
                                                                        tttGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCg                                                                          >  3:159567/1‑149 (MQ=255)
                                                                        tttGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCg                                                                          >  3:298661/1‑149 (MQ=255)
                                                                          tGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCgta                                                                        <  2:74339/149‑1 (MQ=255)
                                                                               atcaACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTaaa                                                                   <  4:250688/149‑1 (MQ=255)
                                                                                 caacaaAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTGGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTaaaaa                                                                 <  2:269741/149‑1 (MQ=255)
                                                                                 caacaaAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTaaaaa                                                                 <  4:200989/149‑1 (MQ=255)
                                                                                  aacaaAGTCATTACAGGCAAAATAGAGCGTATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAAt                                                                >  3:451854/1‑149 (MQ=255)
                                                                                  aacaaAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAAt                                                                <  3:508459/149‑1 (MQ=255)
                                                                                  aacaaAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAAt                                                                >  2:226439/1‑149 (MQ=255)
                                                                                   acaaAGTCATTCCAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTtctc                                                                                    >  1:574908/1‑128 (MQ=255)
                                                                                   acaaAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTtctc                                                                                    <  2:574908/128‑1 (MQ=255)
                                                                                   acaaAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATc                                                               >  1:49827/1‑149 (MQ=255)
                                                                                    caaAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCa                                                              <  3:649781/149‑1 (MQ=255)
                                                                                    caaAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCa                                                              <  3:267863/149‑1 (MQ=255)
                                                                                      aaGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGc                                                            >  1:744456/1‑149 (MQ=255)
                                                                                        gTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGt                                                          >  4:438415/1‑149 (MQ=255)
                                                                                        gTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGt                                                          >  2:749286/1‑149 (MQ=255)
                                                                                         tCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTa                                                         >  2:312489/1‑149 (MQ=255)
                                                                                         tCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTa                                                         <  4:45519/149‑1 (MQ=255)
                                                                                           aTTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTAtt                                                       >  1:589321/1‑149 (MQ=255)
                                                                                            ttACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTAttt                                                      >  4:86114/1‑149 (MQ=255)
                                                                                            ttACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTAttt                                                      >  2:614105/1‑149 (MQ=255)
                                                                                            ttACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTAttt                                                      >  1:629614/1‑149 (MQ=255)
                                                                                            ttACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTAttt                                                      >  4:276580/1‑149 (MQ=255)
                                                                                               cAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTtctc                                                   <  1:424832/149‑1 (MQ=255)
                                                                                                 ggCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGa                                                                                                                          >  2:72606/1‑76 (MQ=255)
                                                                                                 ggCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGa                                                                                                                          <  1:72606/76‑1 (MQ=255)
                                                                                                  gCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTaa                                                >  4:149531/1‑149 (MQ=255)
                                                                                                  gCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTaa                                                >  1:248892/1‑149 (MQ=255)
                                                                                                        tAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTc                                          >  2:149065/1‑149 (MQ=255)
                                                                                                            gCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCAcc                                      <  1:659375/149‑1 (MQ=255)
                                                                                                             cGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCg                                     <  1:226439/149‑1 (MQ=255)
                                                                                                                aTGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGAc                                  >  3:208634/1‑149 (MQ=255)
                                                                                                                    tGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGt                              <  1:312489/149‑1 (MQ=255)
                                                                                                                    tGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGt                              <  2:355730/149‑1 (MQ=255)
                                                                                                                    tGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGt                              <  1:204207/149‑1 (MQ=255)
                                                                                                                    tGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCACCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGt                              <  4:523359/149‑1 (MQ=255)
                                                                                                                       tAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGaa                           >  4:238291/1‑149 (MQ=255)
                                                                                                                         aaGCCCATCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCGCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACg                         <  1:596586/149‑1 (MQ=255)
                                                                                                                               gTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGAt                   >  1:508450/1‑149 (MQ=255)
                                                                                                                               gTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGCGAACGGGTGAt                   >  2:774616/1‑149 (MQ=255)
                                                                                                                                 catcaCCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGc                 <  3:718026/149‑1 (MQ=255)
                                                                                                                                    caCCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGa                    >  4:519246/1‑143 (MQ=255)
                                                                                                                                    caCCGACTGGCGATGGCATCGCCTGCGCCAATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGa                    <  3:519246/143‑1 (MQ=255)
                                                                                                                                     aCCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCt             >  1:118577/1‑149 (MQ=255)
                                                                                                                                      ccGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCtc            >  1:254321/1‑149 (MQ=255)
                                                                                                                                       cGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCtca           <  3:405084/149‑1 (MQ=255)
                                                                                                                                         gcaggCGATGCCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCAc                                       <  4:270277/116‑1 (MQ=255)
                                                                                                                                         aCTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCtcatc         <  4:208634/149‑1 (MQ=255)
                                                                                                                                         aCTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCtcatc         <  1:392440/149‑1 (MQ=255)
                                                                                                                                         aCTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCtcatc         <  1:272130/149‑1 (MQ=255)
                                                                                                                                          cTGGCGATGGCATCGCCGGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCtcatca        <  3:734872/149‑1 (MQ=255)
                                                                                                                                            ggCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATa      <  1:774616/149‑1 (MQ=255)
                                                                                                                                            ggCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATa      <  2:248892/149‑1 (MQ=255)
                                                                                                                                             gCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATAc     <  4:71455/149‑1 (MQ=255)
                                                                                                                                               gATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGa   <  2:589321/149‑1 (MQ=255)
                                                                                                                                                aTGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGaa  <  4:159567/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCATCAACAAAGTCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACCGACTGGCGATGGCATCGCCTGCGCCCATAAGCGGTGATGTTCCGGTTTGTTAGCCAGGTCATCCTGTGCTTTCTCGCCCAGCTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAA  >  NC_000913/3857549‑3857841

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