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→TAC) ettAenergy‑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

GTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAA  >  NC_000913/4630058‑4630554
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
gTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAAc                                                                                                                                                                                                                                                         <  1:173131/251‑1 (MQ=255)
   cgtcAGACGTTGCGCTTCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGAACTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTGCAGCTGCGGTTCCTGCGGCA‑GGTAACCaa                                                                                                                                                                                                                                                      <  2:675243/248‑1 (MQ=255)
     gCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACAACTTCAGAAACCGCTTCTTCACTGGACTAACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCACGA‑AACCAATc                                                                                                                                                                                                                                                    >  2:551516/1‑251 (MQ=255)
      cAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCt                                                                                                                                                                                                                                                   <  1:12626/251‑1 (MQ=255)
             aGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTc                                                                                                                                                                                                                                            <  1:504556/251‑1 (MQ=255)
               cgcATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTg                                                                                                                                                                                                                                          <  1:58905/251‑1 (MQ=255)
                 cATCCGCCGCACGCTCCAGCTGTACGTTCATGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCTCTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTGCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGc                                                                                                                                                                                                                                        <  2:563724/251‑1 (MQ=255)
                        cgcACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GATAACCAATCTTGATGTCTGGCTGCGGAc                                                                                                                                                                                                                                 >  2:113674/1‑251 (MQ=255)
                             catcctgCTGTACCTTCCGGTTATGACCGTCGTGCGCCTGAATGATCTCTTCCAGACGGGCTTGTTCAGCGGACATATTTTCAAAAACGTCTACCGGATCGGCGTACAGCGCATACTCTTCATCCAGGCGTTTCAGGGCGTTAACCACCTCAGAAACCGCTTCTTCAATGGACGCACGCTCGGTGTGTTCCGGGTTCAGCTGCGGTCCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTTCGGTCGCGCt                                                                                                                                                                                                                            <  2:245645/245‑1 (MQ=255)
                                   gCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCt                                                                                                                                                                                                                      <  1:51222/251‑1 (MQ=255)
                                      gTACGTTCAGGTTATGACCGTCGTGAGCCTGAATTATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGTCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCtt                                                                                                                                                                                                                     <  1:430502/249‑1 (MQ=255)
                                           ttCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTc                                                                                                                                                                                                              >  1:605704/1‑251 (MQ=255)
                                             caggttATGACCGTCGTGATCTTTAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACCCGTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCCCGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCtt                                                                                                                                                                                                            <  2:124439/247‑1 (MQ=255)
                                             cagggtttgacctcccggAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCTGCCATATTGTCAAAATCGGCATCCGGCTCGGCGTACAGCGCATTCACTTAATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGGTTCTTCAATGGACTTACGCTCGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCtt                                                                                                                                                                                                            <  2:74999/234‑1 (MQ=255)
                                                   aTGATCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCa                                                                                                                                                                                                       >  1:36638/1‑250 (MQ=255)
                                                   aTGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCa                                                                                                                                                                                                       <  1:380473/250‑1 (MQ=255)
                                                    tGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGGATCCCATCTTGATGTCTTGCTGCGAGCGCGCtgtccctttggattctttcacca                                                                                                                                                                                                      >  2:626187/1‑228 (MQ=255)
                                                       ccgccgtGTGCCTTAATGATCTCTCCCAGACGGCTTGGTTCAGCGGCCAGCTTGTCAAAATCTGCATCCGGATCGGCGTACAGCGCATACCCTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACTCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGcc                                                                                                                                                                                                  <  2:331740/247‑1 (MQ=255)
                                                        cggcctgatcccgAATGATCGCTTCCAGACGGCCTTGTTCAGCGGCCATCTTGTCAAAATCTGCATCCGTATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACGTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCTGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGccc                                                                                                                                                                                                 <  2:244712/239‑1 (MQ=255)
                                                                       tGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCgcagc                                                                                                                                                                                 >  1:280636/1‑251 (MQ=255)
                                                                                  aGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCt                                                                                                                                                                                                                                           >  1:556361/1‑183 (MQ=255)
                                                                                  aGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAACCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCt                                                                                                                                                                                                                                           <  2:556361/183‑1 (MQ=255)
                                                                                                     ccaccTTGTCAAACCCGGCACCCGGACCGGCGTACGGCGCATCCCCTCCACCCGGGCGTTCCAGGGGTTTACCCACTCCAGAACCCGCTTCTTCAAGGGCCTCACGCAGGGTGTGTTCCGGGTTCAGCTGCGTTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGCCGCGCTTCCCCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGAccc                                                                                                                                                   <  1:531106/247‑1 (MQ=255)
                                                                                                            gTCAAAATCGGCCTCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCATTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAAACAATCTTGATGTCTGGCTGCGGCCGCGCCTCACCTTCGATGTCTTTCTCAATCCCC‑CCCATAATGCGACGCAGGGTAGAATTGCCCGCGCCGTTCAGAcaccgggccc                                                                                                                                            >  2:202344/1‑242 (MQ=255)
                                                                                                                            ggtgcGGCGTACAGCGCATACCCTTCACCCAGGCTTTTCAGGGCTTTAACCACTTCAAAACCCGTTTCTTCAATGGACTCCCGCACGGTGTGTTCCGGGTTCAGCGGCGTTTCCTGCGGCA‑GTTAACCAATCTTTATGTCTGGCGGCGGCCCCGTTTCCCCTTCGATTTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAgg                                                                                                                            <  1:202344/247‑1 (MQ=255)
                                                                                                                               ccGGCGTACAGCGCAACCCCTTCATCCAGGCGTTCCAGGGCGTTACCCCCTCCAAAACCCTCTTCTTCAATGGACCCCCGCAGGGTGTGTTCCGGGTTCAGCTGCGTTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAAGGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGgaa                                                                                                                         <  1:618337/250‑1 (MQ=255)
                                                                                                                                    gTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGAAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAAc                                                                                                                    >  1:55247/1‑251 (MQ=255)
                                                                                                                                     tACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCATTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGCAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGACGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTCGACATTCCCCCACCCATCAAACCCAGGCCCCCAAATTTTACCCCAggggaaaaacc                                                                                                                   >  2:363447/1‑243 (MQ=255)
                                                                                                                                       cAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGGAACCAATCTTGATGTCTGGCTGCGGAAGCGCTTCACCTTCGATGTCTTTATTAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCCCCGTTCAGACCCAGGACACCCATTTTTTCCCCCGGGGAGAAACTcc                                                                                                                 >  2:329570/1‑250 (MQ=255)
                                                                                                                                             aTACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGGCTTTTTCAATGCCC‑GCCATAATGCGCACCAGGGTAGGCTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCAAGGGAAGAAACTCcgag                                                                                                              >  1:275853/1‑248 (MQ=255)
                                                                                                                                             aTACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGGAAACAATCTTGATGTCTGGCTGCGGACGCGCCTCACCTTCGGTGGCTTTATCAATGCCC‑GCCATAATGCGCCGCAGGGTAGACTTACCCCCGCCGTTCAGGCCCAGGACACCAATTTTTGTCCCAGGGAAGaaaatccgagggg                                                                                                           >  2:342870/1‑241 (MQ=255)
                                                                                                                                                  cTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt                                                                                                      >  1:141965/1‑251 (MQ=255)
                                                                                                                                                    ccATCCAGGCGTTTAAGGGCGTTAACCCCTACAGAACCCGCTTCTTCATGGGCCTCACGCGCGGTGTGTTCGGGTTTCGGCTGCGGTTCCTGCGGCA‑GATAACCATTCTTGAGGTCTGGCTGCGGACGCGCTTCACCTTGGATGTTTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGtttt                                                                                                     <  1:328420/249‑1 (MQ=255)
                                                                                                                                                    ccATCCAGGCGTTACAGGGCGTAAACCGCTTAAGAAACCTCTTCTTCAGGGGCCTCACGCACGGTGTGTCGCGGGTTCAGCTGGGGTTCCTGGGGCA‑GGTAACCAATCTGTATGTCGGGCTGCGGCCGCGCTTCACCTTCGTTGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt                                                                                                      <  1:367624/248‑1 (MQ=255)
                                                                                                                                                               tttCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCCATGCCC‑GCCATAATGCgcagca                                                                                                                                                                                >  1:517519/1‑164 (MQ=255)
                                                                                                                                                                ttcggggCGTAAACCACTTAAGAACCCGTTTCTTCAATGGACTCCCGCCGGGTTTGTTCCGGTTTCAGCTGCGTTCCCGGCGGCA‑GGTAACAAATTTTGGTGTCGGGCTGCGGACGCGCTCCCCCTCCGATGTCTTTATCAAGGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGAc                                                                                         <  1:631452/246‑1 (MQ=255)
                                                                                                                                                                                  taaaaaaCCGCTTCTTCAATGGACTCCCCCCGGGTGTGTCCCGGGTTCAGCGGCGGTTCCTGCGGCA‑GTTACCCAATTTTGATGTCTGGCGGCGGCCGCCCTTCCCCTTCGATGTCTTAATCAAGGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACtt                                                                       <  1:613994/246‑1 (MQ=255)
                                                                                                                                                                                       aaCCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCg                                                                  >  1:598542/1‑250 (MQ=255)
                                                                                                                                                                                                          ccccgcgcggggtgtTCGGGGTTCCGCTGGGGTTCCGGCGGCA‑GGTACCCATTCTTGATGTCTGGCTGCGGACGCGCTTCCCCTGCGATGTCTTTA‑CAAAGCCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCACGGTATAAAc                                              <  1:605224/240‑1 (MQ=255)
                                                                                                                                                                                                                     ggtttCCGGTTTCATCTGCGTTCCCTGCGGCA‑GGTACCCAATCTTGTTGTCTGGCGGCGGCCCCTCTTCCCCTTCAATGTCTTTACCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCa                                   <  1:656151/248‑1 (MQ=255)
                                                                                                                                                                                                                             gtttCAGCGGCGGTCCCGGGGGCA‑GTTACCCAATCTTGATGTCGGGCTGCGGCCGCTCTTCCCCTTCGATGTTTTTATAAAGCCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGAc                           <  1:511637/249‑1 (MQ=255)
                                                                                                                                                                                                                               ttCAGCTGCGGTTTCTGCGGCA‑GGTAACCAATCTTGCTGTCTGGCTGAGGTCGCGCCTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCCCGCCCTTCAGACCCAAGACACCAATTTTTGCCCCAGGGCAGAAACTCAAAGAGATGTTTTTCAAAATATGACCTTTCGGCGGAACAACTTTGCCGACACGGTGAGTGATATTAACGCATTACGCCACGTTGGACtt                         >  2:360874/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gCA‑GGTAACCAATCTTGATGTATGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGGGCCACGTTGGACTTCGCCTCTTTGTTTATCGTg      >  2:85201/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       aTAACCAATCTTGTTGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATaa  <  1:113674/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
GTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCA‑GATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAA  >  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: ‑