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→GGCyaaUputative 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

TCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGAC  >  NC_000913/45781‑46273
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
tctcATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCtat                                                                                                                                                                                                                                                    >  1:63072/1‑251 (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)
                                                                                       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)
                                                                                                                    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)
                                                                                                                                                                                                  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)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
TCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGCGAC  >  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: