Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,285,238 G→A intergenic (‑25/+175) yjcH ← / ← acs DUF485 domain‑containing inner membrane protein YjcH/acetyl‑CoA synthetase (AMP‑forming)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,285,2380GA95.2% 119.7 / ‑1.9 42intergenic (‑25/+175)yjcH/acsDUF485 domain‑containing inner membrane protein YjcH/acetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (24/16);  total (25/17)
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.42e-01

CCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCG  >  NC_000913/4285005‑4285480
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
ccAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTccc                                                                                                                                                                                                                                   >  1:366880/1‑251 (MQ=255)
 cAACTCCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCATCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCg                                                                                                                                                                                                                                  <  1:318911/251‑1 (MQ=255)
   accccAATTGGAATACCTCGGGTGACGCTGGTTTTCGGATTAAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCt                                                                                                                                                                                                                                <  2:298465/248‑1 (MQ=255)
        aaTTGGAATACCGCGTGGGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCtt                                                                                                                                                                                                                           <  1:28151/251‑1 (MQ=255)
                 cccGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGGGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGaa                                                                                                                                                                                                                  <  1:27906/250‑1 (MQ=255)
                                ggTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCGCgg                                                                                                                                                                                                     <  1:228248/249‑1 (MQ=255)
                                          ttCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAgg                                                                                                                                                                                                                <  1:52198/228‑1 (MQ=255)
                                          ttCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATAAGTAAAATAACGCCGAAATAAACAGCCAGCATAATAATCGACAGGAAGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTATTGTAAGGATTGTTTTAATTCACGCTCCCTTATGGGAGAAgg                                                                                                                                                                                                                >  2:52198/1‑228 (MQ=255)
                                                ggcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGt                                                                                                                                                                                   >  1:244525/1‑251 (MQ=255)
                                                  cgTGCCCAGCCAGCCGTGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAGTCGACAGGATGGTGGCAAACAGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATACGCTGATAAATAGTGCCATTCATCACCGCTTATCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCATTATGGGAGAAGGTTAACGCTCGTGTAACCCTTGCCGAATGt                                                                                                                                                                                   >  2:242703/1‑249 (MQ=255)
                                                           ccagccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATa                                                                                                                                                                         >  1:181151/1‑250 (MQ=255)
                                                                 gggCGCGAAGGCGAGCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGtt                                                                                                                                                                  <  2:181151/251‑1 (MQ=255)
                                                                             gATCAGTAAAATAAAGCCGATATACACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATTCGCATTGTCTTATATCCGCTGATAAATAGTGCCATTCATCACAGCTTCTCCAGAGGTGATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGGAGGCCGGGTAAGGCGATTACGCCGCATcc                                                                                                                                                      >  1:58847/1‑251 (MQ=255)
                                                                                 aGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGGTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGca                                                                                                                                                  <  1:212189/251‑1 (MQ=255)
                                                                                   tAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGAt                                                                                                                                                                                                                                                >  1:344482/1‑155 (MQ=255)
                                                                                   tAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGAt                                                                                                                                                                                                                                                <  2:344482/155‑1 (MQ=255)
                                                                                     aaaTAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACACTTGCCGAATGTAGGCCGGATACGGCGTTTACGCCGCATCAGGcaatc                                                                                                                                               >  1:188803/1‑250 (MQ=255)
                                                                                                     cctctgcctatctataatCGACAGGATGGTGGCAAACCGTTGCCGTGTTTCGACTAACTCCCTGAATTTCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACATATTCTCCAGAGGGAATGGAAGTATTGTTTTAATTCCCGCTCCCTTATTGGATAATGTTAACGCTCCGGTAACCCTTGCCGAATGTAGGCCGGATAAGTGGGTTACGCCGCATCCGGCAATCAATGCCCGATGCGACGc                                                                                                                              <  2:317352/238‑1 (MQ=255)
                                                                                                               ataataATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCTTCt                                                                                                                    >  1:348717/1‑251 (MQ=255)
                                                                                                                taataaTCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCtt                                                                                                                   >  1:61699/1‑251 (MQ=255)
                                                                                                                  ataatCGACAGGATGGTGGCAAACCGTTACCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATAAATGCCTGATGCGACGc                                                                                                                              >  1:317352/1‑238 (MQ=255)
                                                                                                                        gACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCCTTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGCAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTCCGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTCTTCAGGcc                                                                                                           >  2:115681/1‑251 (MQ=255)
                                                                                                                           aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAc                                                                                                        >  1:104779/1‑251 (MQ=255)
                                                                                                                             gATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAAc                                                                                                                                                                                                <  1:60154/161‑1 (MQ=255)
                                                                                                                             gATGGTGGCAAACCGTTGCCGTTTTTAGACTAACTCCCTGAACTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGCGAAGGTTAACGCTCGGGTAAc                                                                                                                                                                                                >  2:60154/1‑161 (MQ=255)
                                                                                                                              aTGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCAAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGAAACCCTTGCCGAATGTAGACCGGATAAGGCGTGTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCATTCGCGTTCTAGCAGGACTACaa                                                                                                      >  2:288260/1‑250 (MQ=255)
                                                                                                                                    gCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTa                                                                                               <  1:147760/251‑1 (MQ=255)
                                                                                                                                         ccGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCgcc                                                                                          >  2:270801/1‑249 (MQ=255)
                                                                                                                                             tGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTTCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCg                                                                                       <  1:214496/250‑1 (MQ=255)
                                                                                                                                                       cGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAAt                                                                            >  1:191090/1‑251 (MQ=255)
                                                                                                                                                          cTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATtgt                                                                         >  1:209357/1‑251 (MQ=255)
                                                                                                                                                                                         ccGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCACTTGTGGGTTACGATGGCATCGCGATAACCTGCTtc                                          >  2:249935/1‑251 (MQ=255)
                                                                                                                                                                                            cTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAc                                                                                                                                                               >  1:285732/1‑131 (MQ=255)
                                                                                                                                                                                            cTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAc                                                                                                                                                               <  2:285732/131‑1 (MQ=255)
                                                                                                                                                                                             tGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCtt                                      >  1:297738/1‑251 (MQ=255)
                                                                                                                                                                                              aaTAAATAGGGCCATTAATCAAAGATTCTCCAGAGGTAATGTAAGTATTGTTTTAATTCCCGCTCCCTTATGGTAGAAGGTTAACGCTCGGTTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGCCGCTGTCGCGTCTTATCAGGCCTTAAAACCGTTACCGTCTCGCCTCGTGCTATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCtta                                     <  2:27523/250‑2 (MQ=255)
                                                                                                                                                                                                tAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                                                 <  1:108271/173‑1 (MQ=255)
                                                                                                                                                                                                tAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                                                 >  2:108271/1‑173 (MQ=255)
                                                                                                                                                                                                  aaTAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGGAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCATACAAACCGTTACCGACTCGCATCGGGCTATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAagc                                 >  2:34737/1‑251 (MQ=255)
                                                                                                                                                                                                          cactcatcaAAGAGTCTCCAGAGGTAATGTAGGGATTGTGTTAATGCCCGCTCCCTTATTGGAGAAGGGTAACGCGCGGGTAACCCTTGCCGAATGTATGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCATGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACTGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCCGCAAGCAGTTtctc                         <  2:163115/248‑1 (MQ=38)
                                                                                                                                                                                                          cATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTtctc                         >  1:374207/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCg  >  1:156672/1‑251 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
CCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCG  >  NC_000913/4285005‑4285480

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