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

GTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGT  >  NC_000913/45735‑46206
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    
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)
                                                                                                                                                                                                                       aaaaCTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCa        <  1:84966/250‑1 (MQ=255)
                                                                                                                                                                                                                        tagcTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAACTGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGGTGTTTGTTTCATCCCCCGTCGAACTGTTGTGGATGCGGGTACTTATCGGCGTTGTCATCGGTGCAATTTATCCCATCGCACCCTCAATGAtaac       >  2:504944/4‑248 (MQ=255)
                                                                                                                                                                                                                          aCTGGGCGCTGCCTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTGCGTTGGCACATCTCTGTTTGGCTATATTTCCGATAAAGGCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGAGTATCCCATCGCCACCTCAATGATCACCg     <  2:350461/251‑1 (MQ=255)
                                                                                                                                                                                                                             ggACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATCTCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGGGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCTCCTCatagctcacagtgt  >  2:178085/1‑238 (MQ=255)
                                                                                                                                                                                                                                     gacgggATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGGTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTAc                                                           <  2:87788/183‑1 (MQ=255)
                                                                                                                                                                                                                                     gACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTAc                                                           >  1:87788/1‑187 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    
GTTCCGTTACGGCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGT  >  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: