Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,748 G→A S404S (AGC→AGT)  yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009136,7480GA97.6% 116.7 / ‑6.0 41S404S (AGC→AGT) yaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/22);  minor base T (1/0);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTG  >  NC_000913/6505‑6988
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
gaatgagTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCa                                                                                                                                                                                                                                           <  2:168567/251‑1 (MQ=255)
          aTCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAAACTTAGCAAGGGgccgcc                                                                                                                                                                                                                                 >  1:244801/1‑251 (MQ=255)
                 gTTGTGATTACTCCTGCGAAACATCATCCCACGCGCCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGGCGATGAc                                                                                                                                                                                                                           <  1:165391/250‑1 (MQ=255)
                  ttGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAACCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCg                                                                                                                                                                                                                         >  1:278307/1‑251 (MQ=255)
                                   aaaCATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGA‑CCCATCAGCGGAAGACTTAGCAAGGGGCCGCCGATGACCGTTGCGagggggc                                                                                                                                                                                                             >  1:274271/1‑240 (MQ=255)
                                     aCATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGCTATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCCAATTGGTAATAGCCATGCCGGCCCATTA‑GTTATTTGCCCGTTGCACACACCA‑CGGAAAACTTACCAAGGTGCACCCGATGACCGTTGCGaaaggggagaaac                                                                                                                                                                                                       >  2:313547/1‑239 (MQ=255)
                                              cACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGACAAATTGGTAATAGCCAAGCCGGCC‑ATTTTGATATCTGCCAGTTGC‑CACATCAGCGGGAGAATTAGAAAGGGGCCCCCGGTGACCGTTGccgaggggccggaccgccaaaaac                                                                                                                                                                                              >  2:421674/1‑227 (MQ=255)
                                                 gcgTCCGGAGAAAGCTGGCTGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCCAATTGGAAATACCCATTCATGTC‑CTTATGATATCTGCCAGTTGA‑CACATCAGCGGAAGACTTAGCGAGGTGCCGCCGATGACCGTTGCGCAGGTGACGATCCGCAAACAcc                                                                                                                                                                                              >  1:115176/1‑247 (MQ=255)
                                                        gagaAAGCTGGCGGCCGATATCCGGATAACGAAACGGATCAAACACCGGGCGCACGCCTAGGTTACGCTGGCGGAGATAATCACTGGCAATGGTAGGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTGAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTa                                                                                                                                                                                   <  1:496934/251‑1 (MQ=255)
                                                          gaAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGCGGGCAGTAAAATGGCGGTCAAATTGGTAATAGCCCTGCCAGCC‑ATTATGATATCGGCCGATTGC‑CACCTCAACGGAAAACTTCGCAAGCAGCCGCCGAAGACCGTTGCGACGGGGCAGATCCGAAAACACCAGATCGccttggg                                                                                                                                                                                 >  1:207718/1‑245 (MQ=255)
                                                          gaAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCCAACACCAGATCGCTTTAgg                                                                                                                                                                                 >  1:146504/1‑251 (MQ=255)
                                                                 ggcggcCGATATCCGGATAACGCAACGGATCAAACACCGGGCCCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTc                                                                                                                                                                          <  1:427365/251‑1 (MQ=255)
                                                                  gcggccGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCa                                                                                                                                                                         <  1:361885/247‑1 (MQ=255)
                                                                          tatCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAAGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGGGCAGATCCGCAAACACCAGATGGCTTTAGGGTTGGTCAGGCGtca                                                                                                                                                                  >  1:152408/1‑248 (MQ=255)
                                                                                      gCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGGGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGAttt                                                                                                                                                     >  1:122347/1‑250 (MQ=255)
                                                                                               cAAACCCCGGACGCACGCCGAGTGTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTGTCGGc                                                                                                                                             <  2:341734/250‑1 (MQ=255)
                                                                                                     ccGGGCTCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATg                                                                                                                                      <  1:418042/251‑1 (MQ=255)
                                                                                                         gcgcACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGt                                                                                                                                  <  1:349070/251‑1 (MQ=255)
                                                                                                                 cGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACg                                                                                                                          <  1:393761/251‑1 (MQ=255)
                                                                                                                         cGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCCAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGATGGCAACGCTGGAGCt                                                                                                                  >  1:182688/1‑251 (MQ=255)
                                                                                                                                gTAGATAATCTCTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCAAGGCAGGCC‑ATGATGATATCTGCCAGTTGC‑CACACCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCa                                                                                                           <  1:115300/251‑1 (MQ=255)
                                                                                                                                          acgggcaagggtacaacccACAGGCGAGAGCAGTAAAATGGCGTTCAAATGGGTAATGGCCATGCAGACC‑ATTATGAGAGCTGCCAGTGGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGAACCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAATGGGCAACGATGGAGCTGAAGGCAAACAGAATaa                                                                                                 <  2:34701/234‑1 (MQ=255)
                                                                                                                                              gCAAGGGTATGACCCACAGGCGAGAGCAGGAAAATGGCGGTCAACTGGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTGGC‑CACATCAGCGGAAGACTTAGCAAGTTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGTTTTAGGGTTGTTCAGGCGGAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGATCTGAAGGCAAACAGAATAAcccc                                                                                             <  2:47108/251‑3 (MQ=255)
                                                                                                                                              gCAAGGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACcac                                                                                             <  2:300264/251‑1 (MQ=255)
                                                                                                                                               cAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACcacg                                                                                            <  1:421674/251‑2 (MQ=255)
                                                                                                                                                   ggTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAg                                                                                                                    <  1:154632/223‑1 (MQ=255)
                                                                                                                                                    gtagtaaCCACAGGCGAGAGCATTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACAGCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGt                                                                                       <  2:122347/246‑1 (MQ=255)
                                                                                                                                                              caGGCGAGAGCAGGAAAATGGCGGGCAAATTGGGAATAGCCAGGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGACCGCTTTAGGGTTGTTCAGGCGGAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCa                                                                             <  2:278307/251‑1 (MQ=255)
                                                                                                                                                                       gCAGTAAAATGGCGGTCAAATTGGTAATAGCCCTGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGGAGGGGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGGTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACGAACTCAGcccccc                                                                      >  2:183454/1‑249 (MQ=255)
                                                                                                                                                                         aGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGTGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCaa                                                                                                          <  1:70702/211‑1 (MQ=255)
                                                                                                                                                                            aaaaTGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAAcc                                                               <  1:218163/251‑1 (MQ=255)
                                                                                                                                                                                 ggtggTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTACATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACACCCGGAACCCAtt                                                           >  2:153228/4‑250 (MQ=255)
                                                                                                                                                                                            gggTAATAGCCAGGCAGGCC‑AGTATGATATCTGCCAGTTGC‑CACAGCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACTATGGAGCTGAAGGCAAACAGAATAACCACGATGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATc                                               <  2:182688/250‑1 (MQ=255)
                                                                                                                                                                                              gTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGc                                             >  1:377286/1‑251 (MQ=255)
                                                                                                                                                                                               tAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGcc                                            >  1:60518/1‑251 (MQ=255)
                                                                                                                                                                                                   aGCCATGCAGGAC‑ATAATGATATCTGCCAGTTGC‑CACATCAGCGGAAGTCATAGCAAGGTGACGCCGATGGACGTTGCGAAGGTGCAGATCCGCAAACCCCAGATCGCTTTAGGGTTGTTCAGGCGTAGAAAGAAGAGAATGTTTTCCGCGTAAATGTTGCTGGCAACGATGGAGCTGAAGCCATACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGAAACACATTAACACCCGCAtagacctcc                                        >  1:455527/1‑243 (MQ=255)
                                                                                                                                                                                                             gCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTg                              <  1:269406/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTGATATCTGCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACCCCCGCATCGCCTTCTGGATAAGCTGAATAc                         <  1:258388/250‑1 (MQ=255)
                                                                                                                                                                                                                             gCCAGTTGC‑CACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCgg               >  2:190732/1‑251 (MQ=255)
                                                                                                                                                                                                                                     c‑cACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCACACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGACTTCTGGATAAGCTTAATACCTTCCAGCGGCATGTAgg       >  1:114751/1‑251 (MQ=255)
                                                                                                                                                                                                                                           tCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCCTGTAGGTtgtg  >  2:388151/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
GAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCC‑ATTATGATATCTGCCAGTTGC‑CACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTG  >  NC_000913/6505‑6988

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: ‑