Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,028 T→C G74G (GGT→GGC yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,0280TC92.5% 137.4 / 1.4 53G74G (GGT→GGCyaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (18/31);  total (20/33)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

CATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTC  >  NC_000913/45784‑46265
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
cATTGCAAAAAACAGGAATAACCATGCACCCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCtata                                                                                                                                                                                                                                           <  1:265722/249‑1 (MQ=255)
      aaaaaaCAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGAt                                                                                                                                                                                                                                   >  1:214971/1‑251 (MQ=255)
       aaaaaCAGGAAAACCAGGGCATCCGTGCAGATACTTGGTTGATCTCAAATTCTCCTTGATTCACCGCCGGATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGGTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACGGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCTCATCGCTGTTTGGCTATATTTCCGATa                                                                                                                                                                                                                                  <  2:188395/251‑1 (MQ=255)
          aaCAGGAATACCCATGCACCCGGCCAATAACTTTGGCGATCTCAAATTCCCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAg                                                                                                                                                                                                                               <  2:335891/251‑1 (MQ=255)
                  aaaCCATGCAACCTTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCGCCGCCGTATTTTGCTGTGGGGAAGCGGTGGGCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGTGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACgg                                                                                                                                                                                                                       <  1:479600/250‑1 (MQ=255)
                   aaCCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGtcggtcggg                                                                                                                                                                                                                      >  2:443473/1‑250 (MQ=255)
                      cATGCAACCGTCCAGAAACTTTGACGATCTCACATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGTATGGTTATGTACTGGTACTGGTTGGCGTGGCTCTGGAGCAACTGACTCCGGCTCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGTACGCACACCGGGCTGTTCTTTGACCCAGCGCTGTTTGTCTATATTTCCGGAAAAGACGGAccgccac                                                                                                                                                                                                                   >  2:1836/1‑245 (MQ=255)
                              ccGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGATCGTTGGCACATCGCTGTTTGGCTATATTTCCGCTAAAGTCGGACGGCGACAAATGtt                                                                                                                                                                                                            >  2:124481/1‑251 (MQ=255)
                                  ccAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGGTCCTc                                                                                                                                                                                                        >  2:286203/1‑251 (MQ=255)
                                      aaaCTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCCCAACGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTg                                                                                                                                                                                                    >  2:248896/1‑251 (MQ=255)
                                        aCTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGat                                                                                                                                                                                                  <  1:3822/251‑1 (MQ=255)
                                              gcGATCACAAATTCTCCTTTAGACCCCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGATATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGGTCGTTGGCACATCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATatcatc                                                                                                                                                                                            <  2:108180/250‑1 (MQ=255)
                                              aCGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATatcatc                                                                                                                                                                                            >  1:477456/1‑251 (MQ=255)
                                                 aTCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGcc                                                                                                                                                                                         >  1:218949/1‑251 (MQ=255)
                                                             cctcTATTCGCCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGata                                                                                                                                                                             <  1:370361/251‑1 (MQ=255)
                                                               tcgaTTCACCGCCGCCTTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTAATGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCGGGCGCTAAAACTGGACGCTGAGTGGCTTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCGCATCGCTGTTTGGCTGTATTTTCGATAAAGACGGAGGTCGCAGAATGTTCCTCATTGATATCATCGCCATCTGCGTGgtat                                                                                                                                                                            >  2:259050/4‑250 (MQ=255)
                                                                ctATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGatat                                                                                                                                                                            <  1:135750/249‑1 (MQ=255)
                                                                    tCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTgg                                                                                                                                                                      >  2:225618/1‑251 (MQ=255)
                                                                     cACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCAGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAACCTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGTTATCATCGCCATCGGCGTGATATCGGTGGc                                                                                                                                                                     >  2:467328/1‑251 (MQ=255)
                                                                         gccgcATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGcgacg                                                                                                                                                                 >  1:443181/1‑251 (MQ=255)
                                                                                   gCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCTCATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGAtgtttgtt                                                                                                                                                        <  1:38384/250‑1 (MQ=255)
                                                                                           gAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATccccc                                                                                                                                               >  1:55760/1‑251 (MQ=255)
                                                                                                    gtCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTg                                                                                                                                      <  1:481860/251‑1 (MQ=255)
                                                                                                       cGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGGGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGGCGCTGTCTGGATTGGCTTGTTGGTCGCGGGAACGCTCGTCGGGTTGTTCGTTGGCACATCGCTGTGTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCCTCGTCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTg                                                                                                                                   <  1:198042/251‑1 (MQ=255)
                                                                                                            cTGGATGGTTATGTACTTGTAATGATTGTCTTGGCGCTGGAGCAACTGACGTCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTTTTGTGCGCGGGAACGCTCGCCGGGCTTGTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGTCGATGTTTGTTTCATCCCCCGTCg                                                                                                                                           <  2:94003/238‑1 (MQ=255)
                                                                                                            cTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTAGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCg                                                                                                                                           >  1:94003/1‑238 (MQ=255)
                                                                                                              ggATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGa                                                                                                                               <  1:443473/248‑1 (MQ=255)
                                                                                                              ggATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACAACGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATg                                                                                                                             <  1:248896/250‑1 (MQ=255)
                                                                                                              ggATGGTTATGTACTGGTAATGATTGGCGTGGCCCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATg                                                                                                                             <  1:232692/250‑1 (MQ=255)
                                                                                                                      aTGTTCTGGTAATGATTGGGGGGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTTGAGTGGCTTGTTGGGCGCGGGACCGCTCGCCGGGCTGTTCGTTGGCACATCGCTTTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACtt                                                                                                                    <  2:474904/251‑1 (MQ=255)
                                                                                                                          aCTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGTATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCg                                                                                                                <  2:421628/251‑1 (MQ=255)
                                                                                                                          aCTGGTAATGAGTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATc                                                                                                                 <  2:113755/250‑1 (MQ=255)
                                                                                                                                    aTTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCAt                                                                                                      <  1:124481/251‑1 (MQ=255)
                                                                                                                                         cgtggcgCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGGGATGCGGGTACTTATCGGCATTGTCATCGGTg                                                                                                 >  1:357365/1‑251 (MQ=255)
                                                                                                                                             gcgcTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGc                                                                                                >  1:267522/1‑248 (MQ=255)
                                                                                                                                                     gCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCa                                                                                     <  1:225618/251‑1 (MQ=255)
                                                                                                                                                       aaCTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATc                                                                                   <  1:377164/251‑1 (MQ=255)
                                                                                                                                                                              gACGCTGACTGGATTTGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGTACGGCGCAAAATGTTCCTCATTTATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTc                                                                                                                                            <  2:28417/171‑1 (MQ=255)
                                                                                                                                                                              gACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTc                                                                                                                                            >  1:28417/1‑171 (MQ=255)
                                                                                                                                                                              gACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAAt                                                                        <  1:43444/239‑1 (MQ=255)
                                                                                                                                                                              gACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGGACTTATCGGCATTGTCATCGGGGCAGATTATCCCATCGCCACCTCAAt                                                                        >  2:43444/1‑239 (MQ=255)
                                                                                                                                                                                    gACTGGATTGGCTTGTTGGGCGCGGGAACGCTAGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCAGATAAAGCCGGACGGCGCCAAATGTTCCTCATTGCTATCATCGCCATCGGCGTGCTATCGGTGGCGACGATGTTTGTTTCATCACCCGTCGAACTGTTGGCGCTGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCCATGATCACCGAGTTCTcc                                                        >  1:66930/1‑249 (MQ=255)
                                                                                                                                                                                        ggggtggctggttggtggcgGGACCGCCCGCCGGGTTGGTCGTTGGCACATCGCTTTTTGGCTAAATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGGGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGGTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTAc                                                   <  2:352850/233‑1 (MQ=255)
                                                                                                                                                                                         gATTGGCTTGTTGGGCGCGGGAACGCTCGCCGTGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTAccc                                                 <  1:16195/251‑1 (MQ=255)
                                                                                                                                                                                              gCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTGTGGCTATATTTCCGATAAAGTCGGGCGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCGCAATGATCACCGAGTTCTCCAGTACCCGTCAg                                            <  2:443181/251‑1 (MQ=255)
                                                                                                                                                                                                  gttGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCgg                                         <  2:204270/248‑1 (MQ=255)
                                                                                                                                                                                                           gggAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCCCCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCAt                              <  1:302692/251‑1 (MQ=255)
                                                                                                                                                                                                                cGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGGTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCACCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGc                          <  1:467328/251‑1 (MQ=255)
                                                                                                                                                                                                                    cGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTa                      <  1:471974/251‑1 (MQ=255)
                                                                                                                                                                                                                      ccGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTAtc                    <  1:21323/251‑2 (MQ=255)
                                                                                                                                                                                                                                cGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCAtgt          <  1:204287/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cacaTCGCTGTTTGGCTATATTTCCGATAAAGTCGGCCGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGTGGGTCCTTATCGGCCTTGTCATCGTTTCAGTTTTTCCCATCGTCACCTCAATGATCACCGGGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCCGCTTTAtctccgcccagtggtatg    >  2:228538/1‑234 (MQ=255)
                                                                                                                                                                                                                                        caTCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCTTCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTc  <  2:29988/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
CATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTC  >  NC_000913/45784‑46265

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