Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,052,025 C→A 5.5% S98I (AGC→ATC)  ubiI ← 2‑octaprenylphenol hydroxylase, FAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,052,0250CA 5.5% 241.2 / 3.2 72S98I (AGC→ATC) ubiI2‑octaprenylphenol hydroxylase, FAD‑dependent
Reads supporting (aligned to +/- strand):  ref base C (34/34);  new base A (2/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

GGAAGGTTTCATTTTCTCCCCAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTGGCGGCATTGATAGCCGAAACGC  >  NC_000913/3051879‑3052172
                                                                                                                                                  |                                                                                                                                                   
ggAAGGTTTCATTTTCTCCCCAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCtt                                                                                                                                                   <  5:285861/149‑1 (MQ=255)
 gAAGGTTTCATTTTCTCCCCAGGCGACCTGCTGTTATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCttt                                                                                                                                                  >  8:158152/1‑149 (MQ=255)
     gTTTCATTTTCTCCCCAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATc                                                                                                                                              <  5:173533/149‑1 (MQ=255)
      tttCATTTTCTCCCCAGGCGCCCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCACCGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCg                                                                                                                                             >  5:340505/1‑149 (MQ=255)
      tttCATTTTCTCCCCAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCg                                                                                                                                             >  5:357412/1‑149 (MQ=255)
       ttCATTTTCTCCCCAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGt                                                                                                                                            >  3:192223/1‑149 (MQ=255)
           ttttCTCCCCAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCaaa                                                                                                                                        >  4:157809/1‑149 (MQ=255)
             ttCTCCCCAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAAc                                                                                                                                       >  3:10120/1‑148 (MQ=255)
                  cccAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGaaa                                                                                                                                   <  7:323500/147‑1 (MQ=255)
                  cccAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAAtg                                                                                                                                 >  8:77543/1‑149 (MQ=255)
                  cccAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAAtg                                                                                                                                 >  3:316522/1‑149 (MQ=255)
                  cccAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAAt                                                                                                                                  <  1:19038/148‑1 (MQ=255)
                   ccAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAAtgt                                                                                                                                >  3:17682/1‑149 (MQ=255)
                    cAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAAtgt                                                                                                                                >  6:49956/1‑148 (MQ=255)
                      ggCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATGGCCCATGCTTTGATCGTCAAACGAAATGTGAc                                                                                                                             <  4:10515/149‑1 (MQ=255)
                           ccTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAg                                                                                                                        <  4:17682/149‑1 (MQ=255)
                             tgctgTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCa                                                                                                                                          <  6:149764/129‑1 (MQ=255)
                             tgctgTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCa                                                                                                                                          >  5:149764/1‑129 (MQ=255)
                              gctgTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTg                                                                                                                     >  6:45956/1‑149 (MQ=255)
                               ctgTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTg                                                                                                                     >  6:63511/1‑148 (MQ=255)
                                  tAATTCTGCGGTGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCt                                                                                                                  <  2:49513/148‑1 (MQ=255)
                                     ttCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGt                                                                                                              <  5:45956/149‑1 (MQ=255)
                                     ttCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGt                                                                                                              <  1:217837/149‑1 (MQ=255)
                                       cTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAg                                                                                                                        <  5:259306/137‑1 (MQ=255)
                                       cTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAg                                                                                                                        >  6:259306/1‑137 (MQ=255)
                                       cTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTc                                                                                                             <  5:63511/148‑1 (MQ=255)
                                              ggCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACtt                                                                                                     >  2:76261/1‑149 (MQ=255)
                                                cTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCAAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTc                                                                                                    >  4:61091/1‑148 (MQ=255)
                                                  aaCAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCa                                                                                                  >  3:133996/1‑148 (MQ=255)
                                                   aCAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCAt                                                                                                 <  5:303339/148‑1 (MQ=255)
                                                      gagTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCg                                                                                             <  1:76261/149‑1 (MQ=255)
                                                             atCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATaa                                                                                       >  6:111799/1‑148 (MQ=255)
                                                             atCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATaa                                                                                       >  5:162290/1‑148 (MQ=255)
                                                                  acgacTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACa                                                                                     <  6:162290/145‑1 (MQ=255)
                                                                  acgacTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCt                                                                                  <  4:10120/148‑1 (MQ=255)
                                                                  acgacTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTg                                                                                 <  5:111799/149‑1 (MQ=255)
                                                                   ctacTGATGCGCTTTGTTACACAGCTCGTAGTGAATCACTGAATTTTCAACGATATTCCCAAGATGGCTATAGCCCATGCTTTGATCGTAAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTgg                                                                                <  4:181822/147‑1 (MQ=255)
                                                                    gacTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGc                                                                               <  7:77543/149‑1 (MQ=255)
                                                                         aTGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATATTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCt                                                                            <  2:191854/147‑1 (MQ=255)
                                                                          tGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCAGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTAc                                                                          >  2:334753/1‑148 (MQ=255)
                                                                             gcTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACgaga                                                                      <  4:192223/149‑1 (MQ=255)
                                                                             gcTTTGTTCCACAGCGCGGAGTGCAGCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACgaga                                                                      <  6:340505/149‑1 (MQ=255)
                                                                                   ttCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGGCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGaa                                                                 <  5:49956/148‑1 (MQ=255)
                                                                                    tCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATg                                                               >  1:233444/1‑149 (MQ=255)
                                                                                     ccACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGt                                                              >  7:282630/1‑149 (MQ=255)
                                                                                      cacaGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAAAGTGCCCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTc                                                             >  5:234288/1‑149 (MQ=255)
                                                                                       acaGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGATTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTcc                                                            >  6:13544/1‑149 (MQ=255)
                                                                                       acaGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGATTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTcc                                                            >  4:342713/1‑149 (MQ=255)
                                                                                        caGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCt                                                           <  7:220872/149‑1 (MQ=255)
                                                                                          gcgcGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGc                                                         <  4:316522/149‑1 (MQ=255)
                                                                                                   gAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGcc                                                 >  8:293761/1‑148 (MQ=255)
                                                                                                    aaTCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCaa                                               <  7:158152/149‑1 (MQ=255)
                                                                                                         cTGAATTTTCCACGATATGCCCAAGATGGCTATAGCCCATGCGTTGATCGTCAAACGAAATGTGACCAACGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACgg                                          <  3:31784/149‑1 (MQ=255)
                                                                                                          tGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACggg                                         >  5:32269/1‑149 (MQ=255)
                                                                                                          tGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACggg                                         >  7:59486/1‑149 (MQ=255)
                                                                                                            aaTTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGt                                        <  2:233444/148‑1 (MQ=255)
                                                                                                            aaTTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGt                                        <  6:357412/148‑1 (MQ=255)
                                                                                                              ttttCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAg                                     <  3:157809/149‑1 (MQ=255)
                                                                                                                       atatGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTc                            <  1:334753/149‑1 (MQ=255)
                                                                                                                            cccAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTg                        <  3:61091/148‑1 (MQ=255)
                                                                                                                             ccAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTATGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTggc                      >  8:226531/1‑149 (MQ=255)
                                                                                                                             ccAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTggc                      >  3:39800/1‑149 (MQ=255)
                                                                                                                             ccAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACAGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTggc                      <  4:133996/149‑1 (MQ=255)
                                                                                                                              cAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTggcg                     >  7:46311/1‑149 (MQ=255)
                                                                                                                                  aTGGCTATAGCCCATGATTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTGGCGGCAt                 <  5:13544/149‑1 (MQ=255)
                                                                                                                                  aTGGCTATAGCCCATGATTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTGGCGGCAt                 <  3:342713/149‑1 (MQ=255)
                                                                                                                                            cccATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTGGCGGCATTGATAGCCg        >  1:74170/1‑148 (MQ=255)
                                                                                                                                            cccATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTGGCGGCATTGATAGCCg        >  6:227722/1‑148 (MQ=255)
                                                                                                                                             ccATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTGGCGGCATTGATAGCCGa       >  1:310886/1‑148 (MQ=255)
                                                                                                                                                tGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTGGCGGCATTGATAGCCGAAAcg   >  5:157442/1‑149 (MQ=255)
                                                                                                                                                 gCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAg                                                                                    >  6:8118/1‑67 (MQ=255)
                                                                                                                                                 gCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAg                                                                                    <  5:8118/67‑1 (MQ=255)
                                                                                                                                                  cTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTGGCGGCATTGATAGCCGAAAcgc  <  5:227722/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGAAGGTTTCATTTTCTCCCCAGGCGACCTGCTGTAATTCTGCGGGGGCTAACAGAGTGATATCTGACGACTGATGCGCTTTGTTCCACAGCGCGTAGTGAATCACTGAATTTTCAACGATATGCCCAAGATGGCTATAGCCCATGCTTTGATCGTCAAACGAAATGTGACCAAAGCTGTCTTTGTCCCACACTTCCATACCGTGATAACAGCTGGCCCTACGAGAGAGAATGTCCTGCCAGACGCCAAGACGGGTGAGTAATTTTTCGCTGGCGGCATTGATAGCCGAAACGC  >  NC_000913/3051879‑3052172

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