Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,367,944 G→A A40T (GCG→ACG)  yhcG → DUF1016 domain‑containing protein YhcG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,367,9440GA96.8% 191.6 / NA 65A40T (GCG→ACG) yhcGDUF1016 domain‑containing protein YhcG
Reads supporting (aligned to +/- strand):  ref base G (1/1);  major base A (27/33);  minor base C (1/1);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACAC  >  NC_000913/3367697‑3368181
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
gTTTCAGGAAATCTCTTCGCATATCTTGGTGAAAGCAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCACAGTCCATTATTCAACTGAGCTCCAGATAAAGGGACAGATAATGGAATCGATCTCAGAAGGAACCTCAGCAGGCTGCCAGCAAATCCACGACGTTATTAGTCATCTGGGCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTTATGACCACGa                                                                                                                                                                                                                                            <  2:524783/251‑1 (MQ=255)
     aGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTAc                                                                                                                                                                                                                                       >  1:337758/1‑251 (MQ=255)
        aaaGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAg                                                                                                                                                                                                                                    <  1:101764/251‑1 (MQ=255)
                  gCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAACCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTACTGACCACGACGTACCAGGAAATTGGcc                                                                                                                                                                                                                          >  1:537267/1‑251 (MQ=255)
                    atatCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCcg                                                                                                                                                                                                                         >  1:402594/1‑250 (MQ=255)
                      aacTTGGGGAAAGGAAATGTAATACGCTGTTCTTAATTGGTCAGTGCGCACCCCCCAGAGTATAGCTCATTCCATTATGCAACGTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCATGCTACCAGCAAAGCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCcgacg                                                                                                                                                                                                                      <  2:57054/249‑1 (MQ=255)
                             ggAAAGTAAAAGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAAGGGAATCGCTCTCAGAAGTAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTc                                                                                                                                                                                                               <  2:845431/251‑1 (MQ=255)
                                   tAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAAttt                                                                                                                                                                                                         >  1:228723/1‑251 (MQ=255)
                                     aaTGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTg                                                                                                                                                                                                        >  1:15976/1‑250 (MQ=255)
                                          aaTACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGCCGGAAACGGTACGTAGCGTTAACGCTTTAATGACCACGCCGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACa                                                                                                                                                                                                    >  2:548256/1‑249 (MQ=255)
                                             aCGCTGTTCTTATTTTTTCAGTGGCCAGCCCCCAAAGTATATCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGAGAATGGAATCGCTCGCTGAAGTAACCACAGCAGGCGACCAGCAAATCCACGTCGGTATTATTCATCTGGCCGATAGCGCCCGGACTGAAACGGTTCGTAGCGTTAACTCGTGAATGACCTCGACGTATCAGGAAATTGGCCGACGGATTGTCGAATTTGAACAAGGTg                                                                                                                                                                                               <  2:585211/251‑1 (MQ=255)
                                             aCGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCa                                                                                                                                                                                                                                     >  2:415358/1‑213 (MQ=255)
                                             aCGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCa                                                                                                                                                                                                                                     <  1:415358/213‑1 (MQ=255)
                                             aCGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGAGACCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGaaacagggg                                                                                                                                                                                               >  2:698811/1‑244 (MQ=255)
                                             aCGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACGTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTg                                                                                                                                                                                               <  2:15976/251‑1 (MQ=255)
                                               gCTGTTCTTACTTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGCCGGAATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGt                                                                                                                                                                                                                                         >  1:846425/1‑207 (MQ=255)
                                               gCTGTTCTTAATTGTTCATTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCACCATATAAAGGGACAGATAATGGAATCGCACTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCGCGACGGTAGTATTCATCGGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGt                                                                                                                                                                                                                                         <  2:846425/207‑1 (MQ=255)
                                               gCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGc                                                                                                                                                                                             >  2:275421/1‑251 (MQ=255)
                                                  gTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGACCAAGGTGGCg                                                                                                                                                                                            >  1:135170/1‑249 (MQ=255)
                                                  gTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCg                                                                                                                                                                                            >  2:61101/1‑249 (MQ=255)
                                                   ttCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAg                                                                                                                                                                                          <  2:769009/250‑1 (MQ=255)
                                                   ttCTTAATTGGTCAGTGCTCAACCCCCAGACTATAGATCATTCAATTATTCAACTTAGCTCCATATCAAGGGACAGCTAATGGAATCGCTCTCAGAAGGTACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTAGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGATTTGTCGAATTTGAACAAGGGGGCAAgg                                                                                                                                                                                         >  2:92943/1‑251 (MQ=255)
                                                       tAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGcc                                                                                                                                                                                       <  1:61101/249‑1 (MQ=255)
                                                               cAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGCCGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCg                                                                                                                                                                              <  1:151833/250‑1 (MQ=255)
                                                                      cACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGTACGGAAACGGTACGTAGTGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTg                                                                                                                                                                       <  1:109546/250‑1 (MQ=255)
                                                                                   ataGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCTGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACTACGTACCAGGAAATTGGCCTACGAATTGTCGAATTTGAACAAGGTGGCGAg                                                                                                                                                                                          <  1:805314/218‑1 (MQ=255)
                                                                                     aGCTCATTCCATTATTCAACGTAGCTCCATATAAAGAGACATATAATGGAATCGCTCTCAGAAGGAAATACAGCAGGCTACCAGCAAATCCACGACGGTATTAGTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATGTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGTTAATCAAGCg                                                                                                                                                       <  2:76044/251‑1 (MQ=255)
                                                                                                 tattCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGTCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAgg                                                                                                                                            <  1:333758/250‑1 (MQ=255)
                                                                                                 tattCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTGTGGTGCGCAGCTAATCAAGCGACTATCAAAggg                                                                                                                                           >  2:831736/1‑250 (MQ=255)
                                                                                                 tattCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGAGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACAAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAggg                                                                                                                                           >  1:85856/1‑250 (MQ=255)
                                                                                                   ttCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGACAACGGTACGTAGCGTTAACGCGTTAATGACCCCGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGGGTATGGAGCGCAGATAATCAAGCGACTATCAAAGGAtt                                                                                                                                         >  2:648050/1‑251 (MQ=255)
                                                                                                     ccaCTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCTCAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGTCTATCAAAGGATTTa                                                                                                                                       <  2:553965/249‑1 (MQ=255)
                                                                                                     cAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGAttt                                                                                                                                        >  1:250017/1‑250 (MQ=255)
                                                                                                     cAACCTAGCTCCATACAAAGGGACAGATAACGGAATCGCTCTCAGAAGGAACCACATCAGTCTACCAGCAAATACACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCAAGGGCTGCGTATGGTGCGCATCTAATCAAGCGACTAACAAAGGATTTa                                                                                                                                       <  2:553972/251‑1 (MQ=255)
                                                                                                               ccATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAgg                                                                                                                              <  1:21814/250‑1 (MQ=255)
                                                                                                                cATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTATCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGta                                                                                                                            <  1:126456/251‑1 (MQ=255)
                                                                                                                cATACAAAGGGACAGAGAATGGAATCGCTCACAGAAGGAACCACAGCATGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATGTTAACAAGGTGGCGAGGCAAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAgcta                                                                                                                            <  2:131013/251‑4 (MQ=255)
                                                                                                                     aaaGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGc                                                                                                                       >  1:682438/1‑251 (MQ=255)
                                                                                                                           aCAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGt                                                                                                                 <  1:172377/251‑1 (MQ=255)
                                                                                                                                    ggAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGGACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTAGAATTTGAACAAGGAGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATAAAGCAACTATCAAGGGATTTATGGCTAAGGGATAAGCGTGGGTTCCCTGCa                                                                                                         >  2:312978/1‑250 (MQ=255)
                                                                                                                                             ctcAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGACCAGGGCTGCGTATGGAGCGCCGCTAATCAAGCGACTATCAAAGGCTTTATGTCTAAGGTATAAGCGAGGGTTCTCTGCCAAAAACTTAc                                                                                               >  2:677484/1‑251 (MQ=255)
                                                                                                                                               ccgAAGGAACCACAGCAGGCTACCAGACAATCCACGACGGTTTTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAAAGCGTTAATGACCACGACGTACCAGGAAATTTGCCGACGAATTGGCGAAATTGAACAAGGTGGCGAGGCCAGTGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGGGGGTTCTTTGCAAAAAACTTAc                                                                                               <  2:766168/247‑1 (MQ=255)
                                                                                                                                                           cagcagGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCtt                                                                                 >  1:226476/1‑251 (MQ=255)
                                                                                                                                                              cagGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTAGAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCttttt                                                                              <  1:312978/251‑1 (MQ=255)
                                                                                                                                                               agGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCtttttt                                                                             <  1:258176/251‑1 (MQ=255)
                                                                                                                                                                    aCCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGACCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTAcc                                                                          <  2:402594/249‑1 (MQ=255)
                                                                                                                                                                      cAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACACGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTa                                                                                                                                                               >  2:277221/1‑162 (MQ=255)
                                                                                                                                                                      cAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTa                                                                                                                                                               <  1:277221/162‑1 (MQ=255)
                                                                                                                                                                        gCAAATCCACGACGGTATTATTCATCTGGTCGATAGAGACCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCtttt                                                                    <  1:698811/251‑1 (MQ=255)
                                                                                                                                                                          aaagccACGAAGGTATTATTCATTTGTTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAAAGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAAGTGTCGAATTTTAACAAGGGGGCGAGGCCAGGGTTGCGTATGGTGCGCATCTAAGCAAGCGACTATCAAAGTATTTATGTCTAAGTTATAAGCGTGGGTTCGCTGCAAAAATCTTACGCCAGATGAGTCTTTTTTCCCTCtcttta                                                                  <  2:537267/247‑2 (MQ=255)
                                                                                                                                                                                        attattCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGaaa                                                     <  1:695907/250‑1 (MQ=255)
                                                                                                                                                                                         ttattCATCTGGACGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGCGGCGAGGCCAGGGCTGCGCATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCACATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAAtt                                                   >  1:576246/1‑251 (MQ=255)
                                                                                                                                                                                           attCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTc                                                  >  1:290440/1‑250 (MQ=255)
                                                                                                                                                                                           attCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTc                                                  <  1:424198/250‑1 (MQ=255)
                                                                                                                                                                                            ttCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCAc                                                >  1:725788/1‑251 (MQ=255)
                                                                                                                                                                                              cATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCa                                              <  1:207668/251‑1 (MQ=255)
                                                                                                                                                                                                   ggTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAg                                             <  1:395172/247‑1 (MQ=255)
                                                                                                                                                                                                          aGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTgg                                  >  1:874283/1‑251 (MQ=255)
                                                                                                                                                                                                                              acttAGCGGTAACGCGTTAATGTCCACGACGTGCCATGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTAGCAAAGGATTTATGGCTAAGGTATAATCGTGGTTTCTCTGCAAAAAACTTACGCCAAATGAGGCTGTGTTACCTCTTTTTTCAACATGTTTAAAATCACCAGACAATGTCTGTCGAATTAACTCCACTGGGa               <  2:552370/247‑1 (MQ=255)
                                                                                                                                                                                                                                  aGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTTTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTcc          <  1:619163/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tatctcgTTAATGATCCCGACGTACCAGGAAAGTGGGCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGACTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTCCGCCAAATGAGGCTGTTTTTCCTCTTTTTTCAACATGGTGAAATTCCCCAGACAATGTCTGGCGAATTAACATCATTGCGAATTCCGCAg      <  2:412085/245‑1 (MQ=255)
                                                                                                                                                                                                                                        aaCGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGa     >  1:795670/1‑250 (MQ=255)
                                                                                                                                                                                                                                          cgcgTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGacac  >  1:93978/1‑251 (MQ=255)
                                                                                                                                                                                                                                               taaggACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGTCCAGGGCTGCGTATGGTTCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACAGGTTGAAAt                                                    <  2:488112/192‑1 (MQ=255)
                                                                                                                                                                                                                                               tAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAAt                                                    >  1:488112/1‑196 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
GTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACAC  >  NC_000913/3367697‑3368181

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