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,5020CT94.7% 113.2 / ‑0.4 38T128I (ACT→ATT) crpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (25/11);  total (26/12)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-01

AGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTC  >  NC_000913/3486266‑3486734
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                        
aGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTgg                                                                                                                                                                                                                            >  2:676082/1‑251 (MQ=255)
         cAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGCTGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTgg                                                                                                                                                                                                                   >  1:273663/1‑251 (MQ=255)
           aaGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCg                                                                                                                                                                                                                 >  1:310120/1‑251 (MQ=255)
             gACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTAGGCCTGTTTGAAGAGGGCCAGGGACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAACTTTCGTACAAAAAATTTCGCCAATTGTTTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGTTGGCGCGTCGTCTGCAATTAATTTCAGAGAAAGTGGGCAACGTGGCGtt                                                                                                                                                                                                               >  2:189821/1‑251 (MQ=255)
                                     ctctCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGc                                                                                                                                                                                       >  1:494642/1‑251 (MQ=255)
                                            aTCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGGGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGAct                                                                                                                                                                                >  1:560090/1‑251 (MQ=255)
                                                      gggTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCGGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGAGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGACGTTCCTAGACGCGACGGGCCGCATTGCACCAACTCTGCTGAATc                                                                                                                                                                      >  1:435720/1‑251 (MQ=255)
                                                                 ttGGCGAACTGGGCTTGTTGGAAGAGGCCCAGGAACGTAGCGCTTGGGTCGGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAaaca                                                                                                                                                            <  1:306543/250‑1 (MQ=255)
                                                                               cTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGAc                                                                                                                                             >  1:514607/1‑251 (MQ=255)
                                                                                 gtgtGAAGAGGGCCAGAACCGTAGGGCGGGGGTACGGGCAAAACCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACt                                                                                                                                            <  1:300371/247‑1 (MQ=255)
                                                                                               aGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAACAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCAAAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACATGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACAAGACGCTATGACTCACCCGGCCGGTagg                                                                                                                             >  2:381843/1‑249 (MQ=255)
                                                                                                  aaggtagctgtggggcactggaaaaaaCCGCCTGGGAAGTGGCTGAAATTTCGTACAAAAAATTTCTCCAATTGATGCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGTGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCaa                                                                                                                          <  2:634181/228‑1 (MQ=255)
                                                                                                     gTAGCGCAGGGGTACGTGGAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATc                                                                                                                       <  1:801950/251‑1 (MQ=255)
                                                                                                       aGCGCATGGGTAGGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAACCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTAATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATATGGCAAAACAACAAGACGCTATGACTCACCAGGACGGTAGGCAAATCa                                                                                                                      >  2:29348/1‑250 (MQ=255)
                                                                                                                gTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCg                                                                                                            <  1:377987/251‑1 (MQ=255)
                                                                                                                       caacccccGCCTGTGAAGTTGCTGAATTTTCGTACAAAATATTTCGCCAATTGATTCAGGTAATCCCGGACATTCTAATGCGTTTGTCTGCACAGTTGGCGCGTCTTCTGCAAGTCATTTCAGAGAAAGTGGGCACCCTGGCGTTCCTCGACGCGACGGGCCGCATTGCCCAGACTCTGCTGAATCTGTCAACACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGaa                                                                                                     <  2:569965/245‑1 (MQ=255)
                                                                                                                            aCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAACTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTgg                                                                                                <  1:12383/251‑1 (MQ=255)
                                                                                                                              cGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTc                                                                                              >  1:414629/1‑251 (MQ=255)
                                                                                                                                ccTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTc                                                                                              >  1:252886/1‑249 (MQ=255)
                                                                                                                                    ggAAGTGGCTTAAATTTCTTAAAAAAAATTTCGCCAATTGATACAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTGCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTg                                                                                        <  2:298476/250‑1 (MQ=255)
                                                                                                                                                      gTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGt                                                                                                           <  1:342619/214‑1 (MQ=255)
                                                                                                                                                      gTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGt                                                                                                           >  2:342619/1‑214 (MQ=255)
                                                                                                                                                                       aaTTGATTCAGGTAAACCCGGCCATTCTGATGCGTTTGTCTGCTCAGATGGCGCGTCGTCTGCACGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGCCTCTGATGAATCTGGCCAAACAACCAGACGCTATGACGCACCCGCACGGTATGAAAATCAAAATTACCCGTCAGGCAATACGTCAGATTGTCGGCTGTTCGCGGGAAACCGTGGGCCGCATTCTg                                                     >  2:703611/1‑251 (MQ=255)
                                                                                                                                                                             ttCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGCCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGAATTCTGAAGATg                                               >  2:343847/1‑251 (MQ=255)
                                                                                                                                                                              tCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGc                                              >  1:681720/1‑251 (MQ=255)
                                                                                                                                                                                       cccGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGCCGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTAAGATTGTGGGCTGTTCTCGTGAACCCGTGGGCCGCATTCTGAAGATgatgaa                                          >  2:571843/1‑241 (MQ=255)
                                                                                                                                                                                       cccGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCAGCCGGTATGCAAATCAAGATTACACGTCCGGCAATTGGACAGCTTGTCGACTGTTCTCGTGAACCAGAGGGACGCACTCTGaacgttcaggaagatt                                     >  2:720931/1‑237 (MQ=255)
                                                                                                                                                                                            aCATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGACAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAAc                                >  2:511014/1‑251 (MQ=255)
                                                                                                                                                                                            aCATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAAc                                >  2:244056/1‑251 (MQ=255)
                                                                                                                                                                                              aTTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGAAAATTGGTCAGCTTGTCGGCTGTTCTCGTGAAACCGGGGGACGCATTCTGAAGATGCTGGGAGATCAGGACCt                              >  1:774336/1‑251 (MQ=255)
                                                                                                                                                                                                        gTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGACGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAACTTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGc                     >  1:96301/1‑250 (MQ=255)
                                                                                                                                                                                                                     gATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTTCCCGTCAGGAAAGTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCa       <  1:809087/251‑1 (MQ=255)
                                                                                                                                                                                                                      aTGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCAt      <  1:244056/251‑1 (MQ=255)
                                                                                                                                                                                                                         gcgcGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCAtcgt   >  1:10250/1‑251 (MQ=255)
                                                                                                                                                                                                                          cgcgTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCAtcgtc  >  1:803441/1‑251 (MQ=255)
                                                                                                                                                                                                                          cgcgTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACAAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCAtcgtc  >  2:149933/1‑251 (MQ=255)
                                                                                                                                                                                                                             gtcgtcTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGAt                                                <  2:257589/202‑1 (MQ=255)
                                                                                                                                                                                                                             gtcgtcTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGCGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGAt                                                >  1:257589/1‑202 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                        
AGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTC  >  NC_000913/3486266‑3486734

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: ‑