Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,486,502 C→T T128I (ACT→ATT)  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,5020CT86.7% 114.1 / 13.0 45T128I (ACT→ATT) crpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base T (23/16);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 9.16e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-01

AGTGCTGATCAAAGACGAAGAGG‑GTAAAGAAATGATCCTCTCCTATCTG‑AATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGG‑CCAGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACT  >  NC_000913/3486266‑3486746
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
aGTGCTGATCAAAGACGAAGAGG‑GTAAAGAAATGATCCTCTCCTATCTG‑AATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGG‑CCAGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTGCAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAgagg                                                                                                                                                                                                                                        >  1:320550/1‑248 (MQ=255)
 gggCTGATAAAAGACGA‑GTGGTGTAAAAAAATGATCCTCCCCT‑TTTGTATTCAGGGTTATTTTTTTGGCGAACGGGGCCTGTTTGAAGAGGC‑CCAGGAACTT‑AGCGCTTGGGTCCGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTggg                                                                                                                                                                                                                                       <  1:190874/249‑1 (MQ=255)
           aaGACAAAGAGG‑GAAAAGAAATGATCCTCTCAATTCTA‑AATCAGGGTGTTTTTATTGGAGAAGTGG‑CCTTGTTGAAGGGGGGCCAAAAACGT‑AGCGCAGGGGTACGTACAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCg                                                                                                                                                                                                                             <  1:184250/251‑1 (MQ=255)
                   gagG‑GTAAAGAAATGATCCTCTCCTATCTG‑AATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGG‑CCAGGAACGT‑AGCGCATGGGTACGTGCGAAACCCGCCTGTGAAGTGGATGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAACCCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGCAAGTGGGCAACCTGGCGTTCCTCgc                                                                                                                                                                                                                     >  1:230656/1‑250 (MQ=255)
                   gagG‑GTAAAGAAAGGATCCGTGCATATCTG‑AGGCAGGGTGATTTTGTTGGCGAGCTGGCCTTGTTTGAAGGGGG‑CCAGAAACGT‑AGCGCATGGAGACGTGCAAAAACCGCCTGTGAAGTTGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGa                                                                                                                                                                                                                     <  2:151598/251‑1 (MQ=255)
                            aaaaaTGATCCCCTCTTATCGG‑AATAAGGGTGATTTTATGGGCGAATGGGGCCTTTTTGAAGAGGG‑CCAGAAACGT‑GGCGCAGGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACgg                                                                                                                                                                                                             <  1:9639/249‑1 (MQ=255)
                             gAAATGATCCTCTCCTATCTG‑AATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGG‑CCAGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCATGTGAAGTGGCTGAACTTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACggg                                                                                                                                                                                                            >  1:86559/1‑251 (MQ=255)
                                 tGATCCTCTCCTATCTG‑AATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGG‑CCAGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGc                                                                                                                                                                                                        >  1:356343/1‑251 (MQ=255)
                                  gATCCTCTCCTATCTG‑AATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGG‑CCAGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGc                                                                                                                                                                                                        >  2:481960/1‑250 (MQ=255)
                                  gATCCTCTCCTATCTG‑AATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGG‑CCAGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGc                                                                                                                                                                                                        >  2:362097/1‑250 (MQ=255)
                                       tctcCTATCTG‑AATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGG‑CCAGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGc                                                                                                                                                                                                   >  2:111462/1‑250 (MQ=255)
                                       tccccTATCTG‑AATCAGGTGGATTTTTTGGGCGAAGTGGCCTTGTTTGAAGGGGC‑CCAGGAA‑GTGAGCGAAGGGGGACGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTg                                                                                                                                                                                                    <  2:30735/246‑1 (MQ=255)
                                                    aTCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGG‑CCAGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTctgctg                                                                                                                                                                                      >  1:140320/1‑251 (MQ=255)
                                                           tGATTTTATTGGAGAACTGGGCCTGTTTGAAGAGGG‑CCAGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTgg                                                                                                                                                                               >  2:226818/1‑251 (MQ=255)
                                                                     ggCGAACTGGGCCTGTTTGAAGAGGG‑CCAGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAAcc                                                                                                                                                                     >  1:430699/1‑251 (MQ=255)
                                                                          aCTGGGCGTGTTTGAAGAGGG‑CCAGAAACGT‑AGCGCTTGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACg                                                                                                                                                                <  1:123386/251‑1 (MQ=255)
                                                                              ggCCTTTTTGAAGAGGG‑CCAGGAACGG‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTAt                                                                                                                                                            <  1:293604/251‑1 (MQ=255)
                                                                                              gacagggaccGT‑AGCGGAGGGGTACTTGCAAAAACCGCCTGTGAAGTGGCTGAAATTGCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACgg                                                                                                                                             <  2:311514/242‑1 (MQ=255)
                                                                                                 cAGGAACGT‑AGCGCATGGGTCCGTGGAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTAt                                                                                                                                          <  1:12178/251‑1 (MQ=255)
                                                                                                 cAGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAg                                                                                                                                                                                               >  1:455700/1‑198 (MQ=255)
                                                                                                  aGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATg                                                                                                                                         >  1:490150/1‑251 (MQ=255)
                                                                                                    gAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGATTCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGc                                                                                                                                        >  1:202464/1‑250 (MQ=255)
                                                                                                                    gTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCCGACGCTATGACTCACCCGGACGGTATGAAAATCAAAATTAcgcg                                                                                                                        >  1:300080/1‑248 (MQ=255)
                                                                                                                      aagtgaaaaaaCCGCATGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTc                                                                                                                      <  2:33346/244‑1 (MQ=255)
                                                                                                                         ggcaaaaaCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAACCCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCTGGACTGTATGCAAATCAAAATTTCCCGTcc                                                                                                                     <  1:491726/245‑2 (MQ=255)
                                                                                                                                   gCCTGTGAAGTGGCTGAGAGTGCGATGAAAAAATTTCTACAAGTGATTCAGGTAAACCCTGATGGTCTGATGCGTTTGTCTGCACAGATGGCGTGTCTTCTGCAAGTAATTTCAGAGAAAGGGGGCAACCGGGCGTTCCTCGACGTTACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCGATGACTCACCCGGACGGGATGCAAAGCAAAATTACCCGTCAGGAAATTGGTCa                                                                                                         <  2:49494/251‑1 (MQ=255)
                                                                                                                                              ggCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCt                                                                                              <  1:33628/251‑1 (MQ=255)
                                                                                                                                              ggCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCt                                                                                              <  1:133214/251‑1 (MQ=255)
                                                                                                                                                  gAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGt                                                                                            <  1:471819/249‑1 (MQ=255)
                                                                                                                                                            aCAAAAAATTTCGCCAATTGATTCATGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTTTGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAAc                                                                                 <  1:299421/250‑1 (MQ=255)
                                                                                                                                                                          cAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCt                                                                  >  2:216888/1‑251 (MQ=255)
                                                                                                                                                                                 ttCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATg                                                           >  2:351251/1‑251 (MQ=255)
                                                                                                                                                                                       tAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGCCTCTGCTGAATCTGGCACAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGa                                                      >  2:177047/1‑250 (MQ=255)
                                                                                                                                                                                            ccGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCa                                                <  2:490150/251‑1 (MQ=255)
                                                                                                                                                                                              ggACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGAt                                                            >  2:29255/1‑237 (MQ=255)
                                                                                                                                                                                              ggACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGAt                                                            <  1:29255/237‑1 (MQ=255)
                                                                                                                                                                                                aCATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGGTATCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAAc                                            >  1:126999/1‑251 (MQ=255)
                                                                                                                                                                                                   ttCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTg                                         >  2:133574/1‑251 (MQ=255)
                                                                                                                                                                                                       taTGCGTTGGTCAGCACAGATGGCGCGTCGTCTGTAATTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCGTTGCACAGACTCTGCTTAAGCTGGCAAAACAACCAGACGCGATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTGAGATTGTCGGCTGTTGTCGTGAAACCGTGGGACGCATTCTGATGATGCTGGAAGAATAGAACCTGAtc                                      <  2:230656/249‑1 (MQ=255)
                                                                                                                                                                                                        aTGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGAtctc                                    <  1:226818/251‑1 (MQ=255)
                                                                                                                                                                                                               tGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGcaca                              <  1:249898/250‑1 (MQ=255)
                                                                                                                                                                                                                             gcgcGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCAtcgt               <  2:87464/251‑1 (MQ=255)
                                                                                                                                                                                                                               gcgTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCAtcgtcg             >  1:44783/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGc     >  2:194752/1‑250 (MQ=255)
                                                                                                                                                                                                                                          aaGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACt  >  2:13551/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
AGTGCTGATCAAAGACGAAGAGG‑GTAAAGAAATGATCCTCTCCTATCTG‑AATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGG‑CCAGGAACGT‑AGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACT  >  NC_000913/3486266‑3486746

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