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% 221.9 / NA 63T29K (ACG→AAG) crpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (30/32);  minor base G (0/1);  total (30/33)
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

ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAA  >  NC_000913/3485959‑3486447
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
aCAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCGGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCt                                                                                                                                                                                                                                                   <  1:38578/250‑1 (MQ=255)
             gATGTCCTGAATGGATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATCGCCCCTTCCCAGGAAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTAAAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGACTGGTTCTTGTCTCATTGCCACCTTCATAAGAACCCATCCAAGAGCAAGCTTATTCACCAGGGTg                                                                                                                                                                                                                                     >  2:415228/1‑251 (MQ=255)
              aTGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCTGGGTg                                                                                                                                                                                                                                     >  1:432541/1‑250 (MQ=255)
                gTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCCAGAGCAAGCTTATTCACCAGGGTGaa                                                                                                                                                                                                                                   >  2:362391/1‑250 (MQ=255)
                  aCTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCGGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaaaaa                                                                                                                                                                                                                                <  2:560536/251‑1 (MQ=255)
                     gcagggaTGCAAAGGACGTCACATTACCTTGTAGTTCAGTTTATAGCCCTTTCCCAGTTAGCGGGAAGCATATTTCGGCAATCAAGAGACAGCGGCGTGATCTGGCTATGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTCCTTGGCAAACCGCAAACAGGCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCg                                                                                                                                                                                                                             <  2:354898/245‑1 (MQ=255)
                              cAAAGGCCGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCAAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCACACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCAAATTCATAAGTAACCATCCAAGAGCAAGCTTATTCACCAGGGCGAAAAAGCGGACACGCTg                                                                                                                                                                                                                    >  2:110826/1‑251 (MQ=255)
                               aaaGGACGTCACAATACCGAACAGTACAGATGATAACCCTTTCCCAGGTCGCGGGAAGCATATATCGGTACTCCCGAGACAGCGGCGTTATGTGGCTCTGCAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCATCGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATACATAAGTACCCATACAAGCTCAAGCTTATTCACCAGGGTGAACAAGCGGAAACGCTg                                                                                                                                                                                                                    >  2:226837/1‑250 (MQ=255)
                                      gTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAt                                                                                                                                                                                                            >  1:364837/1‑251 (MQ=255)
                                      gTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAt                                                                                                                                                                                                            >  1:371885/1‑251 (MQ=255)
                                      gTCACATTACCGTGCAGTACAGTTGATAGCCCCTGCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGTCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAt                                                                                                                                                                                                            <  2:509537/251‑1 (MQ=255)
                                        cacaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTATCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAta                                                                                                                                                                                                           <  2:164210/250‑2 (MQ=255)
                                            ttACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTa                                                                                                                                                                                                       >  2:370932/1‑250 (MQ=255)
                                                cggacAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATGTCGGCAAGCCAGAGACAGCGTCGTTATCTGGCGCGGGAAAAAGCTTATAACAGAGGATAACAGCGCATGGTGTTTGGCAAACCGCAATCAGCCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACAGTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACTTCTTTAAAggt                                                                                                                                                                                                  <  2:545726/247‑2 (MQ=255)
                                                  ggcattACAGTTTATAGCCCATTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGATCTGGAGAAAGCTTATAACAGATGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGCCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTTCATCGTTAAAGGctc                                                                                                                                                                                                <  2:560850/246‑1 (MQ=255)
                                                        aCAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGc                                                                                                                                                                                          <  2:621819/251‑1 (MQ=255)
                                                          aGTTGATAGCCCCTTCCCAGGTGGCGGAAAGCATATTCCGGCAAGCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAg                                                                                                                                                                                        <  2:601724/251‑1 (MQ=255)
                                                           gTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGt                                                                                                                                                                                      >  1:542746/1‑251 (MQ=255)
                                                            ttGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTg                                                                                                                                                                                     >  2:507077/1‑251 (MQ=255)
                                                                 aGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGTTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGa                                                                                                                                                                                 <  1:623773/250‑1 (MQ=255)
                                                                 aGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGACAAACCGCAAACACACCCCACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGa                                                                                                                                                                                 <  1:308039/250‑1 (MQ=255)
                                                                    cccgtCCCAGGGAGCGGGAAGCTTATTTCGGCAATCCAGGGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTAAAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCGCGAATGGTTCTTGTCTCATTGCCGCATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAATGCTCTGTGGCAGTACTGATCaa                                                                                                                                                                             <  2:225033/247‑1 (MQ=255)
                                                                      cTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCaa                                                                                                                                                                             >  1:121452/1‑249 (MQ=255)
                                                                                  ggggAAGCATATTTCGGAAATCCAGAGACAGCGGCTTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTa                                                                                                                                                               <  2:229080/250‑1 (MQ=255)
                                                                                       aGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGaa                                                                                                                                                         <  1:240032/251‑1 (MQ=255)
                                                                                       aGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGa                                                                                                                                                          >  1:131427/1‑250 (MQ=255)
                                                                                             tttcgaccatcaagtgaCAGCGGCGTTATCTGGCTCTGGAGCAATCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAACCAGACCCGACTCTCGAATGTTTCTTGTCTCATTGCCACATTCATAAGAACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAACGCGGAAACGATGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCACATACGAAGAGGGTAAAGAAATGATc                                                                                                                                                   >  2:368650/16‑251 (MQ=255)
                                                                                                cggaaaTCCAGAGCCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGAAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCctc                                                                                                                                                <  2:511209/247‑1 (MQ=255)
                                                                                                   cAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTcc                                                                                                                                             <  1:65325/251‑1 (MQ=255)
                                                                                                    aaacaagagACAGCGGAGTTATCTGGCTCTGGAGAAAGTTTATAACAGGGGGTAACCGCGCAGGGTGCTTGGCAAACCGAAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCGGATCAAAGACGAAGAGGGTAATGAAATGATCCTTTCCt                                                                                                                                            <  2:455346/246‑1 (MQ=255)
                                                                                                           agaCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCAATCCAAGAGAAAGCTTATTCACCAGGGTGAAGAAGCGGAAACGCTGTACTACATCGTTAAAGGCTATGGGGAAGTGCTGATAAAAGACGAGGAGGGTCAAGAAATGAACCTCTCCAATCTGa                                                                                                                                      >  2:84352/1‑250 (MQ=255)
                                                                                                                 cggcgTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAggg                                                                                                                               <  1:525689/251‑1 (MQ=255)
                                                                                                                      tgatttggttctGGAAAAAGCTTAAAACAGAGAATAACCGCGCATGGTCCTTGGCAACCCGCAAACATACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGACCCCCTCCTATCTGAATCAGGGTGAt                                                                                                                           <  2:575732/241‑1 (MQ=255)
                                                                                                                       tATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAtt                                                                                                                          <  2:258174/250‑1 (MQ=255)
                                                                                                                                ctGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGc                                                                                                                  <  1:84047/249‑1 (MQ=255)
                                                                                                                                  ggAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGCCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCg                                                                                                                 <  1:516928/248‑1 (MQ=255)
                                                                                                                                    aaaaaGCTTAGAACAGAGGATAACAGCGCATGGTGCTTGGCAAACCGCAATCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAAGGATCCGCTCCTATCTGAATCAGGCTGATTTTATTGGCGAACTg                                                                                                            <  2:170518/249‑1 (MQ=255)
                                                                                                                                     gaAAGCTTATAACAGAGGATAACCGCGCATGTTGCTTGGCAAACCGGAAACCGACCCGTCTCTCGAATGGTTCTTGTCTCAGTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGAGCAAAGACGAAGAGTGTAAAGAAATGATCCGCGCCTATCCGATGCAGGGTGATTTTATTGTCGAACTgg                                                                                                           <  2:472036/251‑1 (MQ=255)
                                                                                                                                               aaCAGAGGATAACCGCGCATGGTGCTTGGCAACCCTCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGTCTGTTTg                                                                                                  <  2:573696/250‑1 (MQ=255)
                                                                                                                                               aaCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTg                                                                                                  >  1:573696/1‑250 (MQ=255)
                                                                                                                                               aaCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTg                                                                                                  >  1:418132/1‑250 (MQ=255)
                                                                                                                                                         aaCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                         >  1:320085/1‑249 (MQ=255)
                                                                                                                                                         aaCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                         >  1:437556/1‑249 (MQ=255)
                                                                                                                                                          aCCGCGCATGGTGCTGGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCGCCTAGCTGAATCAGGGTGATTTTATTGGCGAACTGGGTCTGTTTGAAGAGGGCCAgg                                                                                      <  2:528464/251‑1 (MQ=255)
                                                                                                                                                            cgcgCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGCATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGa                                                                                     >  2:36168/1‑250 (MQ=255)
                                                                                                                                                              ggcATGGTTCTGGGCAAACCTCAAAATAACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACTCATCCAAGAGCAATCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGTCTGTTTGAAGAGGGCCAGGAACg                                                                                  <  2:168936/250‑1 (MQ=255)
                                                                                                                                                               gcATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGt                                                                                 >  1:121515/1‑251 (MQ=255)
                                                                                                                                                                             aaaCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACg                                                                   <  1:101200/251‑1 (MQ=255)
                                                                                                                                                                                        aGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCg                                                        >  1:364398/1‑251 (MQ=255)
                                                                                                                                                                                               aCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCtgtg                                                  >  1:206167/1‑250 (MQ=255)
                                                                                                                                                                                                   ccGAATGGTTCTTGTCCAATTCCCACATCCAAAAGTACCCATCCAAGAGCAGGCTTATTCACCGGGGTGAAAAAGCGGAACCGCGGTACTACATCGTTAAAGGCTTTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGGGGGCCAGGAACGAAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTg                                             <  2:371885/250‑1 (MQ=255)
                                                                                                                                                                                                     gAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGc                                           >  1:624784/1‑251 (MQ=255)
                                                                                                                                                                                                         gtttCTTGTCTCATTGCCACATTCAAAAGTCCCCACCCAAGAGAAAGCTTATCCCCCAGGGTGAAAAAGCGGAAACGCTGTACTACATCTTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGTTAAAGAAATGATCCTCTCCTATCTGAATGGGGGTGATTTTATTGGCGAACGGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGa                                        <  1:370932/248‑1 (MQ=255)
                                                                                                                                                                                                            tCTTGTCCCATTGCCACATTCATAAGTACCCACCCAAGAGCAAGCTTATTCCCCAGGGTAAAAAAGCGGAAACGCTGTACTACATCGTAAAAGGCTCGGTGGCAGTGCTGATCAAAGACGAAGAGGGAAAAAAAATGATCCTCTCCTTTCTGAATCAGGGTGTTTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAAt                                     <  1:482847/250‑1 (MQ=255)
                                                                                                                                                                                                             cTTGTCTCATTGCCACATTCATAAGTACCCATACAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGCAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCCTGGGTACGTGCGAAAACCGCCTGTGAAGTGGCAGAAAttt                                   >  2:67346/1‑251 (MQ=255)
                                                                                                                                                                                                               tGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGACGTGGCTGAAATTTCg                                 >  1:481580/1‑251 (MQ=255)
                                                                                                                                                                                                                 tctcATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGt                                >  1:359662/1‑250 (MQ=255)
                                                                                                                                                                                                                    cagtGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATCCACGAGGGGGAAAAAGCGAAAACGGTGTACTACATTTTTAGAGGTTCTGTGCAGGTGGTAAGAAAAGACGAAGGGGTAAAAGAAATGATCCTCGCTTATCTGAACAAGGGTTATTTTATTGGGGAGCGGGGCTTGTTTGAAGGGGGCGAGGAAGGTAGCGCGTGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTAc                              <  2:439123/246‑1 (MQ=255)
                                                                                                                                                                                                                      ttGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaa                           >  1:489146/1‑250 (MQ=255)
                                                                                                                                                                                                                        gCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaaaa                         >  1:68916/1‑250 (MQ=255)
                                                                                                                                                                                                                                 aTAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCa                >  1:404961/1‑250 (MQ=255)
                                                                                                                                                                                                                                           aTCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAgg     >  1:492818/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cAAGAGCAAGCTTATTCCCCAGGGTGAAAAAGCGGAAACGCTGTACTACATCTTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTAAATCAGGGTGTTTTTATTGGCGACCTGGGCCTGTTTGAGGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTaa  <  1:590995/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAA  >  NC_000913/3485959‑3486447

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