Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,486,205 C→A T29K (ACG→AAG)  crp → DNA‑binding transcriptional dual regulator CRP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,486,2050CA98.4% 210.6 / NA 64T29K (ACG→AAG) crpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (30/33);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

ATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATT  >  NC_000913/3485965‑3486440
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
aTACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGTCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTGATTCAc                                                                                                                                                                                                                                      <  2:210684/251‑1 (MQ=255)
  acatgtaTGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTTCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCa                                                                                                                                                                                                                                    <  2:495856/245‑1 (MQ=255)
  aCATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAccc                                                                                                                                                                                                                                    >  2:37473/1‑250 (MQ=255)
  aCATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAcc                                                                                                                                                                                                                                     >  2:124314/1‑250 (MQ=255)
         tGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGa                                                                                                                                                                                                                              >  1:218208/1‑250 (MQ=255)
                 aTGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCgg                                                                                                                                                                                                                      >  1:316110/1‑250 (MQ=255)
                   gtctGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGTTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGaa                                                                                                                                                                                                                    <  2:450730/247‑1 (MQ=255)
                     aggCAAAGGACGTCACATTTCCGGGCAGTACAGTTGATACCCCCTGCTCAGGTAGTGGGAAGCATATTCCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACg                                                                                                                                                                                                                 <  1:349279/249‑1 (MQ=255)
                         aaaGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGt                                                                                                                                                                                                             >  1:433050/1‑251 (MQ=255)
                          aaGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGtac                                                                                                                                                                                                           >  2:295194/1‑251 (MQ=255)
                          aaGGACGTCACATTACCGTGCAGTAAAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCTAGAGCAAGCTTNTTCACCAGGGTGAAAAAGCGGAAACGCTGt                                                                                                                                                                                                             >  2:360563/1‑250 (MQ=255)
                           aGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGAAAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCCTTGCCACATTCATAAATAACCATCCAAGAGAAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGtac                                                                                                                                                                                                           >  1:370397/1‑251 (MQ=255)
                             gACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGtact                                                                                                                                                                                                          >  1:452840/1‑250 (MQ=255)
                                   acaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGt                                                                                                                                                                                                   >  2:676331/1‑251 (MQ=255)
                                            ggcggtccggTTGATAGCCCCCTCCCAGGGGNGGGGAAGCATGTTTCGGAAATCCAGAGACAGGGGGGTTATCTTGCTCTGGAGAAGGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGctc                                                                                                                                                                                          <  2:365257/242‑1 (MQ=255)
                                              cAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTAATTGCCACATTCATAAGAACCCATCAAAGAGCAAGCTTATCCACCAGGGAGAAAAAGCGGAAACGCTGTACTAAATCGTTAAAGGCTCtg                                                                                                                                                                                        >  2:60588/1‑251 (MQ=255)
                                                   cAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGc                                                                                                                                                                                    >  1:270259/1‑250 (MQ=255)
                                                    aGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGAATCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGc                                                                                                                                                                                    >  2:591074/1‑249 (MQ=255)
                                                    aGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTCCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAg                                                                                                                                                                                  <  2:303222/251‑1 (MQ=255)
                                                                   cccAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGAc                                                                                                                                                                    <  1:154825/250‑1 (MQ=255)
                                                                       ggTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACGGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAaga                                                                                                                                                               <  2:467942/251‑1 (MQ=255)
                                                                        gTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAagag                                                                                                                                                              <  1:124314/251‑1 (MQ=255)
                                                                             gggAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTa                                                                                                                                                          >  1:549288/1‑250 (MQ=255)
                                                                             ggaaaGCAGATTGCGGAAACCCAGAGACAGCGGCTTTATCGTGCTCTGAAGAAAGCTTATAACAGTGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTaa                                                                                                                                                         <  2:279660/248‑1 (MQ=255)
                                                                              ggaaggagatTTCGGCAGTCCGGGGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGGGGATACCCGCGCGGGGTGCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGTCGAAGAGGGTaaa                                                                                                                                                        <  2:412136/243‑1 (MQ=255)
                                                                                   cATATTCCGGCAATCCAGAGCCAGCGGGGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGTATAACGGCTTATGGGGCTTGGAAAACCGCAAACAGCCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAAt                                                                                                                                                   <  1:345802/251‑1 (MQ=255)
                                                                                    atatTTCGGCAATCCAGAGACAGCGGCGTTATCGGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGAAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATg                                                                                                                                                  <  2:482712/251‑1 (MQ=255)
                                                                                               aTCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaaaaa                                                                                                                                                                                                                          >  2:719723/1‑168 (MQ=255)
                                                                                               aTCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaaaaa                                                                                                                                                                                                                          <  1:719723/168‑1 (MQ=255)
                                                                                                  caggaacgccggcgTTATTGGGCTCGGGAGANAGCTTATAACAGAGGATATCCGCGCATGTGGCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCt                                                                                                                                   <  2:319278/242‑1 (MQ=255)
                                                                                                   agagACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTg                                                                                                                                  <  1:518864/251‑1 (MQ=255)
                                                                                                      gaCAGCGGAGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGAATACCCGCGCATGGTGCTTGTCACTCCGCAAACATACCCGACTCTCGAATGGTTCTTGTCTCATTGCCCCATTCATAAGTACCCATCCAAGAGCAAGCGTATGCACCAGGGTGAAAAAGCGGAAACGCGGTTCTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCATAGACGAAGAGGGTAAAGGAATGATCCTCTCTTATCGGAAt                                                                                                                               <  2:491773/251‑1 (MQ=255)
                                                                                                       aCAGCGGCGCTATCTGGCTCTGGTGAAAGCTTATAACAGAGGATAAGCGCGCATGGGACTTGGCAAACCGCAAACAAAACCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGAACCCATCCAAGAGCAAACTTATTCCCAAGGGAGAAAAAGCGGAAACGCGTTACTACATCGTTAAAGGACCTGCTGCAGTGCTGATCAACGACGAAGAGGGGAAAAAAATGATCATCGCCTAGCTGAAt                                                                                                                               >  2:485490/1‑250 (MQ=255)
                                                                                                                  aTCTGGTTCTGGAGACAGCTTATAACAGAGGATAACCGCGCCTGGTGCTTGGCAAACAGCAATCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAACCGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAACGAAATGATCCTCTCCTATCTGAATCAGGGTGAtttt                                                                                                                   >  2:645425/1‑251 (MQ=255)
                                                                                                                  aTCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCAGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAtttt                                                                                                                   >  1:22878/1‑251 (MQ=255)
                                                                                                                      ggCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAtttt                                                                                                                   <  1:135683/247‑1 (MQ=255)
                                                                                                                            ggAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAAc                                                                                                         >  1:61287/1‑251 (MQ=255)
                                                                                                                             gagcaaGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACAGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACt                                                                                                        >  1:230396/5‑251 (MQ=255)
                                                                                                                              agaAAGCTTATATCAGAGGATAACCGCGCATGGGGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTg                                                                                                       <  2:514582/251‑1 (MQ=255)
                                                                                                                              agaAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCa                                                                                                                             >  2:518358/1‑229 (MQ=255)
                                                                                                                              agaAAGCTTATAACAGAGGATAACCGCGCAGGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCa                                                                                                                             <  1:518358/229‑1 (MQ=255)
                                                                                                                               aaaaGCTTAAAACAGAGGATAACCGCGCATGNTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTgg                                                                                                      <  2:316110/250‑1 (MQ=255)
                                                                                                                                aaaGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGGCCCGACTCGCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTTAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTggg                                                                                                     <  2:583152/251‑1 (MQ=255)
                                                                                                                                    cTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCCACTCTCGAATGGTTCTTGACTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGACAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAGTGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCtt                                                                                                 >  1:67783/1‑250 (MQ=255)
                                                                                                                                          aCAGGGGATAACCGCGATTGTTGCTGGGAAAACCGCAAAAATACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTAGCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGaa                                                                                           <  2:22878/251‑1 (MQ=255)
                                                                                                                                           cAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAag                                                                                          <  1:37473/251‑1 (MQ=255)
                                                                                                                                             gagGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAag                                                                                          >  1:296024/1‑249 (MQ=255)
                                                                                                                                                gATAACCGGGCATGGGGCTTGGCAAACCTAAAACAGACCCGACTCTCGAAGGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCGTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAGCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAAGCAGGGTGATTATATTGGCGAACTGGGCCTGTTTGAAGAGGGc                                                                                     <  2:137818/251‑1 (MQ=255)
                                                                                                                                                            tGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAgc                                                                        >  2:504796/1‑251 (MQ=255)
                                                                                                                                                                        atcagtaaacaaacCCGACTCTCGAATGGTTCTTGTCTCAGTGCCACATTCATAAGTACCCATCCAAGATCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACGGGGCCTGTTTGAAGAAGGCCAGGAACGTAGCGCATGGGTACGt                                                            <  1:350248/245‑1 (MQ=255)
                                                                                                                                                                        acccgcaatcacccccGACTCTCGAATGTGTCTGGTCTCATTGCCACATTCAAAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCTGAAACGCTGTACTACAGCGTTAAAGGCTCTGTGGCAGTGCTGAGCAAAGACGAAGAGGGTAAAGAAATGCTCCTCTCCTATCTGAAGCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACTTAGCGCATGGGTACGt                                                            <  2:43653/238‑1 (MQ=255)
                                                                                                                                                                        aaCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGt                                                            >  1:370196/1‑251 (MQ=255)
                                                                                                                                                                               aaCAGACCCGACTCTCGAGTGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGaaa                                                      >  2:606791/1‑250 (MQ=255)
                                                                                                                                                                                 cAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAAcc                                                   <  2:255252/251‑1 (MQ=255)
                                                                                                                                                                                                aaTGGTTCTTGTCTCATTGCCACATTCATACGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTCTACAACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAACACGCAGAGGGTAATGAAATGATCCTCTCCTATCTGCATCAGGGTGATTTTAGTGGCGAAGTGGACCTGCTTGATGAGGGCCGGGACCTACTCGCAAGGGTACGTGCGAAAACCGCCTGTGacgcggcg                                    >  2:76438/1‑244 (MQ=255)
                                                                                                                                                                                                         tttcccATTGCCACATTCATAAGTACCCACCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCGCTGGGGCAGTGCTGATAAAAGAAGAAGAGGTAAAAGAAATTATCCTCTCTTATCGGAATCAGGGTGATTTTATTGGCGATTGGGGCCTGTTTGAAGAGGGCCAGGAGGTTAGCGCGTGGGTACGGGCGAAAACCGCCTGTGAAGTGGCTGAAATTTc                           <  1:216270/245‑1 (MQ=255)
                                                                                                                                                                                                           tctcATTGCCACATTCATAAGTCCCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGt                         >  1:590430/1‑250 (MQ=255)
                                                                                                                                                                                                              cATTGCCACTTTCATAAGTCCCCACCCAAGAGCAAGCTTATTCCCCAGGGAGAAAAAGGGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATAAAAGACGAAGAGGGTAAAGAAATGATCCTTGCCTACCGGAATCAGGGTGTTTTTATTGGCGAACTGGGCTTGTTTTAAGAGGGCCAGGAACGTAGGGAAGGTGTTGGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTAc                       <  1:19866/249‑1 (MQ=255)
                                                                                                                                                                                                                  tccACATTCATAAGTACACATCCAAGAGCAAGCTTCTTCACCAGGGTGCAACAGCGGCAACGCTTTTCTACATCGTTTAAGGCCATGTGGCAGTGCAGATCAACGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTTAATCATGGTGATTTTATGTTCGAACTGGGCCTTTTTGACGAGGGCCCGGAACGTAGCCCCTGGGACCGTGCGCACACCGCATGTGAATGTGCTGAAATTTCTTACaaaaa                  >  2:256952/2‑250 (MQ=255)
                                                                                                                                                                                                                   ccACATTCATAAGTACCCATCCAAGAGCAAGTTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCTTTAAAGGCTCTGTGGCAGTGCTGATAAAAGACGAAGAGTGTAAAGAAAGGATCCTCTCCTATCGGAATCGGGGTGATTTTATTGGCGACCTGGGCCTTTTTGAAGAGGCCCAGAAACGTAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAAt                <  1:121659/251‑1 (MQ=255)
                                                                                                                                                                                                                       aTTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCtg                                              <  1:366533/218‑1 (MQ=255)
                                                                                                                                                                                                                       aTTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGAGATTTTATTGGCGAAATGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCtg                                              >  2:366533/1‑218 (MQ=255)
                                                                                                                                                                                                                               gTACCCATCCAAGAGCAAGCTTATTCACCAGTGTTAAAAAGCGGAAACGCTGTCCTACATCGTTAAAGGCTCTGGGGCAGGGCTGATAAAAGACGAAGAGGGTAAAGAAATGTTCCTTCCTTATCTGAACAAGGGTGATTTTGTTGGGAAACTGGGCTTTTTTGAAGAGGCCCAGGAAGGTAGCGCAGGGTTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGa    <  1:37375/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aCCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGAtt  >  1:451011/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
ATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATT  >  NC_000913/3485965‑3486440

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