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.2% 211.2 / NA 56T29K (ACG→AAG) crpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (31/24);  minor base G (0/1);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGT‑CACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAC  >  NC_000913/3485959‑3486449
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
aCAGTAATACATTGATGTACTGCATGTATGCAAAGGACGT‑CACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCt                                                                                                                                                                                                                                                       <  1:31544/250‑1 (MQ=255)
   ggacaaCATTGATGTCCTGCCGGTCTGCAAAGG‑CGTACACATTCCCGTGCAGTACAGTTTATAGCCCCTTCCAAGGTCGCGGGAAGAATATTCCGGCAATCCAGGGACATGGGCGTTACCTGGCTCTGGAGAAGGCTTATAACAGAGGATAACCGGGAGTGGTGCTTGGCAAACCGCGAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCGCATTCATAAGTACCCAGCCAAGAGCAAGattatt                                                                                                                                                                                                                                                   <  2:30844/246‑1 (MQ=255)
   gTAATACATTGATGTACTGCATGTATGCAAAGGACGT‑CACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTAtt                                                                                                                                                                                                                                                   >  1:332919/1‑251 (MQ=255)
   gTAATACATTGATGTACTGCATGTATGCAAAGGACGT‑CACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTAtt                                                                                                                                                                                                                                                   <  1:289046/251‑1 (MQ=255)
       tACATTGATGTACTGCATGTATGCAAAGGACGT‑CACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAcc                                                                                                                                                                                                                                               >  1:441106/1‑251 (MQ=255)
         cATTGATGTACTGCATGTATGCAAAGGACGT‑CACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAg                                                                                                                                                                                                                                             <  1:15416/251‑1 (MQ=255)
          aTTGATGTACTGCATGTATGCAAAGGACGT‑CACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAgg                                                                                                                                                                                                                                            <  2:508468/251‑1 (MQ=255)
           ttGATGTACTGCATGTATGCAAAGGACGT‑CACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAggg                                                                                                                                                                                                                                           >  2:518370/1‑251 (MQ=255)
                     gCATGTATGCAAAGGACGT‑CACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGc                                                                                                                                                                                                                                 <  1:103498/250‑1 (MQ=255)
                       aTGTATGCAAAGGACGT‑CACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCgg                                                                                                                                                                                                                               >  1:135643/1‑250 (MQ=255)
                        tGTATGCAAAGGACGT‑CACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGa                                                                                                                                                                                                                              >  1:203351/1‑250 (MQ=255)
                                      gT‑CACATTACCGTGCATTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACAt                                                                                                                                                                                                               <  2:243767/251‑1 (MQ=255)
                                         cacaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAACCCAGAGACAGCGGCTTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCGAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCg                                                                                                                                                                                                             <  1:414154/251‑1 (MQ=255)
                                          aaaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCGGGTAACGGGAAGCATATTTCGGCAATCCAGAGACAGAGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCAACATTCATAAGTACACATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAACAGCGGTACTACAt                                                                                                                                                                                                               >  2:216555/3‑248 (MQ=255)
                                                ccGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAgg                                                                                                                                                                                                     >  2:206904/1‑251 (MQ=255)
                                                  gTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGct                                                                                                                                                                                                   >  1:292888/1‑251 (MQ=255)
                                                     cAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCtg                                                                                                                                                                                                <  2:400549/251‑1 (MQ=255)
                                                           aGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCa                                                                                                                                                                                           >  1:358573/1‑250 (MQ=255)
                                                               gATAGCCCCTACCCAGGTAGCGGGAAGCATATTTCGGCAATCAAGAGCCAGCGGCGTTATCGGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAAAAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCt                                                                                                                                                                                      <  2:13943/251‑1 (MQ=255)
                                                                   gCCCCTTCCCAGTTAGCGGGAAGCATATGTCGGCAATCCAGCGACAGCGGCGTTATCTGGCTCTGGAAAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAATCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATc                                                                                                                                                                                  <  2:99558/251‑1 (MQ=255)
                                                                                 aGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAggg                                                                                                                                                                    <  1:81852/251‑1 (MQ=255)
                                                                                   cGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACTGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTa                                                                                                                                                                  <  1:251283/251‑1 (MQ=255)
                                                                                     ggAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTaaa                                                                                                                                                                <  1:141833/251‑1 (MQ=255)
                                                                                         gCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTGCTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGaa                                                                                                                                                             >  1:356314/1‑250 (MQ=255)
                                                                                                   gCAATCCAGAGCCAGCGGCGTTATCTGGCTCTGGAGAATGCTTATAACAGGGGAAAACCGCGCATGGTGCGTGCCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATGCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAAGGATGCtctc                                                                                                                                                  <  1:514631/251‑1 (MQ=255)
                                                                                                            ggaCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGAAAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGaa                                                                                                                                         <  1:65308/250‑1 (MQ=255)
                                                                                                                aGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATc                                                                                                                                       >  2:33977/1‑249 (MQ=255)
                                                                                                                        tATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAtt                                                                                                                              <  1:221542/250‑1 (MQ=255)
                                                                                                                                   ggAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGaa                                                                                                                                                                         <  1:304374/196‑1 (MQ=255)
                                                                                                                                   ggAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGaa                                                                                                                                                                         >  2:304374/1‑196 (MQ=255)
                                                                                                                                    gagaAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑Ct                                                                                                                >  1:472504/1‑251 (MQ=255)
                                                                                                                                            ttATAACAGAGGATAACCGCGCATGGAGCTTGACAAACCGCAACCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCAACATTCATAAGAACCCATCAAAGAGCAAGCTTATTAACCAGGGT‑GAAAAAGCGGAAAAGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGCGCTGATCAAAGACGAAAAGGGCAAAGAAATGATCCTCTCCTATCTGAATCGTGGTGATTTTATTGGCGAA‑CTGGGCCTGt                                                                                                        >  2:398029/1‑251 (MQ=255)
                                                                                                                                               taACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTg                                                                                                     >  1:425407/1‑251 (MQ=255)
                                                                                                                                                        aTAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGcc                                                                                            <  1:33977/251‑1 (MQ=255)
                                                                                                                                                                gcATGGAGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGt                                                                                   >  2:115349/1‑251 (MQ=255)
                                                                                                                                                                  aTGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAg                                                                                 <  1:245053/251‑1 (MQ=255)
                                                                                                                                                                         ttGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATgg                                                                          >  1:137635/1‑251 (MQ=255)
                                                                                                                                                                            gCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGt                                                                        >  1:493588/1‑250 (MQ=255)
                                                                                                                                                                                  cGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGc                                                                                             >  2:5985/1‑224 (MQ=255)
                                                                                                                                                                                  cGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGc                                                                                             <  1:5985/224‑1 (MQ=255)
                                                                                                                                                                                   gCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCg                                                                 >  1:65992/1‑250 (MQ=255)
                                                                                                                                                                                    cAAACAGACCCGACTCTCCAATGGTTCTTGTCTCATTGCCAAATTCATAAGTACCCATCCAAGCGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGATGATAAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCCAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCGac                                                               >  2:11666/1‑250 (MQ=255)
                                                                                                                                                                                          gACCCGACTCTCGAATGTTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCGAAAACCGc                                                         >  2:114695/1‑251 (MQ=255)
                                                                                                                                                                                                                      aTTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaa                              >  1:392994/1‑250 (MQ=255)
                                                                                                                                                                                                                        tgcaacattaataaggccccacacAAGAGCAGGCTTATTCCCAAGTGT‑GAAAAAGCGGAAACGCTGTAATTAATCTAT‑AAAGGCTGTGTGGAAGTGCGGACCAAAGACGAGGAGGGAAAAAAAATGATCCGCCCTTATCTGACTCAGGGTGTTTTTATTGGAGAA‑ATGGGCCTGTTGGAAGAGGGC‑AGGAAACGAAGGGCTGGGGTAAGAACAAAAGCCGCCTGTGAAGTGGCTGAAATTTCGTACaaaa                            <  2:281143/227‑1 (MQ=255)
                                                                                                                                                                                                                         cccACATTCATAATTACCCACCCAAGAGAAATCTTATTCACCAGGGT‑GAAAAAGCGAAACCGCTGTACTAAATCGTT‑AAAGGCTCGGGGGCAGTGCTGATAAAAGACGAAGAGGGAAAAGAAATGATCCTCTCCTAGCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCTTGTTTGAAGAGGCCCAGGAA‑CGTAGCGCATGGGTGCGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaaaaa                          <  1:212646/250‑1 (MQ=255)
                                                                                                                                                                                                                          ccgaatacaaaattACCCACCAAAGAGAAAGCTTATTCAC‑AGGGTTGAAAAAGCGGAACCGCTTTACTACA‑CGTTTAAAGGCGTTGTGGCAGTTCTGATCAAAGACGAAGAGGGAAAAGAAATTATCCTCTCCTATCTTAATCAGGGTGATTTTAT‑GGGGAATCTGGGCTTTTTTGAAGAGGGCCATGAA‑CGTGCCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAAt                         <  1:55172/238‑1 (MQ=255)
                                                                                                                                                                                                                             caTTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTc                      >  1:241544/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAAt                >  2:287557/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGCTTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAAtt               >  2:15381/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaGTACCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGAGGAAACGATAAACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTAAATCAGGGGGATTTTATTGGCGAA‑CTGGGCATGTTTGAAGAGGGCCAGGAA‑CGAAGCGCATGGGAACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCCCCAAtt               >  1:145696/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gaccccacccaAGCGCAAGCTTTTACACCAGGGT‑GAAAAGGGGGAAACGCTATACAACATCGTT‑AAAGGCTCTGGGGCAGTGCTGATCAAAGACGAAGAGGATAAAGAAATGACCCTCTCTTATCTGAAGCGGGTTGATTTTATGGGCGAA‑CTGGGCCTGTTGGAAGAGGGCCAGGAC‑CTTAGCGGAGGGATACGGGAGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGa             <  2:65992/248‑1 (MQ=255)
                                                                                                                                                                                                                                        aCCCATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGAtt           >  1:514150/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cATCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTc          >  1:93992/1‑249 (MQ=255)
                                                                                                                                                                                                                                            aTCCAAGAGCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAgg       >  1:219439/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 agagCAAGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAc  >  2:459864/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGT‑CACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGT‑GAAAAAGCGGAAACGCTGTACTACATCGTT‑AAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAGAGGGCCAGGAA‑CGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAC  >  NC_000913/3485959‑3486449

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