Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,439,389 G→A 6.3% A244A (GCC→GCT)  yfcJ ← putative arabinose efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,439,3890GA 6.3% 198.9 / 3.7 63A244A (GCC→GCT) yfcJputative arabinose efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base A (2/2);  total (31/32)
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.19e-02

ACAAGCCCACCGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGCCCGACAACGCTCCACAGCGATGGACGTTCTCCCGCCAGGGCCGGTACTTTGCGCACTGTGCCG  >  NC_000913/2439241‑2439536
                                                                                                                                                    |                                                                                                                                                   
acaaGCCCACCGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTg                                                                                                                                                     >  2:534949/1‑149 (MQ=255)
      ccACCGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAg                                                                                                                                               >  3:310324/1‑149 (MQ=255)
       cACCGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAACCGGTCCGGCAGCCAGCCAAACATGACGCGGATCACGACAAATGCGCCGGCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTTCTGGTAAAGt                                                                                                                                              <  2:420753/149‑1 (MQ=255)
       cACCGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGGCGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTg                                                                                                                                                     <  1:190422/142‑1 (MQ=255)
       cACCGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGGCGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTg                                                                                                                                                     >  2:190422/1‑142 (MQ=255)
       cACCGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGt                                                                                                                                              >  1:594377/1‑149 (MQ=255)
        aCCGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCCAGCAAAgga                                                                                                                                             >  3:526604/1‑147 (MQ=255)
        aCCGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCCAGCAAAgga                                                                                                                                             >  1:425066/1‑147 (MQ=255)
          cGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTaga                                                                                                                                           <  3:261991/149‑1 (MQ=255)
          cGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCCAGCAAAggatg                                                                                                                                           <  2:425066/149‑5 (MQ=255)
          cGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCCAGCAAAggatg                                                                                                                                           <  4:526604/149‑5 (MQ=255)
           gTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTagag                                                                                                                                          >  3:321933/1‑149 (MQ=255)
              tCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGcga                                                                                                                                       >  1:708825/1‑149 (MQ=255)
                tACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGcgaaa                                                                                                                                     >  1:507407/1‑149 (MQ=255)
                 aCAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGcgaaac                                                                                                                                    <  2:657403/149‑1 (MQ=255)
                 aCAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGcgaaac                                                                                                                                    <  2:696416/149‑1 (MQ=255)
                    aGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGcgaaacgaa                                                                                                                                 >  1:557820/1‑149 (MQ=255)
                    aGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGcgaaacgaa                                                                                                                                 >  4:219226/1‑149 (MQ=255)
                     gCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGTCCATCCTTTGCTGGCAAAGTAGAGcgaaacgaaa                                                                                                                                >  3:455047/1‑149 (MQ=255)
                     gCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGcgaaacgaaa                                                                                                                                <  1:575890/149‑1 (MQ=255)
                     gCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGcgaaacgaaa                                                                                                                                <  1:620918/149‑1 (MQ=255)
                     gCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGcgaaacgaaa                                                                                                                                <  4:254082/149‑1 (MQ=255)
                     gCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGcgaaacgaaa                                                                                                                                >  3:578129/1‑149 (MQ=255)
                         agagaCAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACAACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTcc                                                                                                                            >  1:109338/1‑149 (MQ=255)
                           agaCAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAg                                                                                                                                               <  4:188512/128‑1 (MQ=255)
                           agaCAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAg                                                                                                                                               >  3:188512/1‑128 (MQ=255)
                               aaTCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTcc                                                                                                                            <  3:399113/143‑1 (MQ=255)
                               aaTCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTcc                                                                                                                            >  4:399113/1‑143 (MQ=255)
                                aTCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAAc                                                                                                                     <  3:667788/149‑1 (MQ=255)
                                   gCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGc                                                                                                                  >  3:466940/1‑149 (MQ=255)
                                     cACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCaa                                                                                                                >  1:578230/1‑149 (MQ=255)
                                       cTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCaaaa                                                                                                              <  4:622660/149‑1 (MQ=255)
                                        tttCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAAc                                                                                                             >  1:564409/1‑149 (MQ=255)
                                         ttCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAATcc                                                                                                            <  4:505695/149‑1 (MQ=255)
                                           cACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCaa                                                                                                          <  3:45624/149‑1 (MQ=255)
                                              gccgccAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGc                                                                                                       <  1:255108/149‑1 (MQ=255)
                                                    aaaaCGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTa                                                                                                 >  1:164636/1‑149 (MQ=255)
                                                    aaaaCGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTa                                                                                                 <  1:63649/149‑1 (MQ=255)
                                                      aaCGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGt                                                                                               >  3:10435/1‑149 (MQ=255)
                                                       aCGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTg                                                                                              <  1:201834/149‑1 (MQ=255)
                                                         ggTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCGTCGCCCATCCTTTGCCGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGcc                                                                                            <  4:298461/149‑1 (MQ=255)
                                                         ggTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGcc                                                                                            <  1:440826/149‑1 (MQ=255)
                                                          gTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCa                                                                                           <  3:240108/149‑1 (MQ=255)
                                                             cGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGAc                                                                                        >  3:171782/1‑149 (MQ=255)
                                                                     aGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCt                                                                                >  3:93394/1‑149 (MQ=255)
                                                                      gCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTg                                                                               <  2:414454/149‑1 (MQ=255)
                                                                         aaaCATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGtt                                                                            <  3:327898/149‑1 (MQ=255)
                                                                         aaaCATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGtt                                                                            >  1:541956/1‑149 (MQ=255)
                                                                           aCATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGAc                                                                                        >  2:629777/1‑135 (MQ=255)
                                                                           aCATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGAc                                                                                        <  1:629777/135‑1 (MQ=255)
                                                                             aTGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCa                                                                        <  2:378996/149‑1 (MQ=255)
                                                                             aTGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCa                                                                        <  4:131893/149‑1 (MQ=255)
                                                                                aCGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGAt                                                                     <  4:578129/149‑1 (MQ=255)
                                                                                 cgcgCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTcc                                                                         <  3:375166/144‑1 (MQ=255)
                                                                                 cgcgCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTcc                                                                         >  4:375166/1‑144 (MQ=255)
                                                                                      aTCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGccc                                                               <  3:640253/149‑1 (MQ=255)
                                                                                                            aCGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGCCCGACAACGCTCCACAGCGATGGa                                         >  2:567840/1‑149 (MQ=255)
                                                                                                                          aaGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGCCCGACAACGCTCCACAGCGATGGACGTTCTCCCGCCAg                           >  2:414840/1‑149 (MQ=255)
                                                                                                                                      gCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGCCCGACAACGCTCCACAGCGATGGACGTTCTCCCGCCAGGGCCGGTACttt               >  2:403720/1‑149 (MQ=255)
                                                                                                                                      gCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGCCCGACAACGCTCCACAGCGATGGACGTTCTCCCGCCAGGGCCGGTACttt               <  4:455047/149‑1 (MQ=255)
                                                                                                                                        ccATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGCCCGACAACGCTCCACAGCGATGGACGTTCTCCCGCCAGGGCCGGTACTTTgc             <  1:403720/149‑1 (MQ=255)
                                                                                                                                        ccATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGCCCGACAACGCTCCACAGCGATGGACGTTCTCCCGCCAGGGCCGGTACTTTgc             >  3:521380/1‑149 (MQ=255)
                                                                                                                                        ccATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGCCCGACAACGCTCCACAGCGATGGACGTTCTCCCGCCAGGGCCGGTACTTTgc             <  2:564409/149‑1 (MQ=255)
                                                                                                                                             cTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGCCCGACAACGCTCCACAGCGATGGACGTTCTCCCGCCAGGGCCGGTACTTTGCGCACt        >  1:706269/1‑149 (MQ=255)
                                                                                                                                              tttGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGCCCGACAACGCTCCACAGCGATGGACGTTCTCCCGCCAGGGCCGGTACTTTGCGCACtg       >  4:107081/1‑149 (MQ=255)
                                                                                                                                                   tGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGCCCGACAACGCTCCACAGCGATGGACGTTCTCCCGCCAGGGCCGGTACTTTGCGCACTGTGCCg  >  3:533160/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACAAGCCCACCGTTTCTACAAGCAGAGAGACAATCGCCACTTTCACGCCGCCAAAACGGTCCGGCATCCAGCCAAACATGACGCGCATCACGACAAATGCGCCGCCAAACGCGGTAAGAGTAAAGCCCGCCATCGCCCATCCTTTGCTGGCAAAGTAGAGCGAAACGAAAGTCCCGATAACAGCAAAACCAACGCCTTGTAGTGCCAGACCTAACCCTGGTTTCCAGATAAGCCCGACAACGCTCCACAGCGATGGACGTTCTCCCGCCAGGGCCGGTACTTTGCGCACTGTGCCG  >  NC_000913/2439241‑2439536

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: ‑