Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,304 A→G Y73Y (TAT→TAC ettA ← energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,3040AG90.0% 97.2 / ‑0.3 43Y73Y (TAT→TACettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (15/21);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAA  >  NC_000913/4630058‑4630554
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
gTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAAc                                                                                                                                                                                                                                                         <  1:173131/251‑1 (MQ=255)
   cgtcAGACGTTGCGCTTCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGAACTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTGCAGCTGCGGTTCCTGCGGCAGGTAACCaa                                                                                                                                                                                                                                                      <  2:675243/248‑1 (MQ=255)
     gCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACAACTTCAGAAACCGCTTCTTCACTGGACTAACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCACGAAACCAATc                                                                                                                                                                                                                                                    >  2:551516/1‑251 (MQ=255)
      cAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCt                                                                                                                                                                                                                                                   <  1:12626/251‑1 (MQ=255)
             aGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTc                                                                                                                                                                                                                                            <  1:504556/251‑1 (MQ=255)
               cgcATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTg                                                                                                                                                                                                                                          <  1:58905/251‑1 (MQ=255)
                 cATCCGCCGCACGCTCCAGCTGTACGTTCATGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCTCTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTGCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGc                                                                                                                                                                                                                                        <  2:563724/251‑1 (MQ=255)
                        cgcACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGAc                                                                                                                                                                                                                                 >  2:113674/1‑251 (MQ=255)
                             catcctgCTGTACCTTCCGGTTATGACCGTCGTGCGCCTGAATGATCTCTTCCAGACGGGCTTGTTCAGCGGACATATTTTCAAAAACGTCTACCGGATCGGCGTACAGCGCATACTCTTCATCCAGGCGTTTCAGGGCGTTAACCACCTCAGAAACCGCTTCTTCAATGGACGCACGCTCGGTGTGTTCCGGGTTCAGCTGCGGTCCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTTCGGTCGCGCt                                                                                                                                                                                                                            <  2:245645/245‑1 (MQ=255)
                                   gCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCt                                                                                                                                                                                                                      <  1:51222/251‑1 (MQ=255)
                                      gTACGTTCAGGTTATGACCGTCGTGAGCCTGAATTATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGTCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCtt                                                                                                                                                                                                                     <  1:430502/249‑1 (MQ=255)
                                           ttCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTc                                                                                                                                                                                                              >  1:605704/1‑251 (MQ=255)
                                             caggttATGACCGTCGTGATCTTTAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACCCGTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCCCGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCtt                                                                                                                                                                                                            <  2:124439/247‑1 (MQ=255)
                                             cagggtttgacctcccggAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCTGCCATATTGTCAAAATCGGCATCCGGCTCGGCGTACAGCGCATTCACTTAATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGGTTCTTCAATGGACTTACGCTCGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCtt                                                                                                                                                                                                            <  2:74999/234‑1 (MQ=255)
                                                   aTGATCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCa                                                                                                                                                                                                       >  1:36638/1‑250 (MQ=255)
                                                   aTGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCa                                                                                                                                                                                                       <  1:380473/250‑1 (MQ=255)
                                                    tGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGATCCCATCTTGATGTCTTGCTGCGAGCGCGCtgtccctttggattctttcacca                                                                                                                                                                                                      >  2:626187/1‑228 (MQ=255)
                                                       ccgccgtGTGCCTTAATGATCTCTCCCAGACGGCTTGGTTCAGCGGCCAGCTTGTCAAAATCTGCATCCGGATCGGCGTACAGCGCATACCCTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACTCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGcc                                                                                                                                                                                                  <  2:331740/247‑1 (MQ=255)
                                                        cggcctgatcccgAATGATCGCTTCCAGACGGCCTTGTTCAGCGGCCATCTTGTCAAAATCTGCATCCGTATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACGTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCTGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGccc                                                                                                                                                                                                 <  2:244712/239‑1 (MQ=255)
                                                                       tGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCgcagc                                                                                                                                                                                 >  1:280636/1‑251 (MQ=255)
                                                                                  aGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCt                                                                                                                                                                                                                                           >  1:556361/1‑183 (MQ=255)
                                                                                  aGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAACCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCt                                                                                                                                                                                                                                           <  2:556361/183‑1 (MQ=255)
                                                                                                     ccaccTTGTCAAACCCGGCACCCGGACCGGCGTACGGCGCATCCCCTCCACCCGGGCGTTCCAGGGGTTTACCCACTCCAGAACCCGCTTCTTCAAGGGCCTCACGCAGGGTGTGTTCCGGGTTCAGCTGCGTTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGCCGCGCTTCCCCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGAccc                                                                                                                                                   <  1:531106/247‑1 (MQ=255)
                                                                                                            gTCAAAATCGGCCTCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCATTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAAACAATCTTGATGTCTGGCTGCGGCCGCGCCTCACCTTCGATGTCTTTCTCAATCCCCCCCATAATGCGACGCAGGGTAGAATTGCCCGCGCCGTTCAGAcaccgggccc                                                                                                                                            >  2:202344/1‑242 (MQ=255)
                                                                                                                            ggtgcGGCGTACAGCGCATACCCTTCACCCAGGCTTTTCAGGGCTTTAACCACTTCAAAACCCGTTTCTTCAATGGACTCCCGCACGGTGTGTTCCGGGTTCAGCGGCGTTTCCTGCGGCAGTTAACCAATCTTTATGTCTGGCGGCGGCCCCGTTTCCCCTTCGATTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAgg                                                                                                                            <  1:202344/247‑1 (MQ=255)
                                                                                                                               ccGGCGTACAGCGCAACCCCTTCATCCAGGCGTTCCAGGGCGTTACCCCCTCCAAAACCCTCTTCTTCAATGGACCCCCGCAGGGTGTGTTCCGGGTTCAGCTGCGTTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAAGGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGgaa                                                                                                                         <  1:618337/250‑1 (MQ=255)
                                                                                                                                    gTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGAAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAAc                                                                                                                    >  1:55247/1‑251 (MQ=255)
                                                                                                                                     tACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCATTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGCAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGACGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTCGACATTCCCCCACCCATCAAACCCAGGCCCCCAAATTTTACCCCAggggaaaaacc                                                                                                                   >  2:363447/1‑243 (MQ=255)
                                                                                                                                       cAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAACCAATCTTGATGTCTGGCTGCGGAAGCGCTTCACCTTCGATGTCTTTATTAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCCCCGTTCAGACCCAGGACACCCATTTTTTCCCCCGGGGAGAAACTcc                                                                                                                 >  2:329570/1‑250 (MQ=255)
                                                                                                                                             aTACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGGCTTTTTCAATGCCCGCCATAATGCGCACCAGGGTAGGCTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCAAGGGAAGAAACTCcgag                                                                                                              >  1:275853/1‑248 (MQ=255)
                                                                                                                                             aTACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAAACAATCTTGATGTCTGGCTGCGGACGCGCCTCACCTTCGGTGGCTTTATCAATGCCCGCCATAATGCGCCGCAGGGTAGACTTACCCCCGCCGTTCAGGCCCAGGACACCAATTTTTGTCCCAGGGAAGaaaatccgagggg                                                                                                           >  2:342870/1‑241 (MQ=255)
                                                                                                                                                  cTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt                                                                                                      >  1:141965/1‑251 (MQ=255)
                                                                                                                                                    ccATCCAGGCGTTTAAGGGCGTTAACCCCTACAGAACCCGCTTCTTCATGGGCCTCACGCGCGGTGTGTTCGGGTTTCGGCTGCGGTTCCTGCGGCAGATAACCATTCTTGAGGTCTGGCTGCGGACGCGCTTCACCTTGGATGTTTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGtttt                                                                                                     <  1:328420/249‑1 (MQ=255)
                                                                                                                                                    ccATCCAGGCGTTACAGGGCGTAAACCGCTTAAGAAACCTCTTCTTCAGGGGCCTCACGCACGGTGTGTCGCGGGTTCAGCTGGGGTTCCTGGGGCAGGTAACCAATCTGTATGTCGGGCTGCGGCCGCGCTTCACCTTCGTTGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt                                                                                                      <  1:367624/248‑1 (MQ=255)
                                                                                                                                                               tttCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCCATGCCCGCCATAATGCgcagca                                                                                                                                                                                >  1:517519/1‑164 (MQ=255)
                                                                                                                                                                ttcggggCGTAAACCACTTAAGAACCCGTTTCTTCAATGGACTCCCGCCGGGTTTGTTCCGGTTTCAGCTGCGTTCCCGGCGGCAGGTAACAAATTTTGGTGTCGGGCTGCGGACGCGCTCCCCCTCCGATGTCTTTATCAAGGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGAc                                                                                         <  1:631452/246‑1 (MQ=255)
                                                                                                                                                                                  taaaaaaCCGCTTCTTCAATGGACTCCCCCCGGGTGTGTCCCGGGTTCAGCGGCGGTTCCTGCGGCAGTTACCCAATTTTGATGTCTGGCGGCGGCCGCCCTTCCCCTTCGATGTCTTAATCAAGGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACtt                                                                       <  1:613994/246‑1 (MQ=255)
                                                                                                                                                                                       aaCCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCg                                                                  >  1:598542/1‑250 (MQ=255)
                                                                                                                                                                                                          ccccgcgcggggtgtTCGGGGTTCCGCTGGGGTTCCGGCGGCAGGTACCCATTCTTGATGTCTGGCTGCGGACGCGCTTCCCCTGCGATGTCTTTACAAAGCCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCACGGTATAAAc                                              <  1:605224/240‑1 (MQ=255)
                                                                                                                                                                                                                     ggtttCCGGTTTCATCTGCGTTCCCTGCGGCAGGTACCCAATCTTGTTGTCTGGCGGCGGCCCCTCTTCCCCTTCAATGTCTTTACCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCa                                   <  1:656151/248‑1 (MQ=255)
                                                                                                                                                                                                                             gtttCAGCGGCGGTCCCGGGGGCAGTTACCCAATCTTGATGTCGGGCTGCGGCCGCTCTTCCCCTTCGATGTTTTTATAAAGCCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGAc                           <  1:511637/249‑1 (MQ=255)
                                                                                                                                                                                                                               ttCAGCTGCGGTTTCTGCGGCAGGTAACCAATCTTGCTGTCTGGCTGAGGTCGCGCCTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCCCGCCCTTCAGACCCAAGACACCAATTTTTGCCCCAGGGCAGAAACTCAAAGAGATGTTTTTCAAAATATGACCTTTCGGCGGAACAACTTTGCCGACACGGTGAGTGATATTAACGCATTACGCCACGTTGGACtt                         >  2:360874/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gCAGGTAACCAATCTTGATGTATGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGGGCCACGTTGGACTTCGCCTCTTTGTTTATCGTg      >  2:85201/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       aTAACCAATCTTGTTGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATaa  <  1:113674/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
GTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAA  >  NC_000913/4630058‑4630554

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