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

ACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCC  >  NC_000913/3485967‑3486451
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
aCATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAcc                                                                                                                                                                                                                                             >  1:186952/1‑250 (MQ=255)
  attgtgttacgtattggaTGCAAAGGCAGTCCCATGACTGTGCAGTACAGGTGATAGCCCCTTCCCAGGTAGCGGGAAGCAAATTTCGGCACTCTAGAGACAGCGGCGTTATCTGGCTCTGGAAAAAGCTTATAACAGAGGATAACCGCGCTGGGTGCTTGGCGACCCGCAAACAGACCCGACTCTCGAATTGTGCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAgg                                                                                                                                                                                                                                          <  2:487674/234‑1 (MQ=255)
  atgtaTGTTCTGCATGTAGGCAAAGGACGTCACATTACCGTGCAGTACATTGTCTACCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAACCCAGAGACGGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGTATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGCCCCGACTCTCGAATTGTTCTTGTCTCATTGCCACATTCATAAGTCCCCATCCAAGAGCAAGCTTATTCACCAgg                                                                                                                                                                                                                                          <  1:77853/247‑1 (MQ=255)
             gCATGTATGCAAAGGACGGCACATTGCCGTGCAGTCCAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCAGGGTGCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAg                                                                                                                                                                                                                                 <  1:343499/249‑1 (MQ=255)
               aTGTATGCAAAGGACGTCACATTCCCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCgg                                                                                                                                                                                                                              >  1:135848/1‑250 (MQ=255)
                 gtacgCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTGGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCGGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGaa                                                                                                                                                                                                                            <  2:373675/246‑1 (MQ=255)
                          ggACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGtact                                                                                                                                                                                                                  <  2:198174/251‑1 (MQ=255)
                                 acaTTACCGTGCCGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTTGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTTGCTAACCGCAAACCAACCCTCCTCCCGAATGGTTCTTGTCTCATTCCCACATTCATAAGTACCCATCCCAGAGCAAGCTTATTCACCAGGTATAAACAACGTAAACGCGGTACTACATCGt                                                                                                                                                                                                           >  2:38500/1‑251 (MQ=255)
                                   aTTACCGTGCCGGACAGTTGATAGCCCCTTCCCAGGAAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTAGCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCCCAAACAGACCCGCCTCTCGAATGGTTCTTGTCTCCTTGCCACATTCATAAGAACCCATCCCAGAGCATGCTTATTCACACGGGTGAAAACTCGGAAACGCTGTACTCCATCGtt                                                                                                                                                                                                          >  2:369141/1‑250 (MQ=255)
                                      aCCGTGCAGTACAGTTGATAGCCCCTTCCCCGGTAGCGGGACGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGCGANAAAGAGGAAACGCTGTACTAGATCGTTaa                                                                                                                                                                                                        >  2:224937/1‑249 (MQ=255)
                                            cAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCtg                                                                                                                                                                                                <  1:108076/251‑1 (MQ=255)
                                                        aTAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCCGGGTGAAAAAGCGGAACCGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTg                                                                                                                                                                                     >  2:338844/1‑251 (MQ=255)
                                                           gCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATc                                                                                                                                                                                  <  2:135848/251‑1 (MQ=255)
                                                                 tCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACCTCGTTAAAGGCTCTGTGGCTGTGCTGATCAAAGAc                                                                                                                                                                            >  1:252397/1‑251 (MQ=255)
                                                                  cccAGGTAACGGGAAGTATATTTCGGGAAACTAAAGACAGCGGCGTTATCAGGCGGGGGAGAAAGCTTATAACGGAGGAGAACCCCGCATGGTGCTTGGCAAACCCCAAACAAACCCGACTCTCGAATGATTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAGGCGGAAACGCTGTACGACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACg                                                                                                                                                                           <  2:287965/251‑1 (MQ=255)
                                                                      ggTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAc                                                                                                                                                                                                                          >  1:2675/1‑200 (MQ=255)
                                                                      ggTAGCGGGAAGAATATTTCGGAAATCCAAAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTAAAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAc                                                                                                                                                                                                                          <  2:2675/200‑1 (MQ=255)
                                                                              gAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTaa                                                                                                                                                                 >  1:101533/1‑249 (MQ=255)
                                                                                         cGGCAATCCAGAGACAGCGGCGTTATCTGGTNCTGGAGAAAGCTTATAACAGAGAGAAACCGCGCATGGTGCTTGGCAAACCGCAGACAGCCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTGCCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAATGCTCTGTGGCAGTTCTGATCAAAGACGAAGAGGGTAAAGAAATGATCctc                                                                                                                                                    <  2:236960/251‑1 (MQ=255)
                                                                                         cGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCctc                                                                                                                                                    >  2:439388/1‑251 (MQ=255)
                                                                                           gCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCtctc                                                                                                                                                  >  1:30124/1‑251 (MQ=255)
                                                                                             aaTCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGGGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTccg                                                                                                                                                <  2:295100/251‑2 (MQ=255)
                                                                                             aaTCCAGAGACAGCGGCGTTATCGGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGAAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCt                                                                                                                                                <  2:122059/251‑1 (MQ=255)
                                                                                                    agaCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGaa                                                                                                                                         <  1:439388/251‑1 (MQ=255)
                                                                                                              gTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAt                                                                                                                               >  1:169047/1‑251 (MQ=255)
                                                                                                                     ggCTGTGGAGAAAGCTTATAACAGGAGATAACGGCGCATGGTGATTGGCAAACCGCAAACCGGCCCGACTCTCGAATGGTTCTTGTCTCATTGCCGCATTCATAAGTACCCACCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGAGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGATGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTAACTGAATCAGGGTGATTTTATTg                                                                                                                        <  2:355093/251‑1 (MQ=255)
                                                                                                                      gCTCTGGAGAAAGCTTATAACAGAGGCTAACCGAGCATGGTGCGTGGCAAACCGCAAACAGACCCGTCTCTCGAATGGTTCTTGTCTCATGGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTGTATTgg                                                                                                                       <  2:36563/251‑1 (MQ=255)
                                                                                                                                       ataaaagtggatcaccgccgttgttGCTTTGCCAACCGAAAACAGTCCCGACGCTCGAATGGTTCTTGTCTCACTGCCAAAGTCATAAGTACCCATCCAAGAGCAAGCTTATTCTCCAGGGTGAAAAAGCGGAATCGCGGTGCTACATCGATAAAGGCTCTGTGGGAGTGCTGATCAAAGACGAAGAGGATAAAGAAACGATCATCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGtt                                                                                                        <  1:474979/226‑1 (MQ=255)
                                                                                                                                       acccgactctttgccgctctgtGGTGTTTGGCAAACCACAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTCCCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCTTGt                                                                                                         <  2:104759/228‑1 (MQ=255)
                                                                                                                                        aaCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTg                                                                                                      >  2:344437/1‑250 (MQ=255)
                                                                                                                                          cAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAg                                                                                                                                                                                          <  1:253220/164‑1 (MQ=255)
                                                                                                                                          cAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAg                                                                                                                                                                                          >  2:253220/1‑164 (MQ=255)
                                                                                                                                             agGATACCCGCTCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAgg                                                                                                <  2:186952/251‑1 (MQ=255)
                                                                                                                                                 aTAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                             >  1:267857/1‑251 (MQ=255)
                                                                                                                                                  tAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATACAAGAGCAAGCTTATTCACCAGGGTGAAACAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGCGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                             >  2:488376/1‑250 (MQ=255)
                                                                                                                                                    aCCGCGCATGGTGCTTAGCAAACCGCAAACCGACCCGACTCTCGACTGGTTCTTGTCTCATTGCCACATTCATAAGTACCCCTCCAAGAGCAAGCTTATTCACCAGGGGGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCGCTGTGGCAGCGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGCGATTTTATTGGCGAACTGGGCCTGTTTGCAGAGGGcccgg                                                                                          >  2:455004/1‑248 (MQ=255)
                                                                                                                                                           aTGGAGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGACGAACTGGGCCTGt                                                                                                         >  1:104759/1‑229 (MQ=255)
                                                                                                                                                                gCTTGGCAAACCGCAAACTGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGCTCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCAt                                                                              >  2:420897/1‑251 (MQ=255)
                                                                                                                                                                           cGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAagag                                                                                                                                                                      <  1:449366/152‑1 (MQ=255)
                                                                                                                                                                           cGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAagag                                                                                                                                                                      >  2:449366/1‑152 (MQ=255)
                                                                                                                                                                             cAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTgg                                                                                                               <  1:370558/205‑1 (MQ=255)
                                                                                                                                                                             cAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTgg                                                                                                               >  2:370558/1‑205 (MQ=255)
                                                                                                                                                                              aaaCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCaaa                                                                                                                                                                               <  2:207778/140‑1 (MQ=255)
                                                                                                                                                                              aaaCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCaaa                                                                                                                                                                               >  1:207778/1‑140 (MQ=255)
                                                                                                                                                                                aCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAAc                                                              <  1:245159/251‑1 (MQ=255)
                                                                                                                                                                                          ctctCGAATGGTTCTTGTCTCATTGCCACATTCATACGTACCCATCCATGAGCAAGCTTATTCACCAGGGCGAAAAAGCGGAAAAGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGAACAACGACGAAGAGGGTAAAGAAATGCTCCTCTCCTATCTGAATCAGGGCGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTATCGCCTGGGTACGTGCGAAAACAGCCtgtat                                                     >  2:181270/1‑248 (MQ=255)
                                                                                                                                                                                                       tttGTCTCATTGCCACATTCATAAGTACCCACCAAAGAGCAAGCTTATTCCCCAGGGTAAAAAAGCGGAAACGTTGTACTACATCGTTAAGGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGAAAAAGAAATGACCCCCCCCTATCTGAATCAGGGTGTTTTTATTGGCGAACTGGGCTTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTCCGGGCAAAAACCGCCTGTGAAGTGGCTGAAAtt                                        <  1:496104/249‑1 (MQ=255)
                                                                                                                                                                                                        ttGTCTCATTGCCAAATTCATAAGTAACCATCCACGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGATGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGCTCAAAGACGAAGAGGGTAACGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGCACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTc                                      >  2:43417/1‑251 (MQ=255)
                                                                                                                                                                                                           tctcATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGt                                    >  1:433106/1‑250 (MQ=255)
                                                                                                                                                                                                            ctcATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTAc                                  >  2:136222/1‑251 (MQ=255)
                                                                                                                                                                                                               aTTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaa                                >  1:362300/1‑250 (MQ=255)
                                                                                                                                                                                                                ttGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTg                                                                                                                                 <  1:417039/152‑1 (MQ=255)
                                                                                                                                                                                                                ttGCCACATTCATAAGTAACCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTg                                                                                                                                 >  2:417039/1‑152 (MQ=255)
                                                                                                                                                                                                                      caTTCATAAGTACCCACCAAAGAGAAAGCTTATTCACCAGGGTAAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATAAAAGAGGAAGAGGGTAAAGAAATGATCCTCCCTTATCTGAATCAGGGTTATTTTATGGGCGAACTGGGCTGGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTc                        <  1:420897/251‑1 (MQ=255)
                                                                                                                                                                                                                                 acacATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCTCCAATTGAtt             >  2:45864/4‑251 (MQ=255)
                                                                                                                                                                                                                                    cacccAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCGCTGTGGCAGTGCTTATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATc                                                                                                                                       <  2:52258/123‑1 (MQ=255)
                                                                                                                                                                                                                                    cATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGCAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATc                                                                                                                                       >  1:52258/1‑126 (MQ=255)
                                                                                                                                                                                                                                    cATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTc            >  1:330580/1‑249 (MQ=255)
                                                                                                                                                                                                                                           gagCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGATCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGGGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAcc   >  1:319939/1‑251 (MQ=255)
                                                                                                                                                                                                                                            agCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGAAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAccc  >  2:309538/1‑251 (MQ=255)
                                                                                                                                                                                                                                            agCAAGCTTATCCACCAGGGTGAAAAAGGGAAAACGCTGTACTACATCGTAAAAGGCTCTGTGGCAGGGTTGATCAAAGACGAAGAGGTTAAAGAAATGATCCTCTCCTATCTGAAGCGGGGTGATTTTATTGGCGAACTGGGCTTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAccc  <  2:30124/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
ACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCC  >  NC_000913/3485967‑3486451

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: