Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,518,627 C→A P34H (CCT→CAT)  insA7 → IS1 family repressor protein InsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,518,6270CA93.2% 118.2 / NA 45P34H (CCT→CAT) insA7IS1 family repressor protein InsA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base A (28/13);  minor base T (1/2);  total (30/15)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGG  >  NC_000913/4518387‑4518862
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
tGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTaa                                                                                                                                                                                                                                      >  1:87401/1‑250 (MQ=255)
     cTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATg                                                                                                                                                                                                                                >  1:369566/1‑251 (MQ=255)
     cTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATg                                                                                                                                                                                                                                >  1:270452/1‑251 (MQ=255)
                     aGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACt                                                                                                                                                                                                              <  1:87935/251‑1 (MQ=255)
                        aTGACTCATGATGTGGTCTGCTATTATTGACATCCTAACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGAGTTTAAGGTGCTTGCGAGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGAGCGTAACGGCAAAAGCACTGCCGGACATCAACGCTATCTCTGCTCTCATTGCCGTACAACATGGCAACTACAGTTCACTTac                                                                                                                                                                                                           >  1:293742/1‑251 (MQ=37)
                          gACTCATGATGTGGTCTGCTATTATTGACATCCTCACGGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCCTCTTACTGCTTTAGTGTATGATGTAGATTTTAAGGCGCTTGCGCGGCTTCCATTTCAATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGTTGCGTAACGGCAGAAGCACTGCCGGACATCGGCGCTATCTCTGCTCTCTTTGCCGTAAAACATGGACACTACTGTTCACTTaccc                                                                                                                                                                                                         >  2:197985/1‑249 (MQ=255)
                                   tgtgGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCtctc                                                                                                                                                                                                 >  1:185920/1‑250 (MQ=255)
                                       gTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGCGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGccc                                                                                                                                                                                            >  1:90394/1‑251 (MQ=255)
                                          tGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGt                                                                                                                                                                                         >  1:65278/1‑251 (MQ=255)
                                                    gACATCCTCACGCCCCAAAAGGATGGGGATGTCGGTAATGCTGCCAACTGCCTGATTTATTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATAAGATTTCCTTCCTTCTCCGCTACTTAAGCCGTGGTGCGTAACGGAAAAAGCATTGCCGGACATCAGCGTTATCCCTCCTCTCATTGCCTTAAAACATGGCAACTACGGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAg                                                                                                                                                                                <  1:164657/250‑1 (MQ=255)
                                                         ccgcACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTc                                                                                                                                                                                                                 <  2:382990/209‑1 (MQ=255)
                                                         ccTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTc                                                                                                                                                                                                                 >  1:382990/1‑212 (MQ=255)
                                                                      aaGGATGGGGATTTCGGTAATGCTGACAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTACGCTACTGAAGGCGTGGTGCGTAACGGCAAACGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCAATTACACCGCCTCTCAGCCCGCGACGCACCAGAAAATCATTGATATGGcc                                                                                                                                                              >  1:305558/1‑250 (MQ=21)
                                                                       aGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGGGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGAAGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACAGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCAt                                                                                                                                                            >  1:252988/1‑251 (MQ=255)
                                                                                    cGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCgg                                                                                                                                                >  1:300932/1‑250 (MQ=255)
                                                                                            ctgcCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTg                                                                                                                                                                                                                                            >  2:294986/1‑152 (MQ=255)
                                                                                            ctgcCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTg                                                                                                                                                                                                                                            <  1:294986/152‑1 (MQ=255)
                                                                                                       cTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTcgcg                                                                                                                                        >  2:110552/1‑239 (MQ=255)
                                                                                                       cTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTcgcg                                                                                                                                        <  1:110552/239‑1 (MQ=255)
                                                                                                                    tgatgGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTgg                                                                                                               >  1:87125/1‑251 (MQ=255)
                                                                                                                     gatgGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCg                                                                                                                                                                                           >  2:60980/1‑174 (MQ=255)
                                                                                                                     gatgGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCg                                                                                                                                                                                           <  1:60980/174‑1 (MQ=255)
                                                                                                                           gATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGGCCTTCCTGCTACGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGAAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCaa                                                                                                         >  1:33325/1‑250 (MQ=255)
                                                                                                                             tattAAAGTGCTTACGTGGCTTCCATTTCCCTCAGTTGTCCTTGCTGATCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTTCAGTACATAATCCCTATCTCTGCTCTGATTGCAGTACACCATGGCACCTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGCAAATCATTGATATGGACATGAATGTCGTCGGATTACGCGCCAGTGCACGCATTATGGGcataggcatctaccc                                                                                                      >  2:106845/3‑237 (MQ=255)
                                                                                                                               ttgtggggtttgtgTGGCTTCCATTTCCATCAGAGTTCCTTCCTGCTCCGCTACGGAAGGCGGGGGGCGTAACGCAAAAAGCACTGCCGGACACCAGCGCTATCGCTGCTCTCATTCCCGTAAAATATGGCAACTACAGTTCATTTACGCCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACg                                                                                                     <  1:345123/237‑1 (MQ=18)
                                                                                                                                tAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTacg                                                                                                                                                                                      >  1:117316/1‑166 (MQ=255)
                                                                                                                                tAAGGTCCTTGGGTGGCTTCCATTTCCACCAGATGTCCTTCTGGCCCCGCTCCTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACACCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTacg                                                                                                                                                                                      <  2:117316/168‑3 (MQ=18)
                                                                                                                                       cTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGt                                                                                            >  1:33105/1‑251 (MQ=255)
                                                                                                                                       cTTGCGTGGCTTCCATTTCCATAAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACTTAAGCGCTATCTCTGCTCTCAATGCCGTAAAACATGGCACCTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATAATTGATATGGACATGCAGGGCCACTGATGTCGCGCCAATGCACGCATTATGGTCGTTGTCCACATCACGGCTTTActt                                                                                            >  2:119749/1‑249 (MQ=255)
                                                                                                                                                   ccatttccatCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAAc                                                                                >  1:292399/1‑251 (MQ=255)
                                                                                                                                                   ccatttccatCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACAACAGTTCACTTACACCGCCTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCACGCCAGAGACCGAATTATGGGCGTTTGCCTCAACAAGGTGTTACGTCAATTaaaaaa                                                                                 >  2:351027/1‑250 (MQ=17)
                                                                                                                                                     atttccatCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCCGAAAATCCTTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTc                                                                              >  1:147262/1‑251 (MQ=255)
                                                                                                                                                               gATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCTCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGt                                                                    >  1:89836/1‑251 (MQ=255)
                                                                                                                                                                aTGTCCTTCCTGCCCCGCTACTGAAGGGGTGGTGCGTAACGGCAAAAGCCCTGCCGGCCATCAGCGCTATCTCTGCCCTCATTGCCGTAAAACAGGGCAACTACAGTTCCCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTc                                                                   <  1:196850/251‑1 (MQ=17)
                                                                                                                                                                 gttccttcctGCCCCGCTCCTGAAGGGGTGGTGGTTAACGGCAAAACCCCTGCGGACCACCAGCGCCTCCCCTCCTCTATTGGCCTTAAAACATGGCAACTACCGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCg                                                                  <  2:87401/249‑1 (MQ=17)
                                                                                                                                                                                   ggcgAAGGCGTGGTGCGTAGCGTCAAAATCTCGGCCTAACATCAGCGCTACCTCTGTTCCTTTTGCCGTAAAACATGGCAATATCATTTCACTTACACCGCCTCTCAGCCCGGTACGCGCCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACTCGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACaa                                                <  2:219962/248‑1 (MQ=17)
                                                                                                                                                                                   aCTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTCCACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTGGGCCTCAACACGGTTTTACGTCACTTCAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACa                                                 >  1:175323/1‑250 (MQ=255)
                                                                                                                                                                                    cTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTAATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACAAGGTTTTCCGTCACTTAAAAAACTCAGGCCGGAGCCGGTAACCTCGCGCATACAAc                                               >  2:328472/1‑251 (MQ=17)
                                                                                                                                                                                      acaGGCGTGGGGCGTAACGTCAAAACTCGTGCCGGACATCGGCGCTATCCCTGCCCTCATTGCCTAAAAACTTGGAACCTCCGTGTCCCTGCCACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCg                                             <  1:105729/249‑1 (MQ=17)
                                                                                                                                                                                        aGGCGTGGTACGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGAACGCATTATGGGCGTTGGCCTCAACACGGATTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCggg                                           >  2:103567/1‑251 (MQ=255)
                                                                                                                                                                                           cGTGGCGCGTAACGGCAAAAGCCAGGCCGTACATCAGCGCTATCTCTGCGCCCATGGCCGAAAAACAGGGCAACAACAGTTCATTGACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTGAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAg                                        <  2:185920/251‑1 (MQ=17)
                                                                                                                                                                                                          cAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTg                         >  1:162458/1‑251 (MQ=255)
                                                                                                                                                                                                           aaaaGCACTGCCGGACACCAGCGCTATCTCGGCTCTCATTGCCTTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGTCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACTGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTgc                        <  1:16594/251‑1 (MQ=18)
                                                                                                                                                                                                                     gCCGGACATCAGCGCTACCTCTGCTCCCATTGCCGTAAAACAGGGCAACTACAGTCCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGaaa                 <  1:73632/249‑1 (MQ=18)
                                                                                                                                                                                                                                   gcTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTggg  >  1:281601/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
TGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGG  >  NC_000913/4518387‑4518862

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