Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,560,309 T→A P124P (CCA→CCT yjiH ← uncharacterized protein YjiH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,560,3090TA91.7% 98.9 / ‑1.7 36P124P (CCA→CCTyjiHuncharacterized protein YjiH
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (14/19);  total (15/21)
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.22e-01

GACAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGC  >  NC_000913/4560069‑4560552
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
gACAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCca                                                                                                                                                                                                                                            <  1:168147/251‑1 (MQ=255)
         aaaCACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGATCAGATCCCAGGAATGCCCAACAACGGt                                                                                                                                                                                                                                    <  2:322594/250‑1 (MQ=255)
                  gATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGcac                                                                                                                                                                                                                            >  1:44010/1‑249 (MQ=255)
                     gatggCTAAAGGTACGATTACGGACGTGCCCAGAACGGCAAACAACGCGACGCCGGACGAGAAGGAGTTGGTGATGATGGCGCTACCGCTGGATTGATATGCGGCGAAAATCACTTTATCGCGATCGGGAATTTCATCTTCCTGCGCCAGCTCTTTGGACATACCGGAGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCAAGCGAGCAGATCCCAGGGATGGCCAACAACGATTTTAGCACAGgtg                                                                                                                                                                                                                       >  1:400330/4‑251 (MQ=255)
                             aaGGTACGATTACGGACGGGCCCAGAAAGGCAAACCCCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAAc                                                                                                                                                                                                              <  1:334647/251‑1 (MQ=255)
                                              gTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTaa                                                                                                                                                                                             <  1:291977/251‑1 (MQ=255)
                                                            aaaCACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCAt                                                                                                                                                                               <  1:147949/251‑1 (MQ=255)
                                                                  cgcgACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGAGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTa                                                                                                                                                                          >  1:200784/1‑250 (MQ=255)
                                                                      aCGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGAGTTTTGCAGGTTAGCGATCATCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTTTCATTAACTGCTGCGCGGCGCGTAATCAACCAAGACAACCAGTAATTg                                                                                                                                                                      >  2:110669/1‑250 (MQ=255)
                                                                            gAAGAGATGTAGTTGTTGATGAGGGCGCTTCCGCTTGTTGGATATGCGGCGTAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTTTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCCTCGGTAATGGAaataat                                                                                                                                                               <  2:21076/251‑1 (MQ=255)
                                                                               gagaAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATccc                                                                                                                                                            <  1:274467/251‑1 (MQ=255)
                                                                                      agtTGGTGATGATGTCGCTACCGCTGGTTTGATATGCGGCGAAAATCCCTTTTTCGCGGTCGGTAATTTCTCCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGTAGGTTAGCTATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCGGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCagaga                                                                                                                                                      <  1:37210/250‑1 (MQ=255)
                                                                                           gTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCTCTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAAt                                                                                                                                                <  1:387132/251‑1 (MQ=255)
                                                                                                gatgGCGCTACCGCTGGTTTGATATGCGGCGAAAATCCCTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAAtgactg                                                                                                                                           <  2:143855/251‑1 (MQ=255)
                                                                                                   ggCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCCCAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACg                                                                                                                                        >  2:271283/1‑251 (MQ=255)
                                                                                                        tACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGAACATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGcc                                                                                                                                   >  1:131614/1‑251 (MQ=255)
                                                                                                        tACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGAGGCATCGGTGTTTTGCCGGTTAGCGATCCGCGCCAGCCAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGAGCAGCGCGTAATCCACCAAAACCATCGGTAATTGAAATAATCACCAGAGAGAGAATGACTGACggagac                                                                                                                                   >  1:484564/1‑247 (MQ=255)
                                                                                                          ccGCTGGTTTGATATGCGGCGAAAATCGCTTTATCGCGTTCGGTAATTTCGCCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCAGCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCTCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCACCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAg                                                                                                                                 <  2:344032/251‑1 (MQ=255)
                                                                                                          ccGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGAAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGGGTAATCCACCAAGACCATCGGTAATTGAAATAAACCCCAGAGAGAGAATGACTGACGGTGcc                                                                                                                                   >  1:261456/1‑249 (MQ=255)
                                                                                                                  tcgATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAgcg                                                                                                                         <  2:152169/248‑1 (MQ=255)
                                                                                                                                      cTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGt                                                                                                                                                                           <  1:108932/181‑1 (MQ=255)
                                                                                                                                      cTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGt                                                                                                                                                                           >  2:108932/1‑181 (MQ=255)
                                                                                                                                               gTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCcgcg                                                                                            >  1:414866/1‑251 (MQ=255)
                                                                                                                                                             ccTTCATGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGAGCAGCGCCAGCGGGCAGATGCCTGGAATGACCAACAACGGTTTTAGCACAGGTGTCATTAACAGCTGCGCGGCGCGTAATACACCAAGACCATCGGTAGTTGAAGTAATCCCCAGAGAGAGAATGACTACCGGTGCTGGCTCCAGCGCACACAGAAAGCCATCTTTAGCCCCCGCGCCACccgcaccgcg                                                                              >  1:76701/1‑251 (MQ=255)
                                                                                                                                                                     cgcCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAgg                                                                      >  1:353642/1‑251 (MQ=255)
                                                                                                                                                                               ttGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTc                                                             <  1:271283/250‑1 (MQ=255)
                                                                                                                                                                                               cggccgctgagtTTTGCATGTTGGCGATCATCGCCTGCGATCAGATCCCAGGAATGCTCAACCACGGTTGTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCTGAGAGAGAATGACTGACGGTGCCACCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCTCCCGCGCCGCGAAATGTGGTGGTTGCTCCGTTCGCATCAGGCAGt                                            <  2:331864/241‑1 (MQ=255)
                                                                                                                                                                                                         ggtttGCAGGTGATCGATCAGCGACAGTGAGCATGTCCCAGGAAAGCCCAACAACGGGTTTAGCACATGTGTCATTAACTGCTGCGCGGCGCGTAACCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGTCAGTTGCCCAAAct                                  <  2:348010/249‑2 (MQ=255)
                                                                                                                                                                                                            tggcatgTTAGCGATCAGCGAAAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGTGTGCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCAACGGGAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGGTCCAGCTCCAGCGCAACGAGAAAGCCATCTTTAGCCCCTGCGCCGCCCTCGCCGCGAAATGAGTGGGAGGCTCCGTTCGCCCAAGGCAGTTGCCCAAACGAt                                <  2:246146/244‑1 (MQ=255)
                                                                                                                                                                                                                  gTTATCGATCAGCGCCAGCGAGCCGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCTCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCa                         <  2:466184/251‑1 (MQ=255)
                                                                                                                                                                                                                    tAGCGATCAGCGCCAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAGGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGGGGTGGCTCCGTTCGCACCATGCAGTTGCCCAACCGATCCGTTc                          >  1:384065/1‑248 (MQ=255)
                                                                                                                                                                                                                         aTCAGCGCCAGTGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGCAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGAGGTGGCTCCGTTCGCACCAGGCAGTTGACCaacgcacccgttcaaaaagg                   >  2:328479/1‑232 (MQ=255)
                                                                                                                                                                                                                           cAGCGCCAGCGAGCAGATCCCAGGAATACCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGaaa                <  1:344183/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccAGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCa           <  1:435264/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aGCGAGCAGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCaaa         <  1:49170/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aGATCCCAGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACgcgc  >  1:372150/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
GACAAATTTAAACACCAGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGC  >  NC_000913/4560069‑4560552

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