Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,227,339 C→T T369I (ACC→ATC)  ygjI → putative transporter YgjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,227,3390CT89.8% 135.4 / ‑5.1 50T369I (ACC→ATC) ygjIputative transporter YgjI
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base T (28/16);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 1.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCTTTGCGT  >  NC_000913/3227092‑3227584
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
tttGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAtcgc                                                                                                                                                                                                                                                    >  1:266889/1‑251 (MQ=255)
 ttGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAtcgc                                                                                                                                                                                                                                                    <  1:330338/250‑1 (MQ=255)
     cGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATcgccgc                                                                                                                                                                                                                                                 >  1:474243/1‑249 (MQ=255)
      ggTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGc                                                                                                                                                                                                                                              >  1:419242/1‑251 (MQ=255)
        tCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGCCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATACCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGt                                                                                                                                                                                                                                            >  2:348282/1‑251 (MQ=255)
           ggTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGCCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTc                                                                                                                                                                                                                                           <  1:431365/249‑1 (MQ=255)
                ccTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCGGCGTGGAGCCATTGCCTGATCGTCATGCCGCTGATGAGGATTCGGATGCTCGGGTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATAGGATCGCCCGATCGTCCATGCt                                                                                                                                                                                                                                     <  2:267996/250‑1 (MQ=255)
                  tttACCGCGATGTTCGGTTCGCTGCTGCTGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACAGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGCTGCTCGGTTACAATACCGTGCAGGATCTGATGCATACTATTTTTAATATGCACGCCGCAGCGTCCATGCTTcc                                                                                                                                                                                                                                  >  2:333272/1‑251 (MQ=255)
                      ccGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTccgccg                                                                                                                                                                                                                              <  2:266889/251‑1 (MQ=255)
                        gcgATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCACTACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGtt                                                                                                                                                                                                                            >  1:358878/1‑251 (MQ=255)
                                   ttCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGCAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGACGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGc                                                                                                                                                                                                                 >  2:521112/1‑251 (MQ=255)
                                         gctgATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATtc                                                                                                                                                                                                                        <  2:409269/238‑1 (MQ=255)
                                         gctgATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATtc                                                                                                                                                                                                                        >  1:409269/1‑238 (MQ=255)
                                            gATGTGGACCGCAACGCCGGTTAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAACCCGTCGAACGGAACGAAACCGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTAGTCATCATGCTGGCTTAc                                                                                                                                                                                                         <  2:285042/250‑1 (MQ=255)
                                               gtgGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTg                                                                                                                                                                                                      <  1:260782/250‑1 (MQ=255)
                                                  gACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAACCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCGTCCGCCGTTATTCATCATGCTGGCTTACCTGAAtt                                                                                                                                                                                                  <  2:487492/251‑1 (MQ=255)
                                                   aCCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAAttt                                                                                                                                                                                                 >  2:19345/1‑251 (MQ=255)
                                                     cGCAACGCCGGTGAAAATATTCTTCACGAAAATCCCGGAAGGCATCTTTGGTAAGAAACCCGTCGAACTGAACGAAAACGGCGTGCCGTCGCGCGCAGGGTGGATCTAGTTCCTAACCGGCATCCCGTGGATGATAATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCTCGGCAGCGTCCATGTTTTCGCCGTTATTCATCATGCTGGTTTTCCGGAATTTAc                                                                                                                                                                                               <  2:367595/251‑1 (MQ=255)
                                                        aaCGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTAcgcg                                                                                                                                                                                            <  1:96085/251‑1 (MQ=255)
                                                        aaCGCCGGTGAAAATTCTCTCCTCCGAAATCCCGGACGGCATCTTTGGTAAGAAAAACGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCAAGTTACTGATCGTCATCCCGCTGCTGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATATGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTAcgcg                                                                                                                                                                                            >  2:360146/1‑251 (MQ=255)
                                                              ggTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAAt                                                                                                                                                                                      >  1:413419/1‑251 (MQ=255)
                                                                        ttcttcTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCt                                                                                                                                                                            >  2:329316/1‑251 (MQ=255)
                                                                          cttctCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGc                                                                                                                                                                          >  2:198547/1‑251 (MQ=255)
                                                                            tctcCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTg                                                                                                                                                                           <  2:138016/248‑1 (MQ=255)
                                                                                cGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCCCCTGCCAcgcg                                                                                                                                                                    <  1:504840/251‑1 (MQ=255)
                                                                                 gAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCAcgcgt                                                                                                                                                                   >  2:14856/1‑250 (MQ=255)
                                                                                                    tgggTAAAAAAACCGCCGAACTGACCGAAATCGGCGTTCCGGGTCGCGCAGCGTGGATCCAGTTCCTGATCGTCCGTCCGTTGATGATTATCCCGAGGCGCGGGTCCAATACCTTGCAGGATCTGATGAATACTATTATTAATATGATCCCCGCAGCGTCCATGATTCCGCCGTTATTCAGCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGTCACGCGATTGCCGTATGGGCTcc                                                                                                                                                   <  2:328888/246‑1 (MQ=255)
                                                                                                          aGAAAACCGTCGAACTGCACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCCGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGGATGGGCTCCCGCCGCAcc                                                                                                                                          >  1:360/1‑251 (MQ=255)
                                                                                                                       aCTGAACAAAAACGGCGTGCGGGCGCGCGCGGCGTGGTCCCAGTTCCTGATCGTCATCCCGCTGATGATTACCCCGATGCTCGGTCCCAATCCCGTGCGGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCTTGTTTCCGCCGTTATTCATCATGCTGGCTTACCGGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCgttg                                                                                                                             <  1:402794/251‑1 (MQ=255)
                                                                                                                         tGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGCTTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCtt                                                                                                                                                                                                           >  2:292005/1‑171 (MQ=255)
                                                                                                                         tGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCtt                                                                                                                                                                                                           <  1:292005/171‑1 (MQ=255)
                                                                                                                                                  cgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCg                                                                                                  >  1:169161/1‑251 (MQ=255)
                                                                                                                                                      gCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGt                                                                                                                                       >  2:191690/1‑210 (MQ=255)
                                                                                                                                                             tCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTc                                                                                       >  1:399965/1‑251 (MQ=255)
                                                                                                                                                                  tcctggaccgcaaccccGTTGATGATTATCCCGATTCTCGTTTCAAATCCCGTGCAGGATCTGATAAATACTATTATTAATATGATCCCCGCAGCGTCCATGCTTCCGCCGTTATTCATAATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCCCCTGCCACGCGGTTTCCTTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCtt                                                                                   <  1:360857/236‑1 (MQ=255)
                                                                                                                                                                      tGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGAc                                                                               >  2:470284/1‑250 (MQ=255)
                                                                                                                                                                       gtccGTCACCCCGCTGATGTTTACCCCGATGCTGGGTCCAAATCCCGGGAGGTTTCTGATGAATATTATTATAAATAGGACCGCCGCAGCGTCCTTGCTTCCGCCTTTATTCATCATGCTGGCTTACCTGAATTGACGCGCCAAATTAGACCCCCTTCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACAt                                                                             <  2:224198/248‑1 (MQ=255)
                                                                                                                                                                                  cGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCAc                                                                                                                                           >  2:45981/1‑178 (MQ=255)
                                                                                                                                                                                  cGCTGATGATTATCCCGATGCTCGGTGCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCAc                                                                                                                                           <  1:45981/178‑1 (MQ=255)
                                                                                                                                                                                    cTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGgcgc                                                                >  2:41481/1‑251 (MQ=255)
                                                                                                                                                                                              ccccGATGCTTGGTTCCAATCCCGGGAAGGATCTGATAAATACTATTATTAATATGACCGCCGCAGCGTCCATGTTTCCGCCGTTATTCATCATGTTGGCTTACCTAAATTTACGCGCCAAATTAGATCCCCTGCCACGCGTTTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTg                                                      <  1:391133/250‑1 (MQ=255)
                                                                                                                                                                                              ccccGATCCGCGGTTCAAATCCCGTGCAGTGTCTTATAAATATTATATTAAATATCATCGCCGCAGCGCCTCTGCTCCCGCTGTTATTCATCTTGCTGGCTTACCTGAATTTACGCGCCAAATAAGATCACCTTCCCCGCGTTTTCCGTTTGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTg                                                      <  2:208629/250‑1 (MQ=255)
                                                                                                                                                                                                        cGGTCCCAATGCCGTGCGGGGTCTGATGAATATTATTATTAATAGGACCCCCGCGGCGTCCTGGCTCCCGCCTTTATTCATCATGCTGGTTTCCCTGAATTTACGCGCCAAATTAGATCCCCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCAttt                                            <  2:474243/251‑1 (MQ=255)
                                                                                                                                                                                                           ttCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCAtttttt                                         >  1:526437/1‑251 (MQ=255)
                                                                                                                                                                                                            tCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCAtttttt                                         >  2:358802/1‑250 (MQ=255)
                                                                                                                                                                                                                aTACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAAc                                    >  2:363249/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gatgaATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGtt                     >  1:509695/1‑250 (MQ=255)
                                                                                                                                                                                                                                    gaATACTATTATTAATATGATCGCCGCCGAGTCCATGCTTCCGCCGTTATTCATCTTGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTCTGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCCTTTTTTATAACGTCGGCGGGATTGTTATCt                 >  2:299455/1‑251 (MQ=255)
                                                                                                                                                                                                                                           attattAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGTACCGGGATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGATTATTCTGACCATCATTTTTTATAACGTCGACGGTATTGTTATCTTGCTTgg          >  2:328080/1‑251 (MQ=255)
                                                                                                                                                                                                                                               taaaTATGACCGCCGCAGCGTCCTGGCTTCCGCCTTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCtt       <  1:358802/248‑1 (MQ=255)
                                                                                                                                                                                                                                                  atatGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCTTTGCg   >  1:303081/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    tatGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATCACGTCGGCGGTATTGTTATCTTCCTTGGCTTTGCgt  >  1:414338/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
TTTGTCGGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCTTTGCGT  >  NC_000913/3227092‑3227584

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