Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 45,950 G→A L48L (CTG→CTA yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091345,9500GA96.5% 167.6 / ‑2.5 58L48L (CTG→CTAyaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (26/29);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

GGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCT  >  NC_000913/45718‑46191
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
ggccccTTTGTTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACATAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAAACGTCCAGACACTTTGACGATCTCCAATTCTCCTCTATTCCACGCCGCATTTTGCAGTGGGGACGCGATGGTCCCTGTCTGGATGGTTATGTGCTGGTAAAGGTTGACGGTGCGCTGTATCACATGACGGACGCGCGAAAACTGGACGCTGACTgg                                                                                                                                                                                                                                  >  2:539730/5‑251 (MQ=255)
      gttgtgggtgctctCCGTTACGGTTCAGGTTGGTTTGATACGTAACGCTCCACTGACTCTCATTGCAAAAAACAGGAATAGCCATGCATCCGCCCAGAACCTTTTACGTTCTCAAATTCTCCCCTATTCACCGCCGCATTTTGCTGTTGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGTCTGGATTg                                                                                                                                                                                                                             <  2:467414/236‑1 (MQ=255)
          gggtggatgcctgTTACGGCTGATGTTAGTTTGAACCGTAACGCCGCTCTGACCTTCATGGCAAAAAACAGGAATAACCCTGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACGGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTgt                                                                                                                                                                                                                       <  2:469500/239‑2 (MQ=255)
                 gTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTACAACTGGACGCTGACTGGATTGGCTTGTTGGgcgc                                                                                                                                                                                                                >  1:120487/1‑251 (MQ=255)
                          cggctgagtggggTTTGATTCGTAACGCCGAACTGACTCTCATTGCAAAAACCAGGAATAACCAGGCAACCGTCCAGAAACTTGGACGACCTCAAATTCCCCTCTATTCTCCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCt                                                                                                                                                                                                       <  2:639907/240‑1 (MQ=255)
                           gtcTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTc                                                                                                                                                                                                      >  1:212958/3‑251 (MQ=255)
                           ggCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTc                                                                                                                                                                                                      >  1:350461/1‑251 (MQ=255)
                              tgctGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGTTTGGCTTGTTGGGCGCGGGaagtgccgcc                                                                                                                                                                                                   >  2:209393/4‑243 (MQ=255)
                                      gTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAATTTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGtt                                                                                                                                                                                           <  1:187179/251‑1 (MQ=255)
                                             aCGTACCGCCGCATTGACTCTCATTGCAAAAAACAGGAATACCCAGACACCCGTCCAGAAACTTTAACGATCTCAAAGTCTCCTTTATTCACCGCCGCATTTTGCTGTGGGGCAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTATAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCTGGAACGCTCGCCGGGCTTTTCGTTGGc                                                                                                                                                                                    <  2:277773/251‑1 (MQ=255)
                                              cGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGCATAACCATGCAACCGTCCAGACACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGCTGGATATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAACACTGGACGCTGTCTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGcc                                                                                                                                                                                   >  1:51375/1‑250 (MQ=255)
                                                   cgcggcACTGCCTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCATAAACTTTGACGATCTCAAATTCTCCTCTATTCTCCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGTCTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCg                                                                                                                                                                              <  1:176663/247‑1 (MQ=255)
                                                               tctcATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCtat                                                                                                                                                                  >  1:63072/1‑251 (MQ=255)
                                                                ctcATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTggcta                                                                                                                                                                   <  1:204474/249‑4 (MQ=255)
                                                                  cATTGCAAAAAACAGGCATAACCATGCAACCGTCCAGAAACTTTGACGATATCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCACCTGACGCCGGCGCTAAAACTGGAAGCTGCCTGGATTGGCATGTTGGGAGCGGGAACGCTCGACAGGCTGTTCGTTGGCACCACGCTGTTTGGCTATAtt                                                                                                                                                               >  2:658930/1‑251 (MQ=255)
                                                                      gCAAAAAAAGGGAATACCCAGGCAACCTCCCAGAAATTTTGTCTATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAATCTGGACGCTGACTGGATTGGCTTGTTTGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGGACATCGCTGTTTGGCTATATTTcc                                                                                                                                                            <  2:263750/250‑1 (MQ=255)
                                                                      acAAAAAACAGGAATCACCATGCAACCGGCCAGAAACTGGGGCGATCCCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAACCTGGACGCTGACTGGATTGGCGTGTTGGGCGCTGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCg                                                                                                                                                           <  2:337842/250‑1 (MQ=255)
                                                                       cAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCCGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGa                                                                                                                                                          <  1:260723/251‑1 (MQ=255)
                                                                        aacaaaCAGGAATAACCATGCAACCGTACAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTACTGTGGGGAAGCGGTGGTCAGTTTCTGGATGGTTATGTACTGGTAATGATTGACGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGTCGCGGGAACGCGCGCCGGGCTGTTCGTTGGCACATCGCTGTGTGGCTATATTTCCGAt                                                                                                                                                         >  2:265377/4‑251 (MQ=255)
                                                                              aagaaaTAACCCTACAACCTACCAGAAATTGTGACGATGTAAAATTCCCCTCTAGTCACCGCCGCATTGTGCTGTGGTGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACGTGTAATGATTGTCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGTATTGGCTTGTTGGGCGCGTGAACGCTCGCCGGGCTGTTCGGTGGCACATCGCTGTTTGGCTATATTTCCGATAAAggc                                                                                                                                                   <  2:329163/247‑3 (MQ=255)
                                                                                     aaCCATGCACCCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCTGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGgc                                                                                                                                            <  1:647807/251‑1 (MQ=255)
                                                                                      aCCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGgcg                                                                                                                                           <  1:267224/251‑1 (MQ=255)
                                                                                       ccATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGgcgc                                                                                                                                          >  1:328098/1‑251 (MQ=255)
                                                                                                cGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGCCGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACCTCGCTGTTTGGCTATATTTCCGTTAAAGTCGGCCGGGTCCAAATGTTc                                                                                                                                 >  2:109080/1‑251 (MQ=255)
                                                                                                   ccagtaatgttGAAGATCTCAAATTCTCCTTTTTTCACTGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGGACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGTATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTc                                                                                                                              <  2:407434/242‑1 (MQ=255)
                                                                                                            ttGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATatc                                                                                                                     >  1:640623/1‑251 (MQ=255)
                                                                                                                gtgacacaatgtctcCTTTATTCACCGCCGTTTTTTGATGGGGGGAGGCGGTGGCCCGTTTCTGAATGGTTATGTACTGGTAATGATTGGCGTGGCGTTGGAGCAACTGACGCCGGCGCTAAAACTTGTCGCTGACTGGATTGGCTTGTTGGGCGCGGGAGAGCTCGTGGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAATGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCg                                                                                                                 <  2:227070/240‑1 (MQ=255)
                                                                                                                cGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGCGGTCCGTTTCTGGATGGTTATGTACTGGAAATGATTGACGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATACAGTCGGACGGCGCAAAATGTTCATCATTGATATCATCg                                                                                                                 >  1:664253/1‑251 (MQ=255)
                                                                                                                  aTCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGGCGGCGCAAAATGTTCCTCATTGATATCCTCGcc                                                                                                               >  2:48141/1‑251 (MQ=255)
                                                                                                                     tcAAATTATCCTCTATTCACCGCCGAGTTTTGCTGTGGGGACGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGCAATGTTTGGAGTGGCGCTGGCGCAAATGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGATCGGACCGGGCACAATCGCCGGGCTGTTCGTGGGCACATGGCTGTTTGGCTCAATTTCAGTCTAAGTCGGGCGGCGCGAAATGTTCCTCATTGTTatcttctacatc                                                                                                            >  2:627187/1‑242 (MQ=255)
                                                                                                                      cAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATc                                                                                                            <  1:64839/250‑1 (MQ=255)
                                                                                                                       aaaTTCTCCTCTATTCCCCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCgg                                                                                                          <  1:673696/251‑1 (MQ=255)
                                                                                                                                  tATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGatattgg                                                                                               >  1:603062/1‑248 (MQ=255)
                                                                                                                                    ttCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCAGCGGCGTGATATCGGt                                                                                              <  1:594180/250‑1 (MQ=255)
                                                                                                                                           ccgcATTTTGCTGTGGGGGAGCGGTGGACCGATTCTGGATGGTTATGAACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTAGACGCTGACTGGATTGGCTTGTTGGGCGCGGGGACGCTCGCCGGGCTGTTCGTTGACCCATCCATGTTTGACTATATTTCCGATAAAGTCGCACGGCGCACAATGTGCCTTATTGATAACATCGCCATCGGCGTGATATCGGTggagctgg                                                                                      >  2:428597/1‑245 (MQ=255)
                                                                                                                                             gcATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTTGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGAtg                                                                                    <  1:635243/251‑1 (MQ=255)
                                                                                                                                             gcATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCAGCGCCATCGGCGTGATATCGGTGGCGACGAtg                                                                                    <  1:536161/251‑1 (MQ=255)
                                                                                                                                                ttttGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTGCTGGTAATGATTGGCGTGGCGCTGCAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGcgcg                                                                                                                                                                                                               <  1:660018/125‑1 (MQ=255)
                                                                                                                                                      tgtgGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGGGATATCGGTGGCGACGATGTTTGTTTc                                                                            >  1:607070/1‑250 (MQ=255)
                                                                                                                                                             aaGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCg                                                                    <  1:48141/251‑1 (MQ=255)
                                                                                                                                                                  gtgtgccttttCTGGATGGTGATGTGCTGGGAATGATTGTAGTGGCGCTGTAGCAACTGACGCCGGCGCTAAAACTGGACGCTGTCTGGATTGGCTTGTTGGGCGCGGGAACTCTCGCCTGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGTACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATGGGTGGCGACGATGTTTGTTTCATCCCCCGTCg                                                                 <  2:470229/241‑1 (MQ=255)
                                                                                                                                                                    ggtCCGTTTCTGGATGGTTATGGACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGCTATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGacct                                                             >  1:24812/1‑248 (MQ=255)
                                                                                                                                                                     gtCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTg                                                            <  1:486230/251‑1 (MQ=255)
                                                                                                                                                                           ttCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGGTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCTCCATCGGCGTGATATCGGTGGCGGCGATGTTTGTTTCATCCCCCGTCGGACtttttgtt                                                      >  2:178434/1‑244 (MQ=255)
                                                                                                                                                                            tCTGGATGGTTATGTACTGGTAATGATGGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGCTTGGCTTGTTGGGCGCGGGAACGTTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCTGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCTCGTCGTACTGTTgggg                                                      <  2:640623/250‑3 (MQ=255)
                                                                                                                                                                              tggtaggTTATGTTCTGGTAACTATGTGCGTGGCGCTGGAGCAACTTACGCCGGTGCTAAATTGGGACTCGGGCTTGATGTGGTTGTTGGGCGCGGGAACGCTAGCCGGGCTGTTCGTTGGCACTTCGCTGTTTTGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATGGATATCAGCGCCATCGTCGTGATAGCGGTGGCGACGATGGGTGTTTCGGCCCCCGTCGTACTGGTGGTGAt                                                    <  2:278452/245‑1 (MQ=255)
                                                                                                                                                                                 aTGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCt                                                                                                                                                                            <  1:514089/127‑1 (MQ=255)
                                                                                                                                                                                 aTGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGAAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGCTTGGCTTGTTGGGCGCGGGGACGCGCGCCGGGCTGTTCGTTGGCACATCGCt                                                                                                                                                                            >  2:514089/1‑127 (MQ=255)
                                                                                                                                                                                   ggTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCggg                                               <  1:455094/250‑1 (MQ=255)
                                                                                                                                                                                        tGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACtt                                          >  1:415308/1‑250 (MQ=255)
                                                                                                                                                                                                 aaTGATTGGAGTGGCGATGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGCTTGGCTTGTTGGGCGCGGGAACGATCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGACGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGTCGATGTTTGTTTCATCCCCCGTCGGACTGTTGGTGGTGCGGGTACTTATCGGCCTTg                                >  2:84966/1‑251 (MQ=255)
                                                                                                                                                                                                    gATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCa                             <  1:530266/251‑1 (MQ=255)
                                                                                                                                                                                                           gtgacgcTGGCGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGCTTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTACGTTAAAGTCGGTCGGCGCAAAATGATCCTCCTTGATATCATCGCCATCGGCGTGAAATAGGTGGCGACGATGTTTGTTTTATCCCCCCTCGTACTGTTGTGGATCCGGGTACtttttagcattgtgctaggtgc                      >  2:77043/5‑231 (MQ=255)
                                                                                                                                                                                                           ggggcgCGGGATCAACTGACGCCGGCGGTAAAGCGGGACGCTGACGGGATTGGCTTGTTGGTCGCGGGATCGCTCGCTGGGCTGTGCGTTGGCACATCGCTTTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGTCGATGTTGGTTTCATCCCCCGGCGAACTGTGGGTGATGCGGGTACTTAACGGCATTTTCATCGGTGc                      <  2:239256/249‑1 (MQ=255)
                                                                                                                                                                                                                         aCTGACGCGGTCGCTAAAACGGGACGCTGACTGGATTTGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTGCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGGTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCg        <  2:328098/251‑1 (MQ=255)
                                                                                                                                                                                                                             aCGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCAc    >  1:98193/1‑251 (MQ=255)
                                                                                                                                                                                                                               gCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGTTATCATCGCCATCGGCGTGATATCGGGGGCGGCGATGTTTGTTTTCTCCCCCGTCGAACTGTTGGTGCTGCGGGGACTTATCTGCCTAGTCATCGGTACAGCTTATCCCATcccccccc  >  2:356385/1‑244 (MQ=255)
                                                                                                                                                                                                                               gCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCt  >  1:478774/1‑251 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
GGCACCTTTGGTGTGGAGTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCT  >  NC_000913/45718‑46191

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