Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,850 A→G Y370Y (TAT→TAC)  yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009136,8500AG98.0% 147.2 / NA 50Y370Y (TAT→TAC) yaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (32/17);  minor base T (0/1);  total (32/18)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGC  >  NC_000913/6608‑7096
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
gggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTa                                                                                                                                                                                                                                                 <  2:236954/251‑1 (MQ=255)
 ggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATgg                                                                                                                                                                                                                                                  >  2:155584/1‑248 (MQ=255)
  gcgcACGCCGAGTTTACGCTGGCGCAGATAATCACTGGCAATGGTATGAACCACAGGCGAGATCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGTCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTTCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGt                                                                                                                                                                                                                                               <  1:123991/251‑1 (MQ=255)
         cctaGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGAGATCTGCCAGTTGCCACATCATCGGAAGACTTAGCAAGGTTCCGCCGATGACCGTTGCGAAGGTGCAGAGCCGCAATCACCAGATCTCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAAc                                                                                                                                                                                                                                        <  1:402158/248‑1 (MQ=255)
           gAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGTAGATCCGGAAACACCAGATCGCTTTAGGGTTGTTCAGGGGGAAAAAGAAGAGATTGTTTTC‑GGCGTAAAAGTAGTTGGAAACg                                                                                                                                                                                                                                       >  2:290581/1‑250 (MQ=255)
                  cGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAACTTTGTAATAGCCATGCAGGCCATTATGATATCTGCCCGTTGCCACATCAGCGGAAGAATTAGCAAGGGGCCGCCGAAGAACGTTGCGAAGGTGCAAATCCGCAAACACCAGATCGCGTTCGGGATGTTCAGGCGTAAAAAGAAGAGAATGATTTC‑GGCGTTAATGTAGTTGGCAAAGATGGAGc                                                                                                                                                                                                                                >  2:190573/1‑250 (MQ=255)
                     gggCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGATAGCAGTAAAATGGCGGTCAAATTGGTAAGAGCCATGCAGGCCATTATGATATCTGCCATGTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGaa                                                                                                                                                                                                                            <  2:205916/250‑1 (MQ=255)
                        cGTAGCTAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGGAATAGCCATGCACGCCATTATGATATCTGCCCGTTTCCCCATCAGCGGAATACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGGGCAGATCCGCAAACACCAGATCGGTTTAGGGATGTTCAGGCGGAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGGTGGCAACGATGGAGCTGAAgg                                                                                                                                                                                                                          >  1:42441/1‑250 (MQ=255)
                            gATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCCGTTGCCCCATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGAGAAGGGGCAGATCCGCAAACACCCGATCGCTTTAGGGGTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGGTGGCAACGAAGGAGCTGAAGGCAAAc                                                                                                                                                                                                                     >  2:345721/1‑251 (MQ=255)
                                    cTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGGAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGGTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCGGAAGGCAAACAGAATAAc                                                                                                                                                                                                             >  1:337547/1‑251 (MQ=255)
                                      ggCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAAcc                                                                                                                                                                                                            <  2:457917/250‑1 (MQ=255)
                                                        aGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTc                                                                                                                                                                                           >  1:230600/1‑249 (MQ=255)
                                                            gagaGCAGTAAAATGGCGGTCAAATTGGTAAGAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCAcc                                                                                                                                                                                     <  2:75257/251‑1 (MQ=255)
                                                                     aaaaTGGCGGTCAAATTGGTAATAGCCCTGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCCGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCa                                                                                                                                                                                                                                          >  2:257558/1‑189 (MQ=255)
                                                                     aaaaTGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCa                                                                                                                                                                                                                                          <  1:257558/189‑1 (MQ=255)
                                                                      aaaTGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAAccc                                                                                                                                                                           >  2:542506/1‑251 (MQ=255)
                                                                       aaTGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGATGGTAACAAACTCAGCACCCCAGGAACCCa                                                                                                                                                                          <  2:308851/251‑1 (MQ=255)
                                                                          ggCGGTAAAATTGGTAATAGCCATGCAGGCCCTTATGATATCGGGCAGTTAACACATCAGCGGAAGTATTCGCCTGATGACGCCGATGAACGTTGCGAAGGTGCAGATCCGCAAACACAAGATAGCTTTAGGGTTGTTCTGGCGTACACAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAATAGAATAACCACGAGGCTGACAAACTCAGCACCCCAGGAACCCAtt                                                                                                                                                                        >  2:59482/1‑250 (MQ=255)
                                                                           gCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGACGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑a                                                                                                                                                                      >  1:272819/1‑250 (MQ=255)
                                                                             ggTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCACAGGAACCCATT‑aac                                                                                                                                                                    >  1:51876/1‑248 (MQ=255)
                                                                             ggTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGGTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAAc                                                                                                                                                                                                 >  1:426211/1‑222 (MQ=255)
                                                                                            aGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGTGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTc                                                                                                                                                     <  1:550022/250‑1 (MQ=255)
                                                                                            aGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCAACCCAGGAACCCATT‑AACACCCGCATCGCCTTCt                                                                                                                                                    >  1:19329/1‑251 (MQ=255)
                                                                                            aGCCATGCAGGCCATTATGATAGCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCAACGCCTTc                                                                                                                                                     <  2:2197/250‑1 (MQ=255)
                                                                                                   cAGGCCATTACGATCTCTGCCGGTTGCCACATCAGCGGAATACTTAGCAAGGTGACGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGCTCGCTTTAGGGATGTTCCGGCGTACAAAGCAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGAAAACAGAATAACCACGAGGGGAACAAACTCAGCACCCCAGGAACCCATT‑AACACCGGCATCGCCTTCTGGATa                                                                                                                                               >  2:304487/1‑249 (MQ=255)
                                                                                                                 atCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCAGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGAAAAAAACTAAGCACCCCAGGA‑CCCATTTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTcc                                                                                                                                >  1:572259/1‑250 (MQ=255)
                                                                                                                     gCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGAGCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCgg                                                                                                                           <  1:399113/251‑1 (MQ=255)
                                                                                                                       cAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCa                                                                                                                         <  1:155584/251‑1 (MQ=255)
                                                                                                                        aGTTGCCACATCAGCGGACGCCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCAt                                                                                                                        >  1:472022/1‑251 (MQ=255)
                                                                                                                        aGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCAt                                                                                                                        >  2:394866/1‑251 (MQ=255)
                                                                                                                                 aTCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTtgt                                                                                                               >  1:136287/1‑251 (MQ=255)
                                                                                                                                            acTTAGAAAGGCGCCGCCGCTGACCGTAGCGCAGGGGCAGAACCGCAAAGACCCGCTCGCTTCAGGGTTGTTCAGGCATAAAAAGAAGAGATTGCTTTA‑GGCGTAAATGAAGTTGGCAACGATGGAGATGAAGGCAAACAGAATAACCACGAGGGAAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGATTGTGCCGTTAACCg                                                                                                    >  2:410081/2‑251 (MQ=255)
                                                                                                                                                  cAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGt                                                                                               >  2:219594/1‑250 (MQ=255)
                                                                                                                                                           gccgATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGc                                                                                      >  1:306763/1‑250 (MQ=255)
                                                                                                                                                                  aCCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGAATGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCt                                                                              >  1:169411/1‑251 (MQ=255)
                                                                                                                                                                     gTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGc                                                                           >  1:224047/1‑251 (MQ=255)
                                                                                                                                                                        gCGAAGGGGCAGACCCGCAACCAACAGCGCGTTTTGGGGTTTTTCAGGGGTAAAAAGAAGAGATGGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCAAAAGGCAAACAAAATAACCACGAGGGTAACAAACTCACCCCCCCATGACCCCATT‑ACCCCCCGCATCGCCTTCGGAAAAAGCTGAATCCCTTCCAGCGGCAGGTAGGTTGTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGt                                                                        <  1:582244/251‑1 (MQ=255)
                                                                                                                                                                           aaGGTGCAGATCCGAAAACACAGGTTCTTTTTGGGGTTGTCCAGGCGAAAAAAGAAGAGTTTTTTTGC‑GGCTTAAATGTAGTTGGCAACGATGGAGCGGAGGGAAAACAGAATAACCCCGGGGGTACCAAACCCAGCCCCCCAGGACCCCTTT‑ACCCCCCGCAGCGCCTTCTGGATAAGCTGAATCCCTCCCAGCGGCATGTAGGTTGGGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACa                                                                     <  1:394866/251‑1 (MQ=255)
                                                                                                                                                                                        gCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAACAGAAGAGCTTGTTTTC‑GGCGTCAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGACTAACCACGAGGGTAACACACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGgtg                                                        >  2:394948/1‑251 (MQ=255)
                                                                                                                                                                                           aaCACCAGACCGTTTTAGGGTGGTTCGGGCAAAAAAAAAAGAAATTTTTTCC‑GGCGTAAATGAAGTGGGCAACGATGGAGCGAAAGGCAAAAAAAATACCCACGAGGGTAACAAACCCACCCCCCCAGGACCCCATT‑ACCCCCCGCATCCCTTTCTGGAAAAGCAGAATCCCTCCCAGCGGCATGTGGGTTGGGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCg                                                     <  1:468409/251‑1 (MQ=255)
                                                                                                                                                                                              cccAGACCGTTTTAGGGTTGTTCGGGCAAAAAAAAAAGAGAT‑GTTTTCCGGCGAAAATGTAGTGGGAAACGAGGGAGCGGAAGGAAAACAAAATAACCACGGGGATAACAACCTCAGCACCCCAGAACCCCATT‑ACCCCCCGACCCGCCTCCGGGAAAACCTAAATCCCTCCCGGCGGCATGTGGGTTGTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATa                                                  <  2:230600/250‑1 (MQ=255)
                                                                                                                                                                                                  gATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATaaaa                                               >  1:519143/1‑250 (MQ=255)
                                                                                                                                                                                                           aGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAatc                                      >  1:85980/1‑250 (MQ=255)
                                                                                                                                                                                                             ggTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAatcatc                                   >  1:114360/1‑251 (MQ=255)
                                                                                                                                                                                                                   tCAGGCGAAAAAAGAAGAGATTGTTTGC‑GGCGAAATTGTAGTTGGCAACGAGGGAGCGGAAGGCAAACAGAATAACCACGAGGGTAACAACCTCAGCACCCCAGGAACCCTTT‑AACCCCCGCATCTCCTTCGTGATAAGCTGAATCCCTCCCAGCGGCATGGGGGTTGTGCCTTCACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACTAGGGTGTCGATAAAAATGCCAATCATCTGGAc                              <  1:446226/250‑1 (MQ=255)
                                                                                                                                                                                                                            aaaaGAAGCGATAGTTTTC‑GACGTAAATGTAGTTGCCAACGATGGAGCTGAAGGCACACAGAATAACCACGAGGGTAACAAACTCCGCACCCAAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGCATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACACGCCAGTACTATCAGCCTGGCGCTTGCCGTACAGATGACCAGGGGGCCGATAAAACTGCCAATCATCTGGACAATCCCTTgc                    >  2:319605/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ttGTTTTC‑GGCGTAAATGTAGTTGGCAACGCTGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATg         >  2:259555/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ttGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCCGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGAGTCGATAACAATGCTAATCATCTGGACAATCCCTTTCGCTGCCGGATg         >  2:221938/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ttttttC‑GGCGTAAATGTAGTGGGCAAGGAGGGAGCTGAAGGCAAAAAGAATAACCACGAGGGAAACAAACTCCGCCCCCCAGGAACCCATT‑AACCCCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCGGCGGAATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATg         <  1:484517/248‑1 (MQ=255)
                                                                                                                                                                                                                                              c‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGCCCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGc  >  1:56730/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
GGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGC  >  NC_000913/6608‑7096

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