Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,372,052 A→G 6.6% K48K (AAA→AAG)  ycjN → putative ABC sugar transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,372,0520AG 6.6% 186.0 / 2.3 60K48K (AAA→AAG) ycjNputative ABC sugar transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (32/24);  new base G (2/2);  total (34/26)
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.91e-01

CTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACGTAGCGGTTCGCTGCCGGAAGTGATCGAAACCAGCCATGACTAC  >  NC_000913/1371906‑1372193
                                                                                                                                                  |                                                                                                                                             
ctcATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAatt                                                                                                                                             <  4:416925/149‑1 (MQ=255)
   aTGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAattgat                                                                                                                                          >  4:418273/1‑149 (MQ=255)
        tAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGccc                                                                                                                                     <  3:82111/149‑1 (MQ=255)
         aaaTCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCg                                                                                                                                    >  2:406369/1‑149 (MQ=255)
                       gCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGaaa                                                                                                                      >  4:554110/1‑149 (MQ=255)
                        cTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGaaaa                                                                                                                     <  3:103913/149‑1 (MQ=255)
                        cTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGaaaa                                                                                                                     <  2:273340/149‑1 (MQ=255)
                            tATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCt                                                                                                                 <  4:182144/149‑1 (MQ=255)
                             aTCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTg                                                                                                                >  2:51616/1‑149 (MQ=255)
                              tCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTgg                                                                                                               <  2:151859/149‑1 (MQ=255)
                                   aCTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCa                                                                                                          >  1:298399/1‑149 (MQ=255)
                                   aCTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCa                                                                                                          >  4:263254/1‑149 (MQ=255)
                                   aCTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCa                                                                                                          >  2:228631/1‑149 (MQ=255)
                                    cTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCAc                                                                                                         >  3:9406/1‑149 (MQ=255)
                                       gTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGt                                                                                                      <  1:243512/149‑1 (MQ=255)
                                          tCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAGTTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTaa                                                                                                   >  3:187830/1‑149 (MQ=255)
                                          tCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTaa                                                                                                   <  2:411058/149‑1 (MQ=255)
                                           cATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCG‑TATCATTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGc                                                                                                 <  2:172947/149‑1 (MQ=255)
                                              gcgcCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCaa                                                                                               <  1:179886/149‑1 (MQ=255)
                                                gcCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGt                                                                                             >  1:156563/1‑149 (MQ=255)
                                                 cccTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTg                                                                                            >  3:160045/1‑149 (MQ=255)
                                                 cacTGATTTCAGGCTGTAAAGAAGGAAATAAAACGAATGTATCCATCGACTTTATGCATTCTTCGGAGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAACACCCAGGCATCAACGTTAAGCAAGTg                                                                                            >  4:410363/3‑149 (MQ=255)
                                                  ccTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGc                                                                                           >  3:349899/1‑149 (MQ=255)
                                                   cTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGt                                                                                                                                   <  1:95650/108‑1 (MQ=255)
                                                   cTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGt                                                                                                                                   >  2:95650/1‑108 (MQ=255)
                                                     gATTTCAGGCTGTAAAGAAGAAAATAAAACGAGTGTATCCATCGAATTTATGCATTCTTCTGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAGGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCg                                                                                        >  2:340520/1‑149 (MQ=255)
                                                       tttCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTg                                                                                      <  1:385098/149‑1 (MQ=255)
                                                        ttCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTgg                                                                                     <  4:160045/149‑1 (MQ=255)
                                                          cAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGgaa                                                                                   <  2:340286/149‑1 (MQ=255)
                                                             gCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGgaagaa                                                                                <  3:554110/149‑1 (MQ=255)
                                                                 tAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATg                                                                            <  4:349899/149‑1 (MQ=255)
                                                                    agaagaAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCt                                                                         >  1:77476/1‑149 (MQ=255)
                                                                     gaagaaAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCta                                                                        >  3:145478/1‑149 (MQ=255)
                                                                      aagaaAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCtat                                                                       >  2:119925/1‑149 (MQ=255)
                                                                        gaaAATATAACGACTGTATCTATCCCCTTTATGCGTTCTTCGGTGCAGCAGGAGCGCCAGCCCGTCATCAGCAAATTGATTGCCCGTATCGAAAGCTAAAACCCTGCCATCACCCTAACGCAAGTGCCCGTGGAAGAAGCCGCCTATaa                                                                     >  2:46854/1‑149 (MQ=255)
                                                                        gaaAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAGTTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATaa                                                                     >  2:49257/1‑149 (MQ=255)
                                                                        gaaAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATaa                                                                     >  3:137932/1‑149 (MQ=255)
                                                                             tAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTa                                                                >  4:557470/1‑149 (MQ=255)
                                                                             tAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTa                                                                >  2:236287/1‑149 (MQ=255)
                                                                               aaaCGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTcaa                                                              <  2:375166/149‑1 (MQ=255)
                                                                               aaaCGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTaaa                                                              >  4:147646/1‑149 (MQ=255)
                                                                                    aaTGTATCCATCGAATTTATGGATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCAt                                                         <  4:410084/149‑1 (MQ=255)
                                                                                     aTGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCAtt                                                        <  1:319125/149‑1 (MQ=255)
                                                                                         aTCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTActc                                                    >  2:293211/1‑149 (MQ=255)
                                                                                         aTCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTActc                                                    >  2:207275/1‑149 (MQ=255)
                                                                                             aTCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTc                                                >  2:15815/1‑149 (MQ=255)
                                                                                              tCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCa                                               >  1:255617/1‑149 (MQ=255)
                                                                                              tCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCa                                               >  2:286493/1‑149 (MQ=255)
                                                                                              tCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCa                                               >  4:249411/1‑149 (MQ=255)
                                                                                                gAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACg                                             >  1:82863/1‑149 (MQ=255)
                                                                                                     tATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAGTTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACGTAGCg                                        <  4:187830/149‑1 (MQ=255)
                                                                                                      aTGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAac                                                                    >  1:423498/1‑120 (MQ=255)
                                                                                                      aTGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAac                                                                    <  2:423498/120‑1 (MQ=255)
                                                                                                         cATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACGTAGCGGTTc                                    <  1:400651/149‑1 (MQ=255)
                                                                                                         cATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACGTAGCGGTTc                                    <  2:77476/149‑1 (MQ=255)
                                                                                                                  gTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACGTAGCGGTTCGCTGCCGGa                           >  3:433663/1‑149 (MQ=255)
                                                                                                                        cAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACGTAGCGGTTCGCTGCCGGAAGTGAt                     <  4:298140/149‑1 (MQ=255)
                                                                                                                         aGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACGTAGCGGTTCGCTGCCGGAAGTGATc                    >  4:495666/1‑149 (MQ=255)
                                                                                                                               gcCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACGTAGCGGTTCGCTGCCGGAAGTGATCGAAAcc              <  1:384994/149‑1 (MQ=255)
                                                                                                                                ccAGGCCGTTATCAGTAAGTTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACGTAGCGGTTCGCTGCCGGAAGTGATCGAAACCa             <  1:49257/149‑1 (MQ=255)
                                                                                                                                 cAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACGTAGCGGTTCGCTGCCGGAAGTGATCGAAACCAg            >  3:522453/1‑149 (MQ=255)
                                                                                                                                           tCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACGTAGCGGTTCGCTGCCGGAAGTGATCGAAACCAGCCATGACTAc  <  3:88315/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCGTGGAAGAAGATGCCTATAACACTAAAGTCATTACTCTTTCACGTAGCGGTTCGCTGCCGGAAGTGATCGAAACCAGCCATGACTAC  >  NC_000913/1371906‑1372193

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