Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,227,143 A→G 16.2% T304A (ACC→GCC)  ygjI → putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,227,1430AG16.2% 216.4 / 30.3 92T304A (ACC→GCC) ygjIputative transporter
Reads supporting (aligned to +/- strand):  ref base A (47/30);  new base G (9/6);  total (56/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGGCATGGCGGCGTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAG  >  NC_000913/3227044‑3227242
                                                                                                   |                                                                                                   
cACGGCATGGCGGCGTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGa                                                                                                     >  1:729694/1‑100 (MQ=255)
 aCGGCATGGCGGCGTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGTTTCCCTGCTGATGTGGa                                                                                                     <  2:1900443/99‑1 (MQ=255)
  cGGCATGGCGGCGTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTGCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGAcc                                                                                                   <  2:1260108/100‑1 (MQ=255)
  cGGCATGGCGGCGTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGGcc                                                                                                   <  2:277119/100‑1 (MQ=255)
  cGGCATGGCGGCGTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGAcc                                                                                                   >  1:491445/1‑100 (MQ=255)
  cGGCATGGCGGCGTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGAcc                                                                                                   <  1:886671/100‑1 (MQ=255)
     cATGGCGGCGTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCa                                                                                                >  2:535948/1‑100 (MQ=255)
        ggcggcGTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACg                                                                                             >  1:1320114/1‑100 (MQ=255)
           ggcgTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCg                                                                                          <  1:2290807/100‑1 (MQ=255)
           ggcgTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCACTGCTGATGTGGAACGCAACGCCg                                                                                          >  2:2235746/1‑100 (MQ=255)
             cgTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGcc                                                                                           <  2:354926/97‑1 (MQ=255)
              gTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTg                                                                                       <  1:89439/100‑1 (MQ=255)
                 ttttGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGaaa                                                                                    >  2:494024/1‑100 (MQ=255)
                  tttGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGaaa                                                                                    >  1:197209/1‑99 (MQ=255)
                       tCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGGCCGCAACGCCGGTGAAAATTttc                                                                              <  1:1109611/100‑1 (MQ=255)
                       tCTACCGGAAGCGCTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTttc                                                                              >  1:1581265/1‑100 (MQ=255)
                        cTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTttct                                                                             >  1:1510131/1‑100 (MQ=255)
                        cTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTttct                                                                             <  1:482447/100‑1 (MQ=255)
                          aCCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTttcttc                                                                           >  2:1822556/1‑100 (MQ=255)
                          aCCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTttcttc                                                                           >  1:664649/1‑100 (MQ=255)
                          aCCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTttcttc                                                                           <  2:729691/100‑1 (MQ=255)
                          aCCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTttcttc                                                                           >  1:2120727/1‑100 (MQ=255)
                             ggAAGCGTTGATGAATCGGTTTGTCGGTCTGGTGTCGTTTGCCGCGATGTTCGGTTCGCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTcc                                                                        <  2:1366744/100‑1 (MQ=255)
                                 gCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGGCCGCAACGCCGGTGAAAATTTTCTTCTCCGaaa                                                                    <  1:2287216/100‑1 (MQ=255)
                                  cGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGGCCGCAACGCCGGTGAAAATTTTCTTCTCTGAAAt                                                                   >  1:2115626/1‑100 (MQ=255)
                                    ttGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGGCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATcc                                                                 >  1:1348965/1‑100 (MQ=255)
                                     tgatgaATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATccc                                                                >  1:1948957/1‑100 (MQ=255)
                                     tgatgaATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATccc                                                                >  2:620893/1‑100 (MQ=255)
                                      gatgaATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGGCCGCAACGCCGGTGAAAATTTGCTTCGCCGAAATCCCg                                                               >  1:2205814/1‑100 (MQ=255)
                                      gatgaATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCg                                                               >  2:1516034/1‑100 (MQ=255)
                                       atgaATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCgg                                                              <  1:2040267/100‑1 (MQ=255)
                                        tgaATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGa                                                             >  1:1031222/1‑100 (MQ=255)
                                         gaATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGaa                                                            >  2:657670/1‑100 (MQ=255)
                                         gaATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGaa                                                            <  2:117211/100‑1 (MQ=255)
                                          aaTCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGGCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAg                                                           >  1:241841/1‑100 (MQ=255)
                                           aTCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAgg                                                          <  2:1060577/100‑1 (MQ=255)
                                           aTCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAgg                                                          <  2:1704002/100‑1 (MQ=255)
                                            tCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAg                                                           <  2:1504547/98‑1 (MQ=255)
                                             cGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGGCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCa                                                        <  2:1820621/100‑1 (MQ=255)
                                             cGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCa                                                        <  1:1605728/100‑1 (MQ=255)
                                              gCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCAt                                                       >  2:899250/1‑100 (MQ=255)
                                              gCTTTGTCGGTCTGGTGTCCTTTACCGCGATGCTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCAt                                                       >  2:234921/1‑100 (MQ=255)
                                                  tGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCTCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCttt                                                   >  1:1658760/1‑100 (MQ=255)
                                                   gTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTg                                                  >  2:1842113/1‑100 (MQ=255)
                                                    tCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTgg                                                 <  2:441311/100‑1 (MQ=255)
                                                    tCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTgg                                                 <  2:493546/100‑1 (MQ=255)
                                                    tCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACAGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTgg                                                 <  1:1447530/100‑1 (MQ=255)
                                                     cGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCCGAAGGCATCTTTGGt                                                >  1:839438/1‑100 (MQ=255)
                                                      ggTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGGCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTa                                               >  2:1981506/1‑100 (MQ=255)
                                                      ggTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTa                                               <  1:1453675/100‑1 (MQ=255)
                                                      ggTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTa                                               >  1:907713/1‑100 (MQ=255)
                                                        tCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGGCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAg                                             <  2:2205795/100‑1 (MQ=255)
                                                        tCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAg                                             <  1:1636044/100‑1 (MQ=255)
                                                          tGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCGCCGAAATCCCGGAAGGCTTCTTTGGTAAgta                                           <  1:1585399/100‑3 (MQ=255)
                                                           ggTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGaaa                                          >  2:210825/1‑100 (MQ=255)
                                                            gtgtCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAATGCCGGTGAAAATTTTGTTCTCCGAAATCACGGAGGGCATCTTTGGTAAGaa                                           >  1:467987/1‑98 (MQ=255)
                                                                ccTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGt                                     >  1:1941184/1‑100 (MQ=255)
                                                                ccTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGt                                     >  2:200763/1‑100 (MQ=255)
                                                                 cTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTc                                    >  2:1053518/1‑100 (MQ=255)
                                                                 cTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTc                                    <  1:89416/100‑1 (MQ=255)
                                                                     aCCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAAc                                <  1:1896000/100‑1 (MQ=255)
                                                                      ccGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACt                               >  1:1224421/1‑100 (MQ=255)
                                                                      ccGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACt                               >  2:1112438/1‑100 (MQ=255)
                                                                      ccGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACt                               >  2:1207379/1‑100 (MQ=255)
                                                                      ccGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACt                               >  1:64040/1‑100 (MQ=255)
                                                                        gcgATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGa                             >  2:12334/1‑100 (MQ=255)
                                                                         cgATGTTCGGTTCCCTGCTGATGTGGGCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGaa                            >  1:1173560/1‑100 (MQ=255)
                                                                         cgATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGaa                            >  2:1032229/1‑100 (MQ=255)
                                                                         cgATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGaa                            <  2:1530854/100‑1 (MQ=255)
                                                                          gATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAAc                           <  2:2086173/100‑1 (MQ=255)
                                                                          gATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAAc                           <  1:1195753/100‑1 (MQ=255)
                                                                           aTGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACg                          <  1:1846442/100‑1 (MQ=255)
                                                                            tGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGa                         <  2:771468/100‑1 (MQ=255)
                                                                            tGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGa                         <  1:458758/100‑1 (MQ=255)
                                                                             gTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGaa                        >  1:1914499/1‑100 (MQ=255)
                                                                                cGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAAc                     >  2:1105249/1‑100 (MQ=255)
                                                                                cGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAAc                     >  2:1105126/1‑100 (MQ=255)
                                                                                cGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAAc                     >  2:397964/1‑100 (MQ=255)
                                                                                 ggTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACg                    >  1:933601/1‑100 (MQ=255)
                                                                                  gTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACgg                   >  2:1182962/1‑100 (MQ=255)
                                                                                    tCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCg                 >  2:1866163/1‑100 (MQ=255)
                                                                                     cccTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGt                >  1:2171740/1‑100 (MQ=255)
                                                                                      ccTGCTGATGTGGGCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGtt               >  1:1687042/1‑100 (MQ=255)
                                                                                      ccTGCTGATGTGGGCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGtt               >  1:2190890/1‑100 (MQ=255)
                                                                                      ccTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGtt               >  1:392184/1‑100 (MQ=255)
                                                                                          ctgATGTGGGCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCg            <  1:526241/99‑1 (MQ=255)
                                                                                           tgATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGgc          <  2:335802/100‑1 (MQ=255)
                                                                                            gATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAGCCGTCGAACTGAACGAAAACGGCGTTCCGgcg         <  1:2235765/100‑1 (MQ=255)
                                                                                               gtgGGCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGgcgcgc      >  1:1909168/1‑100 (MQ=255)
                                                                                                 ggACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGgcgcgcgc    >  2:1754241/1‑100 (MQ=255)
                                                                                                 ggACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGgcgcgcgc    >  1:1096572/1‑100 (MQ=255)
                                                                                                 ggACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGgcgcgcgc    >  1:1858867/1‑100 (MQ=255)
                                                                                                 ggACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGgcgcgcgc    >  1:1025739/1‑100 (MQ=255)
                                                                                                 ggACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGgcgcgcgc    >  2:154996/1‑100 (MQ=255)
                                                                                                   aCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAg  >  1:1378131/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CACGGCATGGCGGCGTATTTTGGTCTACCGGAAGCGTTGATGAATCGCTTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAG  >  NC_000913/3227044‑3227242

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: