Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,943 T→C S379S (AGT→AGC yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,9430TC98.2% 174.7 / NA 55S379S (AGT→AGCyaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (24/30);  minor base G (0/1);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCG  >  NC_000913/46709‑47188
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
gtgGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCAt                                                                                                                                                                                                                                         <  1:341949/249‑1 (MQ=255)
   gtgaaaCACTGCCGTACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACTCTGTCGCTGGCGGTTTTGGGTCTGATCCCGGATATGGGGTTCTGGCTGGTAGTGATGGCTTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCGCCATTGttt                                                                                                                                                                                                                                    <  1:23149/248‑1 (MQ=255)
     tAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCGGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTc                                                                                                                                                                                                                                   <  1:618116/250‑1 (MQ=255)
             gCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTgg                                                                                                                                                                                                                            <  1:584196/249‑1 (MQ=255)
                gggcGGCGTCTGTTGTTGTTTGGCAGCTTTGCCATGATGGCGTTGGCGGTGGCGGTGGTGGGGCTGTGCCTGGATATGGGGATTGTGATGGTGGTGATGGCCGTTGCGGGGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGGTGTGGCTCTATCCTAATGAACGCTTCCCGACAGATAGCCGCGCCTCTGCCGTGGGCGGGATTATTTCCGGAAGCCGTATTGGCACCATTGTGGCGACCTGGGCACt                                                                                                                                                                                                                       <  2:320848/248‑1 (MQ=255)
                    ggCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGGTCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTTTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATACTTATGAACTCTTTCCGACAGATATCCCCGCCTCTGCCCCGGGCGTGATTATGTCCTTAAGCCGTATTGGCCCCATTGGTTTGACCTGGGCACTTCCg                                                                                                                                                                                                                   >  1:5035/1‑251 (MQ=255)
                         ccATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCtt                                                                                                                                                                                                              <  1:367092/251‑1 (MQ=255)
                                    tGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGGGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGGAATTTGCAGTGGGTCTATCCTTATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGGGGGCGTGATTATTTCCTTAAGCCGTATTTGCCCCATTGTTTCGTCCCTGGCACTGCCGATCTTTTTTAATAAtt                                                                                                                                                                                                   >  1:114375/1‑251 (MQ=255)
                                             gggcaggatggcggtggcgctggcgGTTTTGGGGCTGATCACGGGTATGGGGATCTGGCTGGTGGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGTCTCTATCCTAATGAACTCTTCCTGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCGTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACTGTATc                                                                                                                                                                                           <  2:436444/236‑1 (MQ=255)
                                                   gatgaCGCGGGCGTTGGTGGTTTTGGGGCTGAAGCCGGATTTGGGGATTTGGCTGGTAGTGATGGCCTTTGCGTTGTATGCCTTTTTCTCTGGCGGGCCGGGTATTTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAaca                                                                                                                                                                                    <  2:465329/251‑1 (MQ=255)
                                                        cgctggcgctggCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATg                                                                                                                                                                               <  1:31775/251‑1 (MQ=255)
                                                                  ggcgGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATgggg                                                                                                                                                                      <  1:515448/250‑1 (MQ=255)
                                                                     gtgtttGGGGCTTATTCCGGATATGGGGTTCTGGTTGGGTGTTATGGCCTTTGCTGGGTAGGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCGAATGAACTCTTCCCTACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGTGGGCgg                                                                                                                                                                  <  2:404272/248‑1 (MQ=255)
                                                                           ggggCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTAtct                                                                                                                                                            >  1:176086/1‑251 (MQ=255)
                                                                              gctgttCCCGGATAGGGGGATCTGGTTGGTAGTGATGGCCGTGGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCTTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGc                                                                                                                                                         <  1:632369/246‑1 (MQ=255)
                                                                              gctgaTCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGCACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCCATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGc                                                                                                                                                         >  1:287504/5‑251 (MQ=255)
                                                                                  agttctgatctgggtgtCTGGATGGTAGTGGTGGTCTTTGTGGTGTTTGCCCTTTTCTCTGGCTGCCCTTATAATTTGCCGTGTCGCTATCTTCACGAGCTTTTCCCGTCTAATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTCTTGGTTCTATTGTTTCGACATGGGCACTGCCGATCTTTATGAACAAGTACGATATCAGTAGCACGATGCTAATGTCTGCGGGTAtctaactgt                                                                                                                                                      <  2:502694/234‑7 (MQ=255)
                                                                                                       gcgggttgttttggtcttggCGGTGTATGCCTTTTTCGTGTGCGGGTCGGGTAATTTGCTGTGGATCGATCATAAGGAACTCTTCCAGAAAGATTCCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAATGCGTATTGGCCCCATTGTTTCGACCTGGGCGCTGCCGATCTTTATCAATAATTGCTGTATCAGTATCACGATGCTAATGGGGGAGGGTATCTCGCTGTTTGGCTTGTTGATTTCTGTAg                                                                                                                                 <  2:213120/232‑1 (MQ=255)
                                                                                                          ggTAGTGATGGCCTTTTCGGTGTATGCCTTTTCCTCTGGCGGGCCGGGTAAGTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGt                                                                                                                                                                <  1:494381/217‑1 (MQ=255)
                                                                                                             aGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTg                                                                                                                           >  1:11697/1‑251 (MQ=255)
                                                                                                                gctGGCCATTGCGGTGTATGCCTTTTTCTCTGGCGGGACGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATAACCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGGATTGACTCCATTGTTTCGACCTGGGCACTGCCGATCCTTATCAATACTTACGGTATCAGTAACACGGTGCTAATGGGGGAggg                                                                                                                                                                 >  2:12168/3‑210 (MQ=255)
                                                                                                                gATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCggg                                                                                                                                                                 <  1:12168/210‑1 (MQ=255)
                                                                                                                   gtcgctgccggtggagtccTTGTTCTCTGGCGGGCGGGGGAATTGGCATGTGCTCTATCTCAAGGAACTCTTCCCGACAGATATCCGAGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCATCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTCCGGTAGCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCTTAGCGTTTGcgccc                                                                                                                      <  2:90872/233‑5 (MQ=255)
                                                                                                                         tGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTc                                                                                                               >  1:678129/1‑251 (MQ=255)
                                                                                                                               gtATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAgggg                                                                                                         >  1:304951/1‑250 (MQ=255)
                                                                                                                                  tGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGGATCTCGCTGTTTGGCTTGTTGATTTCCGTAACGTTTGCCCCGGGGGCTCGAGGGATgg                                                                                                      >  1:142978/1‑250 (MQ=255)
                                                                                                                                   gCCTTTTTCGCTGGCGGGCCTGGTAATTTGCAGTGGCTCGATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGt                                                                                                      <  1:632653/250‑1 (MQ=255)
                                                                                                                                                  cGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAAccc                                                                                       >  2:115025/1‑250 (MQ=255)
                                                                                                                                                   gggCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCCCTGGCGCAAAcc                                                                                        >  1:670789/1‑249 (MQ=255)
                                                                                                                                                       cGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCt                                                                                                                                                                             <  1:584956/160‑1 (MQ=255)
                                                                                                                                                       cGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGCTTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGGTGTCACTGGCGCAAACCAGCact                                                                                  >  2:268390/1‑249 (MQ=255)
                                                                                                                                                       cGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGAGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCt                                                                                                                                                                             >  2:584956/1‑160 (MQ=255)
                                                                                                                                                             aTTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCtt                                                                                                                                                                                                              <  2:587990/121‑1 (MQ=255)
                                                                                                                                                             aTTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCtt                                                                                                                                                                                                              >  1:587990/1‑121 (MQ=255)
                                                                                                                                                                    gTGGCGCTATCATAATGATCTTGTCCCGACAGTTATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGCCCTGGGCACTGCTGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGCCGTTCcgcg                                                                     <  2:243067/251‑1 (MQ=255)
                                                                                                                                                                     tGGCTCTATCCTAATGAACTCTTCCCGACAGAGATCCGAGCCTCTGCCGTGGGCGGGATTATGTCCTTAAGCCGGATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGAATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTTCCCCGGAGACTCGAAGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCgg                                                                    >  1:120723/1‑251 (MQ=255)
                                                                                                                                                                          tctATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAatat                                                                                >  1:104138/1‑235 (MQ=255)
                                                                                                                                                                          tctATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAatat                                                                                <  2:104138/235‑1 (MQ=255)
                                                                                                                                                                             aTCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCagag                                                              <  1:22015/250‑1 (MQ=255)
                                                                                                                                                                             aTCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAgaga                                                             <  1:445783/251‑1 (MQ=255)
                                                                                                                                                                                  aataatatcttcccaccagagttCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGGATTGGCACCATTGTGTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATggg                                                        <  2:520398/229‑1 (MQ=255)
                                                                                                                                                                                  aaTGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTTGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATgg                                                         <  1:614875/250‑1 (MQ=255)
                                                                                                                                                                                    tGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATg                                                          <  1:633316/247‑1 (MQ=255)
                                                                                                                                                                                    tGAACTCTTCCCAAAAGATATCCGCGCCTCTGCCGTGGGCGCGGTTATGTCCTTCAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGTTTTCCGTAACGTTTGACACGGAGAATCGAGGGATGTAACTGGCGCAAACCAGCACTATGACGACGCGCGGGCAGGGGATgggg                                                       >  2:633316/1‑250 (MQ=255)
                                                                                                                                                                                              ccGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCCGCCATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTc                                             >  2:506845/1‑250 (MQ=255)
                                                                                                                                                                                                      tatCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGATAATGGGGGCGGGt                                                                                                                                                                <  1:298522/127‑1 (MQ=255)
                                                                                                                                                                                                      tatCCGCGCATCTGCCGTGGGCGAGATTATGTCATTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCTATAATTACGGTATCAGTAACCCGATGCTAATGGGGGCGGGt                                                                                                                                                                >  2:298522/1‑127 (MQ=255)
                                                                                                                                                                                                       atCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTtctct                                   >  1:487794/1‑251 (MQ=255)
                                                                                                                                                                                                                      gTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAat                    >  2:483332/1‑251 (MQ=255)
                                                                                                                                                                                                                      gTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAat                    >  1:239336/1‑251 (MQ=255)
                                                                                                                                                                                                                       tGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGCGACTCGAgg                                                                                                           >  2:344893/1‑163 (MQ=255)
                                                                                                                                                                                                                       tGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAgg                                                                                                           <  1:344893/163‑1 (MQ=255)
                                                                                                                                                                                                                        gggCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAata                   <  1:569073/250‑1 (MQ=255)
                                                                                                                                                                                                                         ggCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTCTCAATAATTACGGAATCAGTAACACGATGCTAATGTGGGCGGGGATCTCGCAGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGTGTCTCGAGGGATGTCACTGGCGAAACCCAGAATTATGAACATCCGAGGGCTGGGATTGGGGGAATTTTTTCAGGTTTCTCCCTTTTCTTAAtacataat                 >  2:619842/1‑244 (MQ=255)
                                                                                                                                                                                                                                         tAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCcg  <  1:483332/250‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
GTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCG  >  NC_000913/46709‑47188

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