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

GGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGC  >  NC_000913/4616425‑4616913
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
ggTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCt                                                                                                                                                                                                                                                <  2:83153/251‑1 (MQ=255)
   ttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGGTCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTaa                                                                                                                                                                                                                                              <  1:281841/250‑1 (MQ=255)
   ttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGCCGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATACTTCTAAc                                                                                                                                                                                                                                             >  2:25321/1‑251 (MQ=255)
            tCACAGATGTTCTTCGCCACTGCCAGCAGCAGGTCAGCAGGGGTGATTTCAGGATCGTAGATAAAGTTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAAGGGGCGAATCAGACGGGCCGATATTGGCGGGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCATAAACGTTTTAGTCGAACAGCGATCTCGTCCTGTGGTAGAATTCTAACATACGGtt                                                                                                                                                                                                                                     <  2:420612/250‑1 (MQ=255)
             cacaGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGc                                                                                                                                                                                                                                   <  1:294274/251‑1 (MQ=255)
                aGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACGCTGCAGCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTAGTCGAACAGCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaac                                                                                                                                                                                                                                <  1:242240/251‑1 (MQ=255)
                           gCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATcc                                                                                                                                                                                                                       >  2:245082/1‑249 (MQ=255)
                                       aggagGTCAGCAGGGGTGATTTCAGGATCGGAGATAAAGGTCAGGTGCGGTGAAACCAGCTTCACCTCTGCATCTGCACTTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTa                                                                                                                                                                                                         <  2:157577/248‑1 (MQ=255)
                                       agcagGTCATCAGGGGTGATTTCAGGAGCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCGGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTa                                                                                                                                                                                                         <  2:293315/251‑1 (MQ=255)
                                              catcaGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACATGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGc                                                                                                                                                                                                  >  1:416894/1‑251 (MQ=255)
                                               atcaGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGc                                                                                                                                                                                                  >  1:49160/1‑250 (MQ=255)
                                                  aGGGGTGATTTCAGGATAGCAGATAAAGGTCAGGTTCGGTGAACCCTGCTTCAACTCTGCATCTGCACGTAAGATAGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCCGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCCTCCAGTTGCCCCAGGTAGACCGGCATCg                                                                                                                                                                                              >  2:413720/1‑251 (MQ=255)
                                                                gATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACg                                                                                                                                                                                                                            >  2:61894/1‑207 (MQ=255)
                                                                gATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACg                                                                                                                                                                                                                            <  1:61894/207‑1 (MQ=255)
                                                                                caggtttcgTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGa                                                                                                                                                                 >  1:115741/9‑250 (MQ=255)
                                                                                      cGGTGAAACCGGCTGCAACTCTGCATCTGCACGGAAGATCGGGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACAGACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTc                                                                                                                                                           <  1:148143/250‑1 (MQ=255)
                                                                                        gTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagc                                                                                                                                                        <  1:47916/251‑1 (MQ=255)
                                                                                         tGAAACCTGCTTCAACTCTGCAGCTGCGCGTAAGATAGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTACATAAAGGCTTCTGGCCAGTGT‑CTGGCGAGGGAACGGCAGAAACGTTTTATTCTCACATCGATCTCGTCTTTTGTTATAATGCTAACAGACGTTTGCAACAACGAATCAATTTTCCCCAGGTAGACCGGCATCGATGTGACCGACGGTCCGGGGGGGTAAAGGATGGTcagc                                                                                                                                                        >  2:290727/1‑250 (MQ=255)
                                                                                          gAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGAGGGTAATGGGAGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCAGGC‑AGGGATCAGTCGAGGTAACGCAAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGAAGACCGGCATAGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagc                                                                                                                                                        >  1:183690/1‑249 (MQ=255)
                                                                                             aCCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTcagcag                                                                                                                                                      >  2:122324/1‑248 (MQ=255)
                                                                                                      aaCTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGGGGGAAAGAATGGTCAGCAGAGAGAGTGCGtc                                                                                                                                          >  1:291004/1‑251 (MQ=255)
                                                                                                         tctgcatctgcaCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGc                                                                                                                                             >  1:7330/1‑245 (MQ=255)
                                                                                                         tctgcatctgcaCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGc                                                                                                                                             <  2:7330/245‑1 (MQ=255)
                                                                                                                   cccggaaGATCGCGCGGGTAAGGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAATGCGTCTGGC‑TGGGTTCTGTCTAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACTTACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACTGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTTATCAAGATCTTTc                                                                                                                             <  2:198609/246‑1 (MQ=255)
                                                                                                                   caCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAACGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTc                                                                                                                             >  2:387270/1‑251 (MQ=255)
                                                                                                                   caCGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGCCCGCTAATGGCGTGAATCAAGACGTCTGGC‑AGGGTTCTGTCGCGGTTCCGCCAGAAACGTTTTATTCGAACCTCGATCTCGTCTTGTGTTATAATTCTAACCTACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGAAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGCTCTTTc                                                                                                                             >  1:18730/1‑251 (MQ=255)
                                                                                                                                               aTCAGACGGGCAGCAATTGGCGAACATAGAGGCGTCTGGA‑ATGGTTCTGTCGCGGAAACGCCCGAAACGTTTTATGAGAACATCGATATCGTCTTGGGTTATACTTCTAACATAAGGTTGAAACAACGCATCAAGTTGCCACAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGGGGTAGAGAATGGTCAGCAGAGAGCGGGCGTCATCAGGATCTTTCGCGCCTTCCAGCGACAGCCATTCGGAAc                                                                                                 >  2:38539/1‑251 (MQ=255)
                                                                                                                                                      gggCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTTCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGaa                                                                                                  >  1:277532/1‑243 (MQ=255)
                                                                                                                                                      gggCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTAGTCGAACATCGACCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTCGCCCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGAGCTTTCGCGCCTTCAAGCTCCAGCCATTCGGaa                                                                                                  <  2:277532/243‑1 (MQ=255)
                                                                                                                                                       ggCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCACCCTCGCATCCTGTTGCCCCAGGTAGTCCGGCATCGATGTGACCGACGGTTCGTGGTGGTAAAGATTGGTCAGCAGAGATAGTGCGTCATCAAGCTCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGcc                                                                                         <  1:165250/251‑1 (MQ=255)
                                                                                                                                                       ggCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGcc                                                                                         >  2:325055/1‑251 (MQ=255)
                                                                                                                                                        gCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGcc                                                                                         >  1:149985/1‑250 (MQ=255)
                                                                                                                                                           gATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGa                                                                                      >  2:281134/1‑250 (MQ=255)
                                                                                                                                                            atatTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGGCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGaaa                                                                                    <  1:14762/251‑1 (MQ=255)
                                                                                                                                                                 ggCGTGCATAAAGGCGTCTGAC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTGTTGTTATAATTCTAACATACGGTTGCAACAACGCATCCTGTTTCCCCAGGTAGACCGGCATCGATGTGGCCGACGGTACGTGTTTGTAAAGAATGGTCAGCATAGAGAGTTCGTCATCAAGATCTTTCGCGCCTTCCAGCACCAGCCATTCGGAACCGTTCGCCAGAAAACgg                                                                                <  2:392654/250‑1 (MQ=255)
                                                                                                                                                                   cGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTg                                                                                                                                                                                                                  >  1:122569/1‑118 (MQ=255)
                                                                                                                                                                   cGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTg                                                                                                                                                                                                                  <  2:122569/118‑1 (MQ=255)
                                                                                                                                                                           aaGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATCCGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCggg                                                                     >  1:357881/1‑251 (MQ=255)
                                                                                                                                                                                    ggC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCGCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATa                                                            <  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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGGCGTGCATAAAGGCGTCTGGC‑AGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGGTCAGCAGAGAGAGTGCGTCATCAAGATCTTTCGCGCCTTCCAGCTCCAGCCATTCGGAACCGTTCGCCAGAAAACGGGCGTAATCGGGTAAGACATAGCGCGGTTTGTACGGCGCATGACCTTCAAACATATCGCAGATTACACCTTCATCCAGC  >  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: ‑