Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,148,049 C→G 8.3% G90A (GGC→GCC)  yijO ← AraC family putative transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,148,0490CG 8.3% 195.0 / 7.7 60G90A (GGC→GCC) yijOAraC family putative transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (25/30);  new base G (3/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

AGCGCGGAGGCATAGCGTTCGTCGATATAATCGCGAATAGCTTCAAATAGTGCCTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCAGATAGCTGACGTCGTGATACATGTTATTCCCTTAGGTATGTA  >  NC_000913/4147761‑4148337
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                 
aGCGCGGAGGCATAGCGTTCGTCGATATAATCGCGAATAGCTTCAAATAGTGCCTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaa                                                                                                                                                                                                                                                                                         <  1:78038/299‑2 (MQ=255)
   gcgGAGGCATAGCGTTCGTCGATATAATCGCGAATAGCTTCAAATAGTGCCTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCt                                                                                                                                                                                                                                                                                      >  1:518191/1‑299 (MQ=255)
       tctcATAGCGTTCGTCGATATAATCGCGAATAGCTTCAAATAGTGCCTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACAGCaactcggc                                                                                                                                                                                                                                                                                   <  1:167715/295‑8 (MQ=255)
         gCATAGCGTTCGTCGATATAATCGCGAATAGCTTCAAATAGTGCCTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAgg                                                                                                                                                                                                                                                                                 <  1:38319/298‑1 (MQ=255)
           aTAGCGTTCGTCGATATAATCGCGAATAGCTTCAAATAGTGCCTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGGCAAACa                                                                                                                                                                                                                                                                                             >  2:147188/1‑284 (MQ=255)
             aGCGTTCGTCGATATAATCGCGAATAGCTTCAAATAGTGCCTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGc                                                                                                                                                                                                                                                                                           <  2:495073/284‑1 (MQ=255)
                          ataATCGCGAATAGCTTCAAATAGTGCCTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGcc                                                                                                                                                                                                                                                               >  1:372013/1‑299 (MQ=255)
                                         ttCAAATAGTGCCTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAAtt                                                                                                                                                                                                                                                  <  1:337114/297‑1 (MQ=255)
                                         ttCAAACAGTGCCTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTcc                                                                                                                                                                                                                                                <  1:48283/299‑1 (MQ=255)
                                             aagaGTGCCTGGCTGCGTGGGTCGGTGTGGATTTGGCGTCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                                                                                                                                                                          <  2:116135/282‑1 (MQ=255)
                                             aaTAGTGCCTGGCTGCGTGAGGCGGGCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATtgt                                                                                                                                                                                                                                                          <  2:257352/285‑1 (MQ=255)
                                                 gggCCTGGCGGCGGGAGGCGGTCTGGATTTGGCTGCCAAGCAAAGCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGa                                                                                                                                                                                                                                                       <  2:439618/282‑1 (MQ=255)
                                                    ccTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcagc                                                                                                                                                                                                                                      <  1:585779/298‑1 (MQ=255)
                                                     cTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGAAGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCACTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTTCTGTTTGCCAAACAGCCCCCTAAAGGTATAAGCTGGGGATTGCCCTTGTGGAACAta                                                                                                                                                                                                                                                  >  2:509962/1‑284 (MQ=255)
                                                                 ggCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcagcagagatcggaaga                                                                                                                                                                                                                         >  1:52744/1‑287 (MQ=255)
                                                                   cGGTCTGGAGTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACcagca                                                                                                                                                                                                                                     <  2:45287/284‑1 (MQ=255)
                                                                       cTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGaa                                                                                                                                                                                                                                 >  2:242155/1‑284 (MQ=255)
                                                                         ggATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGt                                                                                                                                                                                                                              <  2:503425/285‑1 (MQ=255)
                                                                         ggATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCg                                                                                                                                                                                                                <  1:332663/299‑1 (MQ=255)
                                                                                            aTCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAg                                                                                                                                                                                             <  1:529545/299‑1 (MQ=255)
                                                                                             tCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCCGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAgc                                                                                                                                                                                            >  1:438694/1‑299 (MQ=255)
                                                                                                       ggCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATc                                                                                                                                                                                  <  1:434245/298‑1 (MQ=255)
                                                                                                               aggatggagatcaTAAGCCTTGGCTGTTGCTGCTCCTGCGGCTGCATTTGCAATTCATTGAGCGTTAGTAGGAGAAAAGAGCCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATGTTGAGATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGcc                                                                                                                                                                                        <  2:452348/272‑1 (MQ=255)
                                                                                                                             aaGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACaa                                                                                                                                                           >  1:269510/1‑299 (MQ=255)
                                                                                                                                  ttGCCGTTTGCTGCTCCTGCGGCTGCCTTTGCATTTCATTGAGCGTTTGTAGAAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGCGTTGCTGTTTGCCAACCAGCACGCTAAAGGTAGTAGCGGGGGCTTGCAATGGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCGGGGAAATAACGTCGCTCCGGCGCCGGTATAAACAAACTCACCTTCAGGCACAATTTca                                                                                                                                                      >  1:441758/1‑298 (MQ=255)
                                                                                                                                         ttGCTGCTCCTGCGGCTGCATTTGCATTTCATTGGGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGcac                                                                                                                                                             <  2:326238/285‑1 (MQ=255)
                                                                                                                                         ttGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGt                                                                                                                                                <  1:558237/298‑1 (MQ=255)
                                                                                                                                           gctgctCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGa                                                                                                                                             >  1:42613/1‑299 (MQ=255)
                                                                                                                                            ctgctCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGCAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCAACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCACGCACGGTATCACAAAACTCACCTTCCAGCACAAt                                                                                                                                                          >  2:436267/1‑285 (MQ=255)
                                                                                                                                            ctgctCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCCACAAACTCACCTTCCAGCACACTTTCCAGCCGTggca                                                                                                                                            >  1:373507/1‑297 (MQ=255)
                                                                                                                                                       gCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTg                                                                                                                                               >  2:334368/1‑285 (MQ=255)
                                                                                                                                                          gcatttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATa                                                                                                                               >  1:90561/1‑298 (MQ=255)
                                                                                                                                                            atttgcatttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGaaaa                                                                                                                                          <  2:77522/285‑1 (MQ=255)
                                                                                                                                                                   tttCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAATGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCg                                                                                                                     >  1:132452/1‑299 (MQ=255)
                                                                                                                                                                     tCATTGAGCGTTTGTAGCAGAAAAGAACCTAGGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTTATATTGTTTTCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCggg                                                                                                                   <  1:523596/299‑1 (MQ=255)
                                                                                                                                                                                       aGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTActg                                                                                                   <  1:381434/297‑1 (MQ=255)
                                                                                                                                                                                                           gCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCa                                                                             >  1:272878/1‑299 (MQ=255)
                                                                                                                                                                                                             cgcgTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAg                                                                            <  1:274524/298‑1 (MQ=255)
                                                                                                                                                                                                                           tgtTGCTTGTCCAGATTTTTATATTGGTTGCCATCCCCTTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAGCAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGt                                                              <  1:81476/298‑1 (MQ=255)
                                                                                                                                                                                                                           gtttGCTTCGCCAGATTGTGTTATTGTTTGCCATCCCAGTGCACGACGCTGTAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAg                                                                            <  2:237554/282‑1 (MQ=255)
                                                                                                                                                                                                                              tGCTTCGCCAGATTTTGATATTGTTTGCCAACCCATTGCACGACGCGGAAGCCGAGTGGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACGTCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGAt                                                          <  1:556992/299‑1 (MQ=255)
                                                                                                                                                                                                                               gCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATc                                                         >  1:73604/1‑299 (MQ=255)
                                                                                                                                                                                                                                   cGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGAACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACATTGCTCCAGCGACGGTCTCAACAAACTCACCTTCCAGCACAATTTCCAACCGTGGAAAATGGACTTGACACGCAAGACCGGGCACCGGGCCGTTACGGCTGGCAAAGTAAAACTGGCGAAGGGACAGCgg                                                                   >  2:17183/1‑285 (MQ=255)
                                                                                                                                                                                                                                       aGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGCACAGGCGGga                                                 >  1:454614/1‑299 (MQ=255)
                                                                                                                                                                                                                                            ttGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACACTTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCCAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGAt                                                          >  2:239947/1‑285 (MQ=255)
                                                                                                                                                                                                                                            ttGATATTGTTGTCAATCCCATTGCACGACGCTTGAGCAGGGTTTCTGGTGGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAGTGTGGAAAATTCCAGCCCCCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAAGCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGa                                                           <  2:170955/284‑1 (MQ=255)
                                                                                                                                                                                                                                               atatTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCAGAt                                         >  1:543557/1‑299 (MQ=255)
                                                                                                                                                                                                                                                tatTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCAGATa                                        >  1:411818/1‑299 (MQ=255)
                                                                                                                                                                                                                                                  ttGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACCGCGGGCCGTTGATCAgg                                                      >  2:513470/1‑283 (MQ=255)
                                                                                                                                                                                                                                                      ttGCCATCCCATTGCACGACGCGGAAGCCGAGTTGCGGTTTGCCAAACAGCACGCTAAAGGGAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAGCGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCAGATAGCTGa                                   <  1:55158/298‑1 (MQ=255)
                                                                                                                                                                                                                                                         ccATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCCACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCCAGATCGGGCACCGGGCCGTTACTGCTTGCACAGTAAATCTGCCCCAGGGACCGATGGCCGTTGATCAGCCgggggccc                                             >  2:314484/1‑280 (MQ=255)
                                                                                                                                                                                                                                                          cATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCa                                            >  2:194212/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                 ttGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCAGTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCAGATAGCt                                     <  2:64000/285‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCAGATAGCTGACGTCGTGATACATGTTAttag              >  1:544123/1‑297 (MQ=255)
                                                                                                                                                                                                                                                                               ctttagcgtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATATGGCGCAGGGACAGGGGGCCGTTGATCAGGCGGGGGAGCAGAAAGCTGAAGTAGTGATACATGTTCTTCCCTTAg         >  1:494636/9‑299 (MQ=255)
                                                                                                                                                                                                                                                                                 ccaaGTTGGTGTTTGCCAAACAGCACGCTAAAGGTAGGAGCGGGGGCTTGCCATTGTGGAAAAGCCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTAGCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCGAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCAGATAGCTGACGTCGTGATACAt                      <  2:452950/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 ccGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTTGAAAACTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGGTACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCAGATAGCTGACGTCGTGATACATGTTATTCCCTTAGGt       <  1:145106/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  cGAGTTGCTGGTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATGGTAGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAAGCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCAGATAGCTGACGTCGTGATACATGTTATTCCCTTAgg        <  1:20676/297‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     gTTGCTGTTGGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGACATTGTGGAAAAGTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCAGATAGCTGACGTCGTGATACATGTTAt                 <  2:474081/285‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      gtGCTGTTTGGCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCAGATAGCTGACGTCGTGATACATGTTATTCCCTTAGGTATGTa  >  1:460906/2‑299 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                 
AGCGCGGAGGCATAGCGTTCGTCGATATAATCGCGAATAGCTTCAAATAGTGCCTGGCTGCGTGAGGCGGTCTGGATTTGGCTGCCAAGCAAATCGCGGCAGTGGCTAAGCAGGCTGGCGACGATAAGCCTTGCCGTTTGCTGCTCCTGCGGCTGCATTTGCATTTCATTGAGCGTTTGTAGCAGAAAAGAACCTATGCGTGGGCCGCGTCGGGCGACGTGTTGCTTCGCCAGATTTTGATATTGTTTGCCATCCCATTGCACGACGCTGAAGCCGAGTTGCTGTTTGCCAAACAGCACGCTAAAGGTAGTAGCGGGGGCTTGCCATTGTGGAAAATTCCAGCCACCAGCAGGAACGTACAGCACATCGCCGGGAACTAACGTTGCTCCAGCGCCGGTATCAACAAACTCACCTTCCAGCACAATTTCCAGCCGTGGAAAATCGACCTGATACGCAAGATCGGGCACCGGGCCGTTACTGCTGGCAAAGTAAATCTGGCGCAGGGACAGCGGGCCGTTGATCAGGCGGGAGAGCAGATAGCTGACGTCGTGATACATGTTATTCCCTTAGGTATGTA  >  NC_000913/4147761‑4148337

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: