Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,965,174 G→A H153Y (CAC→TAC)  rhlB ← ATP‑dependent RNA helicase RhlB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,965,1740GA88.9% 111.9 / 3.0 45H153Y (CAC→TAC) rhlBATP‑dependent RNA helicase RhlB
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (17/23);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGA  >  NC_000913/3964928‑3965418
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
gCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTc                                                                                                                                                                                                                                                   <  1:233437/251‑1 (MQ=255)
          aaTACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCCTTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAACGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGGGGTTAGTGTAGTGCTGCTTGGCGt                                                                                                                                                                                                                                         >  1:218496/1‑251 (MQ=255)
          aaTACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCg                                                                                                                                                                                                                                          <  1:92004/250‑1 (MQ=255)
          aaTACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGt                                                                                                                                                                                                                                         >  1:291648/1‑251 (MQ=255)
                     ccGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGCGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGCTATCTTTAATAAAGCACAGATCGTCAATGCGATCGGCTTC‑GTCCAGATCCAACACCTGTATGGAACCGAGGGTAATGTAGTTCTGCTTGGCgatgtctcttat                                                                                                                                                                                                                               >  2:182722/1‑239 (MQ=255)
                       gtgttCTGTTCCTGGTCCTCTTCAATATAATCTGCATTATTGATATGCCCTACCGCCAGTTCAGGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATAACTTTAATAAAGCCCAGATCGTACATGCGATCGGCGTC‑GTCCAGTTACACCACCTGAATGGCACGGAGGTTAATGTAGTGCTGCTTGGCGT‑AGTCAATTAAACg                                                                                                                                                                                                                           <  2:230682/248‑1 (MQ=255)
                       gTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGCGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGACAGCGGGTATCTTTAATAAAGCCCAAATCGTACCTGCGGTCGGCCTC‑GTCCAGTACCACCACCTGAATGGAACCGAGGTTAATGTAGTTCTCCTTGGCCT‑AGGAAATTAAAcc                                                                                                                                                                                                                           >  2:253347/1‑250 (MQ=255)
                        ttttCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGAAATCTTTAATAAAGCCCAGATCGGACATGCGATCGGCTTC‑GTCCAGTACCACCAACTGAATGGCACCGAGGTTAATGTAGTTCTTCTTGGCGT‑AGACAATTAAActc                                                                                                                                                                                                                          >  2:209445/1‑249 (MQ=255)
                           tCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGTAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGAACGGCTTC‑GACCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGA‑CGGCAATTAAACGcccc                                                                                                                                                                                                                       >  2:133752/1‑251 (MQ=255)
                                        cACTTCAATAGATTCGGCATTGTGCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCa                                                                                                                                                                                                          <  2:107970/251‑1 (MQ=255)
                                                 atatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCaa                                                                                                                                                                                                 <  1:287288/251‑1 (MQ=255)
                                                 atatTCGGCAGTGTTCATCTGGTAGAACGCCAGTTTGTGTGCCCGGTACGAAAGCGTGGTGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACTGCCAGCGGATATCGTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTGCTGATTGGCGT‑AGTCAATTAACCTCGCCGTGGTGCCATTCAGAATGTCa                                                                                                                                                                                                  <  2:218496/250‑1 (MQ=255)
                                                                 aTCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACCGCATGTTGAGGCGCTGGATTGCAGTAGGCATACGGCGGAACAGCCAGCGGATATCTTTAAAAAAGCCCAGAACGTACATGCGATCGGCTGCAGTC‑AATACCACCACCTGAATGGTACCGAGGTTAATTAAGTTCTGCTTAG‑GTGAGTCCATTAAAACCCCCGGGGTGCCAAACAGAATTTAAACGCCCCTTTCCAacac                                                                                                                                                                                 >  2:339692/1‑251 (MQ=255)
                                                                   cTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAACCCAGATCGTACATGCGATCGGCTTC‑GTCCAGCACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGCGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACtt                                                                                                                                                                               >  1:312805/1‑251 (MQ=255)
                                                                     gCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCTCCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACttt                                                                                                                                                                              <  1:253347/250‑1 (MQ=255)
                                                                              ccAGTTCACGTACCAGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGCTCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGGGCCAATCAGAAAGTCAACGCCGCTTTCCAGCCCTTTCAGCTGTTTg                                                                                                                                                                    >  2:160866/1‑251 (MQ=255)
                                                                                   tcatgTACCCCGGACGAAAGCGTGTCGGAGACCAGAAGGTTGAGGGTCGGTTTTGCAGGCGGCATACGGCTGAACAGCCAGCGGAGAACTTTAATAAAGCCCAGATCGGACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGTCACCGAGGTTAATGTAGTTCTGCTTGTCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATTTCAACGCCGCTTTCCAGCTCTTTCAGCTTTTTGTCGt                                                                                                                                                                <  2:223965/246‑1 (MQ=255)
                                                                                         aCCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCAt                                                                                                                                                          <  1:209445/250‑1 (MQ=255)
                                                                                            cGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGTAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCTGCTTC‑GTCCATTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTTTCGTCGCCATcacc                                                                                                                                                      >  1:94845/1‑251 (MQ=255)
                                                                                                     aGCGTGGCGGAGAACTGCAGGTTGGGGCGCTGGTTTGCCGGCGGCATACGGCGGAACAGCCAGCTGATATATTTAATAAAGCCCAGATCGTACATGCGATCGGTTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGGGCCAATCAGAATGTCAACGCCGCTTTCCAGCTCTTTCGGCTGTTGGTCGTAGCCATCACCACCGTAAGc                                                                                                                                             <  2:312805/251‑1 (MQ=255)
                                                                                                      gCGTGTCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGcc                                                                                                                                            <  1:182259/251‑1 (MQ=255)
                                                                                                         tGGCGGAGAACAGCATGTTAAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGGTATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGGTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGCGCCAATCAGACTGACAACGCCGCTTTCCAGCACTTTCACCGGTTTGTCGTAGCCATCACCACCGTAACCCAgc                                                                                                                                         >  2:86127/1‑250 (MQ=255)
                                                                                                             ggAGAACAGCATGTTGAGGCGCTGGTTTGCAGGAGGCATACGGCGGAACAGCCCGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTAC‑GTCCAGTACCACCACCTGAATGGCACCGAGGGTCATGTAGTTCTTCTTGTCGT‑AGTCACTTAAACGCCCCGTGGTGCCAATCAGAATATCAACGCCGCTTTCCCGCACTTTTAGCTGTTTGTCGTAGCCATACCCAACGTAAGcctaaccc                                                                                                                                      >  2:322835/1‑244 (MQ=255)
                                                                                                                   cAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCa                                                                                                                               <  1:128442/251‑1 (MQ=255)
                                                                                                                                                         ggAACCGCCAGCGGATATCTTTAAAAAAGCCCACATCGTTCATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATgg                                                                                         <  1:51655/251‑1 (MQ=255)
                                                                                                                                                           aaCAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGCATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCCTCACCACCGTAAGCCAGCCCCAGCTTCAGGCAAGGAGCTTCCGCCAGCGGTTCTGCGTCGGCTTggc                                                                                        >  2:261793/1‑249 (MQ=255)
                                                                                                                                                                           cTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTc                                                                        >  1:227367/1‑250 (MQ=255)
                                                                                                                                                                               aataaAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCAcgcg                                                                   <  1:176275/251‑1 (MQ=255)
                                                                                                                                                                                    aGCCCAGACCGTACATGCGCTCGGCTGC‑GTACAGTTCCCCCACCGGAATGGCACCGAGGTGAATGTAGTTCTGCTTGGCGT‑AGTCAATTACACGCCCCGTGGTGCCAATCAGAATGGCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGt                                                              <  2:94845/251‑1 (MQ=255)
                                                                                                                                                                                     gCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCCCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACTGCAAGTTCACGCGTCGGTg                                                             <  1:279523/251‑1 (MQ=255)
                                                                                                                                                                                             cGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAAt                                                      <  1:29182/250‑1 (MQ=255)
                                                                                                                                                                                                               c‑GTCCAGTACCACCACCTGAATGGCTCCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGAtt                                    <  1:261793/250‑1 (MQ=255)
                                                                                                                                                                                                                  tCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCAc                                 >  1:139106/1‑251 (MQ=255)
                                                                                                                                                                                                                    cAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCt                               <  1:160866/251‑1 (MQ=255)
                                                                                                                                                                                                                     aGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGACAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCtt                              <  1:86127/251‑1 (MQ=255)
                                                                                                                                                                                                                       tACCACCACCTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCa                                                                                                             <  1:241915/170‑1 (MQ=255)
                                                                                                                                                                                                                       tACCACCACCTGAATGGCACCGAGGTTAATGTAGATCTGCTTGGCGT‑AGTCAATTAACCGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCCGGCCAGTAGCTTCCGCCa                                                                                                             >  2:241915/1‑170 (MQ=255)
                                                                                                                                                                                                                            ccaccTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTg                                                                                                    >  1:193357/1‑174 (MQ=255)
                                                                                                                                                                                                                             caccTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCAt                                                         <  2:109846/216‑1 (MQ=255)
                                                                                                                                                                                                                             caccTGAATGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCAt                                                         >  1:109846/1‑216 (MQ=255)
                                                                                                                                                                                                                                    aTGGCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCGATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGc               <  1:182722/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ggCACCGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCTCGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAg             <  1:293283/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cGAGGTTAATGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGACTGTCAACGCAGCTTTACAGCACTTGCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCCCGGCAAGTTCACGCGTCGATGCCTTAATTAAGGCACGCGGCGGATTCACCTTGCGATCGGCCATCGGAGGATgg        >  2:203237/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  aaTGTAGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATgagagaga  <  1:184303/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   aTGTAGTTCTGCAGGGCGT‑AGTCACTTAAACGCCCCGTGTTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGGCGTAGCCATCACCACCGTAAGCCCGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGAGTCGGCATGGATCTGCACGGCAAGTTCACGACTCGGCGCCATACTTAAGGCCCGCGGCTGAATCCCCTTGAGATCGGCGATCGCAGGATgggagag   >  2:266593/1‑248 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
GCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTC‑GTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGT‑AGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGA  >  NC_000913/3964928‑3965418

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