Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,616,669 G→T L133I (CTA→ATA)  yjjI ← DUF3029 domain‑containing protein YjjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,616,6690GT96.3% 185.8 / ‑6.0 54L133I (CTA→ATA) yjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (30/22);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGC  >  NC_000913/4616425‑4616913
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
ggTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCt                                                                                                                                                                                                                                                <  2:83153/251‑1 (MQ=255)
   ttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGGTCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTaa                                                                                                                                                                                                                                              <  1:281841/250‑1 (MQ=255)
   ttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGCCGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATACTTCTAAc                                                                                                                                                                                                                                             >  2:25321/1‑251 (MQ=255)
            tCACAGATGTTCTTCGCCACTGCCAGCAGCAGGTCAGCAGGGGTGATTTCAGGATCGTAGATAAAGTTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAAGGGGCGAATCAGACGGGCCGATATTGGCGGGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCATAAACGTTTTAGTCGAACAGCGATCTCGTCCTGTGGTAGAATTCTAACATACGGtt                                                                                                                                                                                                                                     <  2:420612/250‑1 (MQ=255)
             cacaGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGc                                                                                                                                                                                                                                   <  1:294274/251‑1 (MQ=255)
                aGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACGCTGCAGCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTAGTCGAACAGCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaac                                                                                                                                                                                                                                <  1:242240/251‑1 (MQ=255)
                           gCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATcc                                                                                                                                                                                                                       >  2:245082/1‑249 (MQ=255)
                                       aggagGTCAGCAGGGGTGATTTCAGGATCGGAGATAAAGGTCAGGTGCGGTGAAACCAGCTTCACCTCTGCATCTGCACTTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTa                                                                                                                                                                                                         <  2:157577/248‑1 (MQ=255)
                                       agcagGTCATCAGGGGTGATTTCAGGAGCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCGGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTa                                                                                                                                                                                                         <  2:293315/251‑1 (MQ=255)
                                              catcaGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACATGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGc                                                                                                                                                                                                  >  1:416894/1‑251 (MQ=255)
                                               atcaGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGc                                                                                                                                                                                                  >  1:49160/1‑250 (MQ=255)
                                                  aGGGGTGATTTCAGGATAGCAGATAAAGGTCAGGTTCGGTGAACCCTGCTTCAACTCTGCATCTGCACGTAAGATAGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCCGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCCTCCAGTTGCCCCAGGTAGACCGGCATCg                                                                                                                                                                                              >  2:413720/1‑251 (MQ=255)
                                                                gATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACg                                                                                                                                                                                                                            >  2:61894/1‑207 (MQ=255)
                                                                gATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACg                                                                                                                                                                                                                            <  1:61894/207‑1 (MQ=255)
                                                                                caggtttcgTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGa                                                                                                                                                                 >  1:115741/9‑250 (MQ=255)
                                                                                      cGGTGAAACCGGCTGCAACTCTGCATCTGCACGGAAGATCGGGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACAGACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTc                                                                                                                                                           <  1:148143/250‑1 (MQ=255)
                                                                                        gTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagc                                                                                                                                                        <  1:47916/251‑1 (MQ=255)
                                                                                         tGAAACCTGCTTCAACTCTGCAGCTGCGCGTAAGATAGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTACATAAAGGCTTCTGGCCAGTGTCTGGCGAGGGAACGGCAGAAACGTTTTATTCTCACATCGATCTCGTCTTTTGTTATAATGCTAACAGACGTTTGCAACAACGAATCAATTTTCCCCAGGTAGACCGGCATCGATGTGACCGACGGTCCGGGGGGGTAAAGGATGGTcagc                                                                                                                                                        >  2:290727/1‑250 (MQ=255)
                                                                                          gAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGAGGGTAATGGGAGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCAGGCAGGGATCAGTCGAGGTAACGCAAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGAAGACCGGCATAGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagc                                                                                                                                                        >  1:183690/1‑249 (MQ=255)
                                                                                             aCCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagcag                                                                                                                                                      >  2:122324/1‑248 (MQ=255)
                                                                                                      aaCTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGGGGGAAAGAATGGTCAGCAGAGAGAGTGCGtc                                                                                                                                          >  1:291004/1‑251 (MQ=255)
                                                                                                         tctgcatctgcaCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGc                                                                                                                                             >  1:7330/1‑245 (MQ=255)
                                                                                                         tctgcatctgcaCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGc                                                                                                                                             <  2:7330/245‑1 (MQ=255)
                                                                                                                   cccggaaGATCGCGCGGGTAAGGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAATGCGTCTGGCTGGGTTCTGTCTAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACTTACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACTGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTTATCAAGATCTTTc                                                                                                                             <  2:198609/246‑1 (MQ=255)
                                                                                                                   caCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAACGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTc                                                                                                                             >  2:387270/1‑251 (MQ=255)
                                                                                                                   caCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGCCCGCTAATGGCGTGAATCAAGACGTCTGGCAGGGTTCTGTCGCGGTTCCGCCAGAAACGTTTTATTCGAACCTCGATCTCGTCTTGTGTTATAATTCTAACCTACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGAAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGCTCTTTc                                                                                                                             >  1:18730/1‑251 (MQ=255)
                                                                                                                                               aTCAGACGGGCAGCAATTGGCGAACATAGAGGCGTCTGGAATGGTTCTGTCGCGGAAACGCCCGAAACGTTTTATGAGAACATCGATATCGTCTTGGGTTATACTTCTAACATAAGGTTGAAACAACGCATCAAGTTGCCACAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGGGGTAGAGAATGGTCAGCAGAGAGCGGGCGTCATCAGGATCTTTCGCGCCTTCCAGCGACAGCCATTCGGAAc                                                                                                 >  2:38539/1‑251 (MQ=255)
                                                                                                                                                      gggCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTTCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGaa                                                                                                  >  1:277532/1‑243 (MQ=255)
                                                                                                                                                      gggCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTAGTCGAACATCGACCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTCGCCCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGAGCTTTCGCGCCTTCAAGCTCCAGCCATTCGGaa                                                                                                  <  2:277532/243‑1 (MQ=255)
                                                                                                                                                       ggCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCACCCTCGCATCCTGTTGCCCCAGGTAGTCCGGCATCGATGTGACCGACGGTTCGTGGTGGTAAAGATTGGTCAGCAGAGATAGTGCGTCATCAAGCTCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGcc                                                                                         <  1:165250/251‑1 (MQ=255)
                                                                                                                                                       ggCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGcc                                                                                         >  2:325055/1‑251 (MQ=255)
                                                                                                                                                        gCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGcc                                                                                         >  1:149985/1‑250 (MQ=255)
                                                                                                                                                           gATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGa                                                                                      >  2:281134/1‑250 (MQ=255)
                                                                                                                                                            atatTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGGCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGaaa                                                                                    <  1:14762/251‑1 (MQ=255)
                                                                                                                                                                 ggCGTGCATAAAGGCGTCTGACAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTGTTGTTATAATTCTAACATACGGTTGCAACAACGCATCCTGTTTCCCCAGGTAGACCGGCATCGATGTGGCCGACGGTACGTGTTTGTAAAGAATGGTCAGCATAGAGAGTTCGTCATCAAGATCTTTCGCGCCTTCCAGCACCAGCCATTCGGAACCGTTCGCCAGAAAACgg                                                                                <  2:392654/250‑1 (MQ=255)
                                                                                                                                                                   cGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTg                                                                                                                                                                                                                  >  1:122569/1‑118 (MQ=255)
                                                                                                                                                                   cGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTg                                                                                                                                                                                                                  <  2:122569/118‑1 (MQ=255)
                                                                                                                                                                           aaGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATCCGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCggg                                                                     >  1:357881/1‑251 (MQ=255)
                                                                                                                                                                                    ggCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCGCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATa                                                            <  2:149985/251‑1 (MQ=255)
                                                                                                                                                                                             tCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCCGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGt                                                     <  1:256942/250‑1 (MQ=255)
                                                                                                                                                                                                      gTAACGACAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCACCAGGTAGACCGGCATCGATGTGCCCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGgcg                                           >  1:333293/1‑251 (MQ=255)
                                                                                                                                                                                                        aaCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGgcgc                                          >  2:101436/1‑250 (MQ=255)
                                                                                                                                                                                                         aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCAt                                        >  2:75966/1‑251 (MQ=255)
                                                                                                                                                                                                               gAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCtt                                  <  1:122324/251‑1 (MQ=255)
                                                                                                                                                                                                                 aaCGTTTTATTCGAACATCGCTCTCGTCTTGTGTTATAATTCTAACATACGTTTGAAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGAGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTc                                 >  2:234329/1‑250 (MQ=255)
                                                                                                                                                                                                                       ttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACat                           >  2:11169/1‑250 (MQ=255)
                                                                                                                                                                                                                                 aTCGATCTCGTCTTTTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGTCCTGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCAGCAAGATCTTTCGCGCCTTCCAGCTCCATCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGGTTGTACGGCGCATGACCTTCAAACATATCGCAGATTa                <  2:94185/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCGAACATATCGCAGAtt                 >  2:312876/1‑250 (MQ=255)
                                                                                                                                                                                                                                 aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTAAGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGt                                                 >  1:201489/1‑218 (MQ=255)
                                                                                                                                                                                                                                 aTCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCAGCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCGTGCGCGCCTTCCATCTCCAGCCATTGGGAACCGTTCGGCAGAAAACTGGCGTAAACGGGTAAGTCATAGCGCGGTTTGt                                                 <  2:201489/218‑1 (MQ=255)
                                                                                                                                                                                                                                 aTCGATCTAGTCTTGTGTTATACTTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACAGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGACAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGAtt                 >  2:397939/1‑250 (MQ=255)
                                                                                                                                                                                                                                         cGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCtt          <  2:297619/249‑1 (MQ=255)
                                                                                                                                                                                                                                            cTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATcc     >  2:411766/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gtgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAgc  <  1:411766/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGC  >  NC_000913/4616425‑4616913

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