Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 266,523 G→A L11L (CTG→TTG)  yafY ← inner membrane lipoprotein YafY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,5230GA97.4% 125.9 / ‑3.6 40L11L (CTG→TTG) yafYinner membrane lipoprotein YafY
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (19/19);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

CACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCGGACACCCCGGTTTTGTTCAGAATGATGACGATGACAACAATGGC  >  NC_000913/266290‑266768
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
cACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCTGCAAGGGGGGGCCGGATAAGTAGTCAACCTGTGACCCTTTTTTGCTGACATCATCCTGGGGGCTGACCACATCACTGAAAAGCGGGGGGGGGGGGGCCTTTAATTTACTGATGGCCTTCAGGGCATCACTTCTTTCATCGCCTGCGGTGAGAAACAGAGTGGGGGGAACCAACGCCAGCCGAgggaggg                                                                                                                                                                                                                                       >  1:333394/1‑246 (MQ=255)
     ggTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTa                                                                                                                                                                                                                                  <  2:333394/251‑1 (MQ=255)
                    cATCCTCGCAGCCATCGGGCGACTGGCTGGTGCGTAAAACCTGCAAGGTGGTACCTGATAAGTACGCATCCTGTGACCATTCTTCGCTGACATCATCCTGGTTGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGTAACAGAGTGTTGGTAACCAACCCCAGCAGAGGCAGTGTTTTCCCTGTCATTAttttt                                                                                                                                                                                                                   <  2:31902/251‑1 (MQ=21)
                           gCAGCCATCGGGCGACTGGCGGGTAAGGAAAACCCGCAAGGGGGGGCCGGAGAAGTAGTACACCCGTGACCCTTCTTTGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGGGGGGGGGGGGCCTTAAATTTACTGATGGCCTTCAGGTCCTCACTTCTGTCCTCGCATGCGATGAGAAACAGAGTGGGGGGAACCCAAGCCCGCAGAGGCAGGGGTTTTCGTTTTATTATTTTTTTTCTg                                                                                                                                                                                                            >  1:140775/1‑250 (MQ=255)
                               ccATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATc                                                                                                                                                                                                       <  1:351014/251‑1 (MQ=255)
                                cATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCa                                                                                                                                                                                                      <  2:140775/251‑1 (MQ=255)
                                    gggCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCATCCTGTGACCATTCTTCGTTGACATCATCCTGGTTGCTGACCACATCACTGAAACGCGGGGGGTTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACg                                                                                                                                                                                                  <  2:108204/251‑1 (MQ=255)
                                        gACTGGCGGGTACGTAAAAACTGCCAGGGGGGACCGGATAAATAGTCAACCTGTGGCCATTCTTCGCTGACATCAACCTGGTGGGTGACCACAACACTGAAACGCGGGGGGGTGGGGGCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAAAGAGTGGGGGGAACCAAAGCCCGCAGAGGGAGGGTTTTACGTTTTATTTTTTTTTTTCTGAAATCAGACggaac                                                                                                                                                                                              >  2:79411/1‑247 (MQ=255)
                                              cGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTGCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTTAAATCAGACGAACCACTTTg                                                                                                                                                                                        <  2:230695/251‑1 (MQ=255)
                                                ggTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGc                                                                                                                                                                                      <  1:79411/251‑1 (MQ=255)
                                                    cGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAg                                                                                                                                                                                  <  1:106146/251‑1 (MQ=255)
                                                               cAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGcc                                                                                                                                                                       <  2:406839/251‑1 (MQ=255)
                                                                                        aCCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCcg                                                                                                                                              >  1:34202/1‑251 (MQ=255)
                                                                                                 cATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGgt                                                                                                                                     >  1:295422/1‑251 (MQ=255)
                                                                                                         gCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATTCCCACAAAAATACCCACCGGCACGCCGACGAACGGCGTCAgcg                                                                                                                              >  2:91647/1‑251 (MQ=255)
                                                                                                          cTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGcgcg                                                                                                                            >  2:299628/1‑251 (MQ=255)
                                                                                                                    ccTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAAGCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATGCCCACCGGCACGCCGACGAACGGTGTCAGCGCGTCACTGgcag                                                                                                                  <  2:2634/251‑1 (MQ=255)
                                                                                                                          ggCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCgg                                                                                                            >  1:288641/1‑251 (MQ=255)
                                                                                                                                ccACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTgcgc                                                                                                      >  1:20269/1‑251 (MQ=255)
                                                                                                                                                                   aaTTTACTGATTGCCTTCAGGTCAGCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGGTTCAGTATGTTGTTCCTTAAATCAGACGAACCACTTTGGCAAAGACATAAATTCCCACGAAAATACCCACCTGCACGCCGACGAACGGTGTCAGCTTGACACTGTGAGCACCGGCTGTGCCACTGCCCGTCAGCAGTGcttgagttggggcaagg                                                                   <  2:167875/251‑17 (MQ=17)
                                                                                                                                                                       tACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGCGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCACGGGTc                                                                >  1:340588/1‑250 (MQ=255)
                                                                                                                                                                           gATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGc                                                           >  1:109953/1‑251 (MQ=255)
                                                                                                                                                                              ggCCTTCAGGTCATCACTTCTGTCATCGCATGCGAAGAGAAACAGAGCGGTAGCAACCAACCCCAGCCGAGGCAGAGTTTTACATTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAACGACATAAATGCCCACGAAAATACCCACCGCCACGCCGACGAACGGTGTCAGCGCGACACTGGACGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCcgc                                                        >  1:177575/1‑251 (MQ=255)
                                                                                                                                                                               gCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCcgcg                                                       <  2:295422/251‑1 (MQ=255)
                                                                                                                                                                                 cTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCcgcg                                                       >  1:291825/1‑249 (MQ=255)
                                                                                                                                                                                      ggTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGt                                                <  2:159167/251‑1 (MQ=255)
                                                                                                                                                                                       gTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTc                                               <  2:73989/251‑1 (MQ=255)
                                                                                                                                                                                            caCTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCGGac                                           <  1:242507/250‑1 (MQ=255)
                                                                                                                                                                                                   gTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCGGACACCCCgg                                    <  2:355095/250‑1 (MQ=255)
                                                                                                                                                                                                   gTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCGGACACCCCg                                     <  2:384777/249‑1 (MQ=255)
                                                                                                                                                                                                         gCATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCGGACACCCCGGTgttgtt                             >  2:53016/1‑251 (MQ=255)
                                                                                                                                                                                                          caggCGATGAGAAACAGAGTGGTGGCAACCACCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGCGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGGCGGACACCCCGGTTTTGTTc                            >  2:214169/4‑251 (MQ=255)
                                                                                                                                                                                                          cATGCGATGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGCCACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCGGAACCCCCGGTTTTGTTc                            >  2:214176/1‑251 (MQ=255)
                                                                                                                                                                                                                attggaAACAGAGTGGTGTGTACCAACGCCAGCAGAGGAAGTGTTGTACCTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTGGGGCAAAGACAGAAATGCCCACGAAAACTCCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGACACTTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGAGCCGCGAGACTGTCGGGCACCCAGGTTTTGTTCAGaagg                      <  2:202652/247‑3 (MQ=255)
                                                                                                                                                                                                                aTGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCg                                              >  2:338034/1‑227 (MQ=255)
                                                                                                                                                                                                                aTGAGAAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCg                                              <  1:338034/227‑1 (MQ=255)
                                                                                                                                                                                                                 tcagaAACAGAGAGGTGGCAACCAGCGCCAGCAGAGGTAGTGTTTTACGTTTCATTATTTTTTCCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATTCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCGGACACCCCGGTTTTGTTCAGAatg                      >  2:366909/3‑250 (MQ=255)
                                                                                                                                                                                                                    gaAACAGAGTGGTGGCAACCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCGGACACCCCGGTTTTGTTCAGAatgatg                   >  2:52146/1‑250 (MQ=255)
                                                                                                                                                                                                                                 ggCAACCAACGCCAGCGGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCGGACACCCCGGTTTTGTTCAGAATGATGACGATGACAACAAt     >  2:264206/1‑251 (MQ=255)
                                                                                                                                                                                                                                   caaccaacGCCAGCAGAGGCAGTGTTTTACGTTTCATTCTTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACAGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCGGACACCACGGTTTTGTTCAGAATGATGACGATGACAACAATgg   >  2:339903/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaCCAACGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGGCAGGGCCGCGAGACTGTCGGACACCCCGGTTTTGTTCAGAATGATGACGATGACACAAATGGc  >  1:121648/1‑251 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
CACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCACGCCGACGAACGGTGTCAGCGCGACACTGGCAGCACCGGCTGCGCCACCTCCCGTCAGCAGTGCGGCAACAGTGGCAAGGGTCAGGGCCGCGAGACTGTCGGACACCCCGGTTTTGTTCAGAATGATGACGATGACAACAATGGC  >  NC_000913/266290‑266768

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