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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,0280TC98.4% 192.2 / NA 62G74G (GGT→GGC) yaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (29/32);  total (29/33)
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.99e-01

TCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGG‑ACGCTGACTGGAT‑TGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGAC  >  NC_000913/45781‑46273
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
tctcATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGG‑ACGCTGACTGGAT‑TGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCtat                                                                                                                                                                                                                                                    >  1:63072/1‑251 (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)
                                                                                       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)
                                                                                                                    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)
                                                                                                                                                                                                  gcttgttggtggctggaactgcttccGGGCTGGTCGTTGGCAAATTGCTTTTTGGTTATATTTCCGATAAAGTCGGTCGGCGCAAAATGTTCATCATTGATTTCATCGCCATCGGCGTGATATCGGTGTCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAg                                                    <  2:581551/227‑1 (MQ=255)
                                                                                                                                                                                                      gttGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATgtttgttt                                                                                                                                                               >  1:338961/3‑140 (MQ=255)
                                                                                                                                                                                                      gttGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATAGCGGTGGCGACGATgtttgttt                                                                                                                                                               <  2:338961/138‑1 (MQ=255)
                                                                                                                                                                                                            cgcgGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGcttt                                          >  2:162936/1‑251 (MQ=255)
                                                                                                                                                                                                             gcgGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATccc                                                                                                                                                         >  2:208643/1‑139 (MQ=255)
                                                                                                                                                                                                             gcgGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATccc                                                                                                                                                         <  1:208643/139‑1 (MQ=255)
                                                                                                                                                                                                                  aaCGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCGAGTACCCGTCAGCGGGCGTTTTCCATCa                                    <  1:490075/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gTTCGTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATcgccgcca                     <  1:209393/251‑3 (MQ=255)
                                                                                                                                                                                                                                     gTTGGCACATCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCAtgtg                 <  1:356385/251‑1 (MQ=255)
                                                                                                                                                                                                                                           acaTCGCTGCTTGACTATATTTCCGCAAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGACATCGGCGTGATATCGGTGGCCCCGCTGTTTGTTTCATCCCACGTAGTACTGTTGGTTATGGGGGTCCTTTTCGGCATTATAATCGGTGCCGATTATCCCATTGCCCCCTCAATGATCCCCGAGTTCTCCACTACACATATGCAGGCGTTTTCCAACAgtcttttcgtcaccctgtggactgt           >  2:284279/1‑227 (MQ=255)
                                                                                                                                                                                                                                             aTCGCTGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCg         >  1:384570/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  tGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGc                                                                                                                                    <  1:263208/123‑1 (MQ=255)
                                                                                                                                                                                                                                                  tGTTTGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCACGTCGGACTGTTGGTGATGc                                                                                                                                    >  2:263208/1‑123 (MQ=255)
                                                                                                                                                                                                                                                     ttGGCTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGGCGCGAc  <  1:600524/250‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
TCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGG‑ACGCTGACTGGAT‑TGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGAC  >  NC_000913/45781‑46273

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