Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 45,978 C→T L58L (CTG→TTG)  yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091345,9780CT96.7% 166.4 / ‑4.0 61L58L (CTG→TTG) yaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (27/31);  total (29/32)
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.45e-01

GTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGG‑ACGCTGACTGGAT‑TGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGT  >  NC_000913/45735‑46206
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    
gTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTACAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGgcgc                                                                                                                                                                                                                               >  1:120487/1‑251 (MQ=255)
         cggctgagtggggTTTGATTCGTAACGCCGAACTGACTCTCATTGCAAAAACCAGGAATAACCAGGCAACCGTCCAGAAACTTGGACGACCTCAAATTCCCCTCTATTCTCCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCt                                                                                                                                                                                                                      <  2:639907/240‑1 (MQ=255)
          gtcTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGG‑ACGCTGACTGGAT‑TGGCTTGCTGGGCGCGGGAACGCTc                                                                                                                                                                                                                     >  1:212958/3‑251 (MQ=255)
          ggCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTc                                                                                                                                                                                                                     >  1:350461/1‑251 (MQ=255)
             tgctGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGTT‑TGGCTTGTTGGGCGCGGGaagtgccgcc                                                                                                                                                                                                                  >  2:209393/4‑243 (MQ=255)
                     gTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAATTTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGG‑ACGCTGACTGGAT‑TGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGtt                                                                                                                                                                                                          <  1:187179/251‑1 (MQ=255)
                            aCGTACCGCCGCATTGACTCTCATTGCAAAAAACAGGAATACCCAGACACCCGTCCAGAAACTTTAACGATCTCAAAGTCTCCTTTATTCACCGCCGCATTTTGCTGTGGGGCAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTATAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCTGGAACGCTCGCCGGGCTTTTCGTTGGc                                                                                                                                                                                                   <  2:277773/251‑1 (MQ=255)
                             cGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGCATAACCATGCAACCGTCCAGACACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGCTGGATATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAACACTGG‑ACGCTGTCTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGcc                                                                                                                                                                                                  >  1:51375/1‑250 (MQ=255)
                                  cgcggcACTGCCTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCATAAACTTTGACGATCTCAAATTCTCCTCTATTCTCCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGTCTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCg                                                                                                                                                                                             <  1:176663/247‑1 (MQ=255)
                                              tctcATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCtat                                                                                                                                                                                 >  1:63072/1‑251 (MQ=255)
                                               ctcATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTggcta                                                                                                                                                                                  <  1:204474/249‑4 (MQ=255)
                                                 cATTGCAAAAAACAGGCATAACCATGCAACCGTCCAGAAACTTTGACGATATCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCACCTGACGCCGGCGCTAAAACTGG‑AAGCTGCCTGGAT‑TGGCATGTTGGGAGCGGGAACGCTCGACAGGCTGTTCGTTGGCACCACGCTGTTTGGCTATAtt                                                                                                                                                                              >  2:658930/1‑251 (MQ=255)
                                                     gCAAAAAAAGGGAATACCCAGGCAACCTCCCAGAAATTTTGTCTATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAATCTGG‑ACGCTGACTGGAT‑TGGCTTGTTTGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGGACATCGCTGTTTGGCTATATTTcc                                                                                                                                                                           <  2:263750/250‑1 (MQ=255)
                                                     acAAAAAACAGGAATCACCATGCAACCGGCCAGAAACTGGGGCGATCCCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAACCTGG‑ACGCTGACTGGAT‑TGGCGTGTTGGGCGCTGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCg                                                                                                                                                                          <  2:337842/250‑1 (MQ=255)
                                                      cAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCCGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGa                                                                                                                                                                         <  1:260723/251‑1 (MQ=255)
                                                       aacaaaCAGGAATAACCATGCAACCGTACAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTACTGTGGGGAAGCGGTGGTCAGTTTCTGGATGGTTATGTACTGGTAATGATTGACGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGTCGCGGGAACGCGCGCCGGGCTGTTCGTTGGCACATCGCTGTGTGGCTATATTTCCGAt                                                                                                                                                                        >  2:265377/4‑251 (MQ=255)
                                                             aagaaaTAACCCTACAACCTACCAGAAATTGTGACGATGTAAAATTCCCCTCTAGTCACCGCCGCATTGTGCTGTGGTGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACGTGTAATGATTGTCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGTAT‑TGGCTTGTTGGGCGCGTGAACGCTCGCCGGGCTGTTCGGTGGCACATCGCTGTTTGGCTATATTTCCGATAAAggc                                                                                                                                                                  <  2:329163/247‑3 (MQ=255)
                                                                    aaCCATGCACCCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCTGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGgc                                                                                                                                                           <  1:647807/251‑1 (MQ=255)
                                                                     aCCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGgcg                                                                                                                                                          <  1:267224/251‑1 (MQ=255)
                                                                      ccATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGgcgc                                                                                                                                                         >  1:328098/1‑251 (MQ=255)
                                                                               cGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGCCGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACCTCGCTGTTTGGCTATATTTCCGTTAAAGTCGGCCGGGTCCAAATGTTc                                                                                                                                                >  2:109080/1‑251 (MQ=255)
                                                                                  ccagtaatgttGAAGATCTCAAATTCTCCTTTTTTCACTGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGGACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGTAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTc                                                                                                                                             <  2:407434/242‑1 (MQ=255)
                                                                                           ttGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATatc                                                                                                                                    >  1:640623/1‑251 (MQ=255)
                                                                                               gtgacacaatgtctcCTTTATTCACCGCCGTTTTTTGATGGGGGGAGGCGGTGGCCCGTTTCTGAATGGTTATGTACTGGTAATGATTGGCGTGGCGTTGGAGCAACTGACGCCGGCGCTAAAACTTG‑TCGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAGAGCTCGTGGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAATGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCg                                                                                                                                <  2:227070/240‑1 (MQ=255)
                                                                                               cGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGCGGTCCGTTTCTGGATGGTTATGTACTGGAAATGATTGACGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATACAGTCGGACGGCGCAAAATGTTCATCATTGATATCATCg                                                                                                                                >  1:664253/1‑251 (MQ=255)
                                                                                                 aTCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGGCGGCGCAAAATGTTCCTCATTGATATCCTCGcc                                                                                                                              >  2:48141/1‑251 (MQ=255)
                                                                                                    tcAAATTATCCTCTATTCACCGCCGAGTTTTGCTGTGGGGACGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGCAATGTTTGGAGTGGCGCTGGCGCAAATGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGATCGGACCGGGCACAATCGCCGGGCTGTTCGTGGGCACATGGCTGTTTGGCTCAATTTCAGTCTAAGTCGGGCGGCGCGAAATGTTCCTCATTGTTatcttctacatc                                                                                                                           >  2:627187/1‑242 (MQ=255)
                                                                                                     cAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATc                                                                                                                           <  1:64839/250‑1 (MQ=255)
                                                                                                      aaaTTCTCCTCTATTCCCCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCgg                                                                                                                         <  1:673696/251‑1 (MQ=255)
                                                                                                                 tATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGatattgg                                                                                                              >  1:603062/1‑248 (MQ=255)
                                                                                                                   ttCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCAGCGGCGTGATATCGGt                                                                                                             <  1:594180/250‑1 (MQ=255)
                                                                                                                          ccgcATTTTGCTGTGGGGGAGCGGTGGACCGATTCTGGATGGTTATGAACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTAG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGGACGCTCGCCGGGCTGTTCGTTGACCCATCCATGTTTGACTATATTTCCGATAAAGTCGCACGGCGCACAATGTGCCTTATTGATAACATCGCCATCGGCGTGATATCGGTggagctgg                                                                                                     >  2:428597/1‑245 (MQ=255)
                                                                                                                            gcATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TTGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGAtg                                                                                                   <  1:635243/251‑1 (MQ=255)
                                                                                                                            gcATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCAGCGCCATCGGCGTGATATCGGTGGCGACGAtg                                                                                                   <  1:536161/251‑1 (MQ=255)
                                                                                                                               ttttGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTGCTGGTAATGATTGGCGTGGCGCTGCAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGcgcg                                                                                                                                                                                                                              <  1:660018/125‑1 (MQ=255)
                                                                                                                                     tgtgGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGGGATATCGGTGGCGACGATGTTTGTTTc                                                                                           >  1:607070/1‑250 (MQ=255)
                                                                                                                                            aaGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCg                                                                                   <  1:48141/251‑1 (MQ=255)
                                                                                                                                                 gtgtgccttttCTGGATGGTGATGTGCTGGGAATGATTGTAGTGGCGCTGTAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGTCTGGAT‑TGGCTTGTTGGGCGCGGGAACTCTCGCCTGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGTACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATGGGTGGCGACGATGTTTGTTTCATCCCCCGTCg                                                                                <  2:470229/241‑1 (MQ=255)
                                                                                                                                                   ggtCCGTTTCTGGATGGTTATGGACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGCTATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGacct                                                                            >  1:24812/1‑248 (MQ=255)
                                                                                                                                                    gtCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTg                                                                           <  1:486230/251‑1 (MQ=255)
                                                                                                                                                          ttCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGGTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCTCCATCGGCGTGATATCGGTGGCGGCGATGTTTGTTTCATCCCCCGTCGGACtttttgtt                                                                     >  2:178434/1‑244 (MQ=255)
                                                                                                                                                           tCTGGATGGTTATGTACTGGTAATGATGGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGCT‑TGGCTTGTTGGGCGCGGGAACGTTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCTGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCTCGTCGTACTGTTgggg                                                                     <  2:640623/250‑3 (MQ=255)
                                                                                                                                                             tggtaggTTATGTTCTGGTAACTATGTGCGTGGCGCTGGAGCAACTTACGCCGGTGCT‑AAATTGGGACTCGGGCTTGATGTGG‑TTGTTGGGCGCGGGAACGCTAGCCGGGCTGTTCGTTGGCACTTCGCTGTTTTGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATGGATATCAGCGCCATCGTCGTGATAGCGGTGGCGACGATGGGTGTTTCGGCCCCCGTCGTACTGGTGGTGAt                                                                   <  2:278452/245‑1 (MQ=255)
                                                                                                                                                                aTGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCt                                                                                                                                                                                           <  1:514089/127‑1 (MQ=255)
                                                                                                                                                                aTGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGAAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGCT‑TGGCTTGTTGGGCGCGGGGACGCGCGCCGGGCTGTTCGTTGGCACATCGCt                                                                                                                                                                                           >  2:514089/1‑127 (MQ=255)
                                                                                                                                                                  ggTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCggg                                                              <  1:455094/250‑1 (MQ=255)
                                                                                                                                                                       tGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACtt                                                         >  1:415308/1‑250 (MQ=255)
                                                                                                                                                                                aaTGATTGGAGTGGCGATGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGCT‑TGGCTTGTTGGGCGCGGGAACGATCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGACGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGTCGATGTTTGTTTCATCCCCCGTCGGACTGTTGGTGGTGCGGGTACTTATCGGCCTTg                                               >  2:84966/1‑251 (MQ=255)
                                                                                                                                                                                   gATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCa                                            <  1:530266/251‑1 (MQ=255)
                                                                                                                                                                                          gtgacgcTGGCGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGCT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTACGTTAAAGTCGGTCGGCGCAAAATGATCCTCCTTGATATCATCGCCATCGGCGTGAAATAGGTGGCGACGATGTTTGTTTTATCCCCCCTCGTACTGTTGTGGATCCGGGTACtttttagcattgtgctaggtgc                                     >  2:77043/5‑231 (MQ=255)
                                                                                                                                                                                          ggggcgCGGGATCAACTGACGCCGGCGGTAAAGCGGG‑ACGCTGACGGGAT‑TGGCTTGTTGGTCGCGGGATCGCTCGCTGGGCTGTGCGTTGGCACATCGCTTTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGTCGATGTTGGTTTCATCCCCCGGCGAACTGTGGGTGATGCGGGTACTTAACGGCATTTTCATCGGTGc                                     <  2:239256/249‑1 (MQ=255)
                                                                                                                                                                                                        aCTGACGCGGTCGCTAAAACGGG‑ACGCTGACTGGAT‑TTGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTGCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGGTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCg                       <  2:328098/251‑1 (MQ=255)
                                                                                                                                                                                                            aCGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCAc                   >  1:98193/1‑251 (MQ=255)
                                                                                                                                                                                                              gCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGTTATCATCGCCATCGGCGTGATATCGGGGGCGGCGATGTTTGTTTTCTCCCCCGTCGAACTGTTGGTGCTGCGGGGACTTATCTGCCTAGTCATCGGTACAGCTTATCCCATcccccccc                 >  2:356385/1‑244 (MQ=255)
                                                                                                                                                                                                              gCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCt                 >  1:478774/1‑251 (MQ=255)
                                                                                                                                                                                                                       aaaaCTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCa        <  1:84966/250‑1 (MQ=255)
                                                                                                                                                                                                                        tagcTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAACTGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGGTGTTTGTTTCATCCCCCGTCGAACTGTTGTGGATGCGGGTACTTATCGGCGTTGTCATCGGTGCAATTTATCCCATCGCACCCTCAATGAtaac       >  2:504944/4‑248 (MQ=255)
                                                                                                                                                                                                                          aCTGG‑GCGCTGCCTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTGCGTTGGCACATCTCTGTTTGGCTATATTTCCGATAAAGGCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGAGTATCCCATCGCCACCTCAATGATCACCg     <  2:350461/251‑1 (MQ=255)
                                                                                                                                                                                                                             gg‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATCTCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGGGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCTCCTCatagctcacagtgt  >  2:178085/1‑238 (MQ=255)
                                                                                                                                                                                                                                     gacgggAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGGTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTAc                                                           <  2:87788/183‑1 (MQ=255)
                                                                                                                                                                                                                                     gACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTAc                                                           >  1:87788/1‑187 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    
GTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGG‑ACGCTGACTGGAT‑TGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGT  >  NC_000913/45735‑46206

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