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

ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAA  >  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)
                                                           gTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑t                                                                                                                                                                                      >  1:542746/1‑251 (MQ=255)
                                                            ttGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑Tg                                                                                                                                                                                     >  2:507077/1‑251 (MQ=255)
                                                                 aGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGTTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGa                                                                                                                                                                                 <  1:623773/250‑1 (MQ=255)
                                                                 aGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGACAAACCGCAAACACACCCCACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGa                                                                                                                                                                                 <  1:308039/250‑1 (MQ=255)
                                                                    cccgtCCCAGGGAGCGGGAAGCTTATTTCGGCAATCCAGGGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTAAAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCGCGAATGGTTCTTGTCTCATTGCCGCATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAATGCTCTGTGGCAG‑TACTGATCaa                                                                                                                                                                             <  2:225033/247‑1 (MQ=255)
                                                                      cTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCaa                                                                                                                                                                             >  1:121452/1‑249 (MQ=255)
                                                                                  ggggAAGCATATTTCGGAAATCCAGAGACAGCGGCTTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTa                                                                                                                                                               <  2:229080/250‑1 (MQ=255)
                                                                                       aGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑Gaa                                                                                                                                                         <  1:240032/251‑1 (MQ=255)
                                                                                       aGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑Ga                                                                                                                                                          >  1:131427/1‑250 (MQ=255)
                                                                                             tttcgaccatcaagtgaCAGCGGCGTTATCTGGCTCTGGAGCAATCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAACCAGACCCGACTCTCGAATGTTTCTTGTCTCATTGCCACATTCATAAGAACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAACGCGGAAACGATGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCACATACGAAGAGGGTAAA‑GAAATGATc                                                                                                                                                   >  2:368650/16‑251 (MQ=255)
                                                                                                cggaaaTCCAGAGCCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGAAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCctc                                                                                                                                                <  2:511209/247‑1 (MQ=255)
                                                                                                   cAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTcc                                                                                                                                             <  1:65325/251‑1 (MQ=255)
                                                                                                    aaacaagagACAGCGGAGTTATCTGGCTCTGGAGAAAGTTTATAACAGGGGGTAACCGCGCAGGGTGCTTGGCAAACCGAAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCGGATCAAAGACGAAGAGGGTAAT‑GAAATGATCCTTTCCt                                                                                                                                            <  2:455346/246‑1 (MQ=255)
                                                                                                           agaCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCAATCCAAGAGAAAGCTTATTCACCAGGGTGAAGAAGCGGAAACGCTGTACTACATCGTTAAAGGCTATGGGGAAG‑TGCTGATAAAAGACGAGGAGGGTCAA‑GAAATGAACCTCTCCAATCTGa                                                                                                                                      >  2:84352/1‑250 (MQ=255)
                                                                                                                 cggcgTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAggg                                                                                                                               <  1:525689/251‑1 (MQ=255)
                                                                                                                      tgatttggttctGGAAAAAGCTTAAAACAGAGAATAACCGCGCATGGTCCTTGGCAACCCGCAAACATACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGACCCCCTCCTATCTGAATCAGGGTGAt                                                                                                                           <  2:575732/241‑1 (MQ=255)
                                                                                                                       tATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGAtt                                                                                                                          <  2:258174/250‑1 (MQ=255)
                                                                                                                                ctGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGc                                                                                                                  <  1:84047/249‑1 (MQ=255)
                                                                                                                                  ggAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGCCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCg                                                                                                                 <  1:516928/248‑1 (MQ=255)
                                                                                                                                    aaaaaGCTTAGAACAGAGGATAACAGCGCATGGTGCTTGGCAAACCGCAATCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAAGGATCCGCTCCTATCTGAATCAGGCTGATTTTATTGGCGAACTg                                                                                                            <  2:170518/249‑1 (MQ=255)
                                                                                                                                     gaAAGCTTATAACAGAGGATAACCGCGCATGTTGCTTGGCAAACCGGAAACCGACCCGTCTCTCGAATGGTTCTTGTCTCAGTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGAGCAAAGACGAAGAGTGTAAA‑GAAATGATCCGCGCCTATCCGATGCAGGGTGATTTTATTGTCGAACTgg                                                                                                           <  2:472036/251‑1 (MQ=255)
                                                                                                                                               aaCAGAGGATAACCGCGCATGGTGCTTGGCAACCCTCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGTCTGTTTg                                                                                                  <  2:573696/250‑1 (MQ=255)
                                                                                                                                               aaCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTg                                                                                                  >  1:573696/1‑250 (MQ=255)
                                                                                                                                               aaCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTg                                                                                                  >  1:418132/1‑250 (MQ=255)
                                                                                                                                                         aaCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                         >  1:320085/1‑249 (MQ=255)
                                                                                                                                                         aaCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                         >  1:437556/1‑249 (MQ=255)
                                                                                                                                                          aCCGCGCATGGTGCTGGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCGCCTAGCTGAATCAGGGTGATTTTATTGGCGAACTGGGTCTGTTTGAAGAGGGCCAgg                                                                                      <  2:528464/251‑1 (MQ=255)
                                                                                                                                                            cgcgCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGCATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGa                                                                                     >  2:36168/1‑250 (MQ=255)
                                                                                                                                                              ggcATGGTTCTGGGCAAACCTCAAAATAACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACTCATCCAAGAGCAATCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGTCTGTTTGAAGAGGGCCAGGAACg                                                                                  <  2:168936/250‑1 (MQ=255)
                                                                                                                                                               gcATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGt                                                                                 >  1:121515/1‑251 (MQ=255)
                                                                                                                                                                             aaaCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACg                                                                   <  1:101200/251‑1 (MQ=255)
                                                                                                                                                                                        aGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCg                                                        >  1:364398/1‑251 (MQ=255)
                                                                                                                                                                                               aCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCtgtg                                                  >  1:206167/1‑250 (MQ=255)
                                                                                                                                                                                                   ccGAATGGTTCTTGTCCAATTCCCACATCCAAAAGTACCCATCCAAGAGCAGGCTTATTCACCGGGGTGAAAAAGCGGAACCGCGGTACTACATCGTTAAAGGCTTTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGGGGGCCAGGAACGAAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTg                                             <  2:371885/250‑1 (MQ=255)
                                                                                                                                                                                                     gAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGc                                           >  1:624784/1‑251 (MQ=255)
                                                                                                                                                                                                         gtttCTTGTCTCATTGCCACATTCAAAAGTCCCCACCCAAGAGAAAGCTTATCCCCCAGGGTGAAAAAGCGGAAACGCTGTACTACATCTTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGTTAAA‑GAAATGATCCTCTCCTATCTGAATGGGGGTGATTTTATTGGCGAACGGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGa                                        <  1:370932/248‑1 (MQ=255)
                                                                                                                                                                                                            tCTTGTCCCATTGCCACATTCATAAGTACCCACCCAAGAGCAAGCTTATTCCCCAGGGTAAAAAAGCGGAAACGCTGTACTACATCGTAAAAGGCTCGGTGGCAG‑TGCTGATCAAAGACGAAGAGGGAAAA‑AAAATGATCCTCTCCTTTCTGAATCAGGGTGTTTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAAt                                     <  1:482847/250‑1 (MQ=255)
                                                                                                                                                                                                             cTTGTCTCATTGCCACATTCATAAGTACCCATACAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGCAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCCTGGGTACGTGCGAAAACCGCCTGTGAAGTGGCAGAAAttt                                   >  2:67346/1‑251 (MQ=255)
                                                                                                                                                                                                               tGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGACGTGGCTGAAATTTCg                                 >  1:481580/1‑251 (MQ=255)
                                                                                                                                                                                                                 tctcATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGt                                >  1:359662/1‑250 (MQ=255)
                                                                                                                                                                                                                    cagtGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATCCACGAGGGGGAAAAAGCGAAAACGGTGTACTACATTTTTAGAGGTTCTGTG‑CAGGTGGTAAGAAAAGACGAAG‑GGGTAAAAGAAATGATCCTCGCTTATCTGAACAAGGGTTATTTTATTGGGGAGCGGGGCTTGTTTGAAGGGGGCGAGGAAGGTAGCGCGTGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTAc                              <  2:439123/246‑1 (MQ=255)
                                                                                                                                                                                                                      ttGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaa                           >  1:489146/1‑250 (MQ=255)
                                                                                                                                                                                                                        gCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaaaa                         >  1:68916/1‑250 (MQ=255)
                                                                                                                                                                                                                                 aTAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCa                >  1:404961/1‑250 (MQ=255)
                                                                                                                                                                                                                                           aTCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAgg     >  1:492818/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cAAGAGCAAGCTTATTCCCCAGGGTGAAAAAGCGGAAACGCTGTACTACATCTTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTAAATCAGGGTGTTTTTATTGGCGACCTGGGCCTGTTTGAGGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTaa  <  1:590995/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAG‑TGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAA  >  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: ‑