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

ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGT  >  NC_000913/3485959‑3486445
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aCAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACTGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCGGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGAAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCtt                                                                                                                                                                                                                                                  <  2:237998/251‑1 (MQ=255)
      aTACATTGATGTACTGCATGTATTCAAAGGACGTCCCATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGTTAGCGGGAAGCATA‑TTGCGGCATTCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATACCAGAGGATAACCGCGCAGGGTGCTTGGCAAACCGCAAACAGTCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTCCCCATCCAAGAGCAAGCTTATTCAc                                                                                                                                                                                                                                            <  1:566806/251‑1 (MQ=255)
        aCATTGATGTACTGTATGTATGCAAAGGTCGTCACATTCCCGTGCAGTACAGTTGATAGCCCCTTCCCAGGGAGCGGGAAGCATA‑TTGCGGCAATCCAGAGACAGCGGCTTTATCGGGCTCTGGAGAAAGCTTATAACAGGGGATAACCTCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAcc                                                                                                                                                                                                                                           <  2:413688/250‑1 (MQ=255)
          aTTGATGTACTGCATGTATGCAAAGGACGTCCCATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAgg                                                                                                                                                                                                                                        <  2:73152/251‑1 (MQ=255)
             gATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTg                                                                                                                                                                                                                                     <  1:387772/251‑1 (MQ=255)
                      cgtGTATGCAAAGGACGTCACGTTACCGTGCAGTACAGTTGATAGCCCCTTCACAGCTACCGGGAAGCATA‑TTTCGGCAATCCAAAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGCGGATAACCGCGCATGGAGCTGGGCAAACCGCCAACCGACCCCACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGACCCCATCCAAGAGCAAGCTTATTCACCAGGGTGGAAAAG‑Cgg                                                                                                                                                                                                                           >  2:51901/3‑251 (MQ=255)
                       aTGTATGAAAAGGACGTCTCATTTCCGTGCAGTACTGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAAACAGCGGCGTTATCGGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGa                                                                                                                                                                                                                          <  1:287426/251‑1 (MQ=255)
                                    aCGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGtactac                                                                                                                                                                                                            <  2:244713/251‑1 (MQ=255)
                                      gTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGAAAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCAATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCAGTACTACAt                                                                                                                                                                                                          >  1:298516/1‑251 (MQ=255)
                                      gTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCACAGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCAGTACTACAt                                                                                                                                                                                                          >  2:219587/1‑251 (MQ=255)
                                          cattaacatGCAGTACAGCTGATAGCCCCTTCCCAGGTAGCGGGACGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTCTCTGACTCTGGAGACAGCTTATAACAGAGGATAACCGCGCATAGCGCTTGGCAAACCGCAAACAGACACGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAATTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGTAA‑CGCTGTACTACATCGtt                                                                                                                                                                                                      >  2:161678/9‑251 (MQ=255)
                                            ttacggcggtgTACAGTTGATGTCCCCTCCCCAGTTGGAGGGAAGCATA‑TTCCGGCAATCCAGAGACAGCGGCGTTATCGGGCGTTGGATAAATCTTATAACAGGGGTAAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGTAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTa                                                                                                                                                                                                     <  2:125371/240‑1 (MQ=255)
                                                      gTACAGTTGATCGCCCCTTCCCAGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCtgt‑g                                                                                                                                                                                         >  1:455705/1‑251 (MQ=255)
                                                          aGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGc                                                                                                                                                                                       >  1:418559/1‑249 (MQ=255)
                                                                  gCCCCTTCCCAGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATc                                                                                                                                                                             >  1:241275/1‑251 (MQ=255)
                                                                        tCCCAGGTAGCGGGAAGCATA‑TTGCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGa                                                                                                                                                                        <  1:569364/250‑1 (MQ=255)
                                                                          ccAGGTAGCGGGAA‑CGTAATTACGGAAATCCAGAGACAGCTGCGTTATCTGGCTCTGGAGAAACCTTATAACAGAGGATAACGGCTCATTGTGCTTGCCACACCGCAACCAGACCCGACTCTCGAATTGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAT‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGa                                                                                                                                                                     <  2:46056/251‑1 (MQ=255)
                                                                            aGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAag                                                                                                                                                                   <  1:187495/251‑1 (MQ=255)
                                                                            aGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAag                                                                                                                                                                   <  1:463865/251‑1 (MQ=255)
                                                                                  cGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGt                                                                                                                                                              >  1:508774/1‑250 (MQ=255)
                                                                                     gAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATTTGGCTCTGGAGAAAGCTTATAACAAGGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATTGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAACACGAAGAGTGTAAAg                                                                                                                                                          <  2:437100/251‑1 (MQ=255)
                                                                                       aGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGaa                                                                                                                                                        <  1:412642/251‑1 (MQ=255)
                                                                                                ccGGCAATCCAGAGACGGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCct                                                                                                                                                <  1:431224/250‑1 (MQ=255)
                                                                                                 cGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCct                                                                                                                                                <  1:213530/250‑1 (MQ=255)
                                                                                                                 gcggcgTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTAGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAgg                                                                                                                               >  1:337114/1‑251 (MQ=255)
                                                                                                                      gTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGAATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑AGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGAACAAAGACGTAGAGGGTAAAGAAATCATCCGCTCCTATCTGAAGCAGGGGGAt                                                                                                                          >  2:114760/1‑251 (MQ=255)
                                                                                                                      gTTATCTGGCACTGGAGAACGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGACTGGATCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTC‑GTCGGCAGTGCGGATCAAAGACGACGCGGGTAAAGAAATGATCCTCTCCTATCGGAATCAGGGTGAt                                                                                                                          >  2:124997/1‑251 (MQ=255)
                                                                                                                             ggCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTg                                                                                                                   >  1:131573/1‑251 (MQ=255)
                                                                                                                             ggCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTg                                                                                                                   <  1:125739/251‑1 (MQ=255)
                                                                                                                                tctGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCg                                                                                                                >  1:152025/1‑251 (MQ=255)
                                                                                                                                                    gagGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAag                                                                                              >  1:267921/1‑249 (MQ=255)
                                                                                                                                                        ataacccagcATGGTGCTTTGCAAACCGCAAACAGACCCGACTCTCGAATTGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTTAAGAGGGcc                                                                                        <  2:204437/243‑1 (MQ=255)
                                                                                                                                                         tAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                        >  2:89063/1‑250 (MQ=255)
                                                                                                                                                             cgcgCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACCTCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGGGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCCGGac                                                                                   >  2:91484/1‑250 (MQ=255)
                                                                                                                                                             cgcgCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGaa                                                                                   <  1:198834/251‑1 (MQ=255)
                                                                                                                                                             agcgcATGGTGCTTGGTAAACCGCAAACAGACCCGACTCTCGCATGGATCTTGTCTCATTGCCAAATGCATAAGTACCCATCCAAGAGCAAGCTTCTTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTAAATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGACGAGGGCAAAGAAATGATCCTCTCCTATCTGAATAAGGGACCTTTTATTGGAGACCTGGGCCTGTTTGAAGTGGGcccaga                                                                                    >  2:20571/2‑246 (MQ=255)
                                                                                                                                                                gcATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGt                                                                                >  1:564939/1‑251 (MQ=255)
                                                                                                                                                                  aTGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTGAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGt                                                                                >  1:194151/1‑249 (MQ=255)
                                                                                                                                                                              aaaCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGTGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAACGGCTCTGT‑GGCAGTGCTGATCAAAGACGaa                                                                                                                                                                    >  1:247624/1‑153 (MQ=255)
                                                                                                                                                                              aaaCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGaa                                                                                                                                                                    <  2:247624/153‑1 (MQ=255)
                                                                                                                                                                              aaaCCGCAAAACGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCATGAACGTAGCGCATGGGT‑ACg                                                                 <  2:91874/251‑1 (MQ=255)
                                                                                                                                                                               aaCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTTAATAGGTCCAGGAACGTAGCGCATGGGT‑ACGt                                                                <  1:198560/251‑1 (MQ=255)
                                                                                                                                                                                                                tGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGa                                                                                                                                                                               <  1:539628/108‑1 (MQ=255)
                                                                                                                                                                                                                tGTCTCATTGCCACATTCATAAGTAACCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGa                                                                                                                                                                               >  2:539628/1‑108 (MQ=255)
                                                                                                                                                                                                                  tctcATGGCCACATACATAAGACCCCACCCAAGAGCAAGTTTATTCAC‑AGGGTGAAAAAGGCGGAAA‑CCTTGTACAACTCCGTAAAAGGCTCTGT‑GGCAGTGCTGATAAAAGACGAAGAGGGAAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACGGGTCCTGTTTGAAGAGGGCCAGAAACGTAGCGCAGGGGT‑ACGGGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGt                              <  1:390499/250‑1 (MQ=255)
                                                                                                                                                                                                                      aTTGCCACAGTCATAAGTACCCATCCAAGATCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGGGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGGGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGTGCCAGGAAAGTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGCAGTGGCTGAAATTTCGTACaa                          >  2:414918/1‑250 (MQ=255)
                                                                                                                                                                                                                                    taGTACACATCCACGAGCACGCTTATTAAACAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGCAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAAAGTAGCGCATGGGT‑ACGTGCGAAAAACGCCCGTGAAGTGGCTGAAATTTCGGATAAAAAATTTCGCCAAtt           >  2:365056/2‑251 (MQ=255)
                                                                                                                                                                                                                                      gTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTg          >  1:306300/1‑250 (MQ=255)
                                                                                                                                                                                                                                         cccacccaAGAGCAAGCTTATTCACGGGGGTAAAAAAG‑‑GGAAAGCTCTGAACTAAACCGTAAAAGGCTTTGT‑GGCAGTGCTGAAAAAAGAGGAAGAGGGAAAAGAAATGTCCCTCCCTTTTCTGAATCGGGGTGATTTTATTGGGGAATGGGGCTTGTTAGAAGAGGGCAAGGAACGTACCG‑GTGGGTTACGTGAGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTc      <  1:313332/251‑1 (MQ=255)
                                                                                                                                                                                                                                              ccAAGAGCAAGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGG‑GGCAGTGCTGATCAAAGAAGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTCGGCATGTTTGAAGATGGCCAGGAACGTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGACAATTGATTCAggc  >  1:266283/1‑249 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATA‑TTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAG‑CGGAAA‑CGCTGTACTACATCGTTAAAGGCTCTGT‑GGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGT  >  NC_000913/3485959‑3486445

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