Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,052,785 A→T 5.2% intergenic (+17/+33) prmB → / ← PP_1828 N5‑glutamine S‑adenosyl‑L‑methionine‑dependent methyltransferase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,052,7850AT 5.2% 243.3 / 2.4 77intergenic (+17/+33)prmB/PP_1828N5‑glutamine S‑adenosyl‑L‑methionine‑dependent methyltransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base T (2/2);  total (39/38)
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.33e-01

CATGTCGAGGCGTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTACGCGCT  >  NC_002947/2052640‑2052930
                                                                                                                                                 |                                                                                                                                                  
cATGGCGAGGCGTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGc                                                                                                                                                 <  2:115686/149‑1 (MQ=255)
 aTGTCGAGGCGTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCg                                                                                                                                                >  2:59728/1‑149 (MQ=255)
    tCGAGGCGTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGcc                                                                                                                                             >  5:373439/1‑149 (MQ=255)
      gAGGCGTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCgg                                                                                                                                               >  3:247236/1‑145 (MQ=255)
      gAGGCGTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGccc                                                                                                                                            <  6:373439/148‑1 (MQ=255)
      gAGGCGGTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCgg                                                                                                                                               <  4:247236/145‑1 (MQ=255)
         gcattGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCTTCGCGGCCCCaaa                                                                                                                                        <  5:162411/146‑2 (MQ=255)
            ttGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCTTCGCGGCCCAAAtagc                                                                                                                                     <  4:27917/149‑4 (MQ=255)
             tGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTcc                                                                                                                                    <  3:291390/149‑1 (MQ=255)
              gTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTccc                                                                                                                                   >  8:42583/1‑149 (MQ=255)
               tACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTcccc                                                                                                                                  <  4:80336/149‑1 (MQ=255)
                 ccccGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCg                                                                                                                                >  5:275316/1‑149 (MQ=255)
                  cccGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCgg                                                                                                                               >  6:82175/1‑149 (MQ=255)
                  cccGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCgg                                                                                                                               >  1:300657/1‑149 (MQ=255)
                  cccGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCAGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCgg                                                                                                                               >  5:144666/1‑149 (MQ=255)
                    cGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCt                                                                                                                             >  3:323401/1‑149 (MQ=255)
                           gACTTTGGCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATTCg                                                                                                                       <  5:6102/148‑1 (MQ=255)
                              tttGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCt                                                                                                                   >  4:104585/1‑149 (MQ=255)
                                tGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCt                                                                                                                   <  8:80584/147‑1 (MQ=255)
                                     ggctggAGGCCGAGCCGCCCCGGCAGCGGGTGGCCATGCTCACCGCCGAGCAGTGCCGCCAGCATCCGGAAGACCCCACGGCCCTCCGTTAGGCTATTCCGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCCGCCGAGGTCTTCGCggg                                                                                                            <  4:115554/149‑1 (MQ=255)
                                      gctggAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCgg                                                                                                             <  1:59728/147‑1 (MQ=255)
                                          gAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCgg                                                                                                             >  3:205879/1‑143 (MQ=255)
                                          gAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCgg                                                                                                             <  4:205879/143‑1 (MQ=255)
                                          gAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAg                                                                                                         <  3:80327/147‑1 (MQ=255)
                                             ttCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATc                                                                                                    <  4:323401/149‑1 (MQ=255)
                                                 aGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGa                                                                                                <  5:82175/149‑1 (MQ=255)
                                                     gggcgggcATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCaa                                                                                            <  7:42583/149‑1 (MQ=255)
                                                       gcgggcATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTATTTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACa                                                                                          <  7:63104/149‑1 (MQ=255)
                                                           gcATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGccc                                                                                      >  4:340894/1‑149 (MQ=255)
                                                            cATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATcc                                                                                                   <  6:127497/135‑1 (MQ=255)
                                                            cATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATcc                                                                                                   >  5:127497/1‑135 (MQ=255)
                                                            cATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCcct                                                                                     <  6:275316/149‑1 (MQ=255)
                                                              tGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCcctgc                                                                                   <  4:149992/149‑1 (MQ=255)
                                                                     ttCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGaa                                                                             >  2:206090/1‑148 (MQ=255)
                                                                     ttCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGaa                                                                             >  3:292950/1‑148 (MQ=255)
                                                                        aTGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCcctgcctg                                                                                >  1:36491/1‑142 (MQ=255)
                                                                        aTGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCcctgcctg                                                                                <  2:36491/142‑1 (MQ=255)
                                                                        aTGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAacac                                                                          <  3:104585/148‑1 (MQ=255)
                                                                        aTGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCAGGAACACg                                                                         >  8:225253/1‑149 (MQ=255)
                                                                         tGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCTCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACgg                                                                        >  6:40809/1‑149 (MQ=255)
                                                                         tGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACgg                                                                        <  6:176305/149‑1 (MQ=255)
                                                                            tCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGATGATGCGGTATTAGCGGGGAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCa                                                                      >  7:67074/1‑148 (MQ=255)
                                                                              aCCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAAc                                                                   >  5:132469/1‑149 (MQ=255)
                                                                              aCCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAAc                                                                   <  6:144666/149‑1 (MQ=255)
                                                                              aCCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGGCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAAc                                                                   >  6:96419/1‑149 (MQ=255)
                                                                                cgccgAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGc                                                                 >  1:233728/1‑149 (MQ=255)
                                                                                  ccgAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCcctgcct                                                                                 <  1:289993/131‑1 (MQ=255)
                                                                                  ccgAGCACTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGAGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCcctgcct                                                                                 >  2:289993/1‑131 (MQ=255)
                                                                                   cgAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGccac                                                             >  2:73366/1‑149 (MQ=255)
                                                                                           gccgccAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACcaggcag                                                      <  3:61122/148‑1 (MQ=255)
                                                                                            ccgccAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGt                                                    <  2:233728/149‑1 (MQ=255)
                                                                                            ccgccAGCATCAGGAATTGTTCAAGGCCCTAGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACcaggcagg                                                     <  7:346696/148‑1 (MQ=255)
                                                                                                cAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCTCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATg                                                <  5:40809/149‑1 (MQ=255)
                                                                                                    aTCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTg                                             >  8:118014/1‑148 (MQ=255)
                                                                                                    aTCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTg                                             >  5:214209/1‑148 (MQ=255)
                                                                                                     tCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGa                                            >  4:86353/1‑148 (MQ=255)
                                                                                                       aGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAAc                                         >  6:253926/1‑149 (MQ=255)
                                                                                                          aaTTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACg                                        <  7:35202/147‑1 (MQ=255)
                                                                                                          aaTTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACg                                        >  8:35202/1‑147 (MQ=255)
                                                                                                                  aaGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATccgccg                              >  6:94851/1‑149 (MQ=255)
                                                                                                                   aGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGt                             <  3:86353/149‑1 (MQ=255)
                                                                                                                   aGGCCCTCGTTTAGGCTATGTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGt                             <  2:86873/149‑1 (MQ=255)
                                                                                                                       ccTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGt                                              >  8:300822/1‑128 (MQ=255)
                                                                                                                       ccTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGt                                              <  7:300822/128‑1 (MQ=255)
                                                                                                                           gTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCgcgc                     >  8:126825/1‑149 (MQ=255)
                                                                                                                               aGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTg                   <  1:73366/148‑1 (MQ=255)
                                                                                                                               aGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTgg                  <  5:96419/148‑1 (MQ=255)
                                                                                                                                  gaaattagggCCGCTTCGCGGCCCCAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCATCGCGCTCGTAtt              >  5:293160/8‑149 (MQ=255)
                                                                                                                                   tATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATccgcc                               >  4:75558/1‑131 (MQ=255)
                                                                                                                                   tATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATccgcc                               <  3:75558/131‑1 (MQ=255)
                                                                                                                                   tATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTAttt             >  8:220112/1‑149 (MQ=255)
                                                                                                                                     tttGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTAttt             <  7:220112/147‑1 (MQ=255)
                                                                                                                                     tttGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTAttt             <  7:118014/147‑1 (MQ=255)
                                                                                                                                        agggCCGCTTCGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGAGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTa        >  3:351571/2‑149 (MQ=255)
                                                                                                                                           gCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTAcgc     >  5:127181/1‑149 (MQ=255)
                                                                                                                                              gcgAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTACGCGCt  >  6:74420/1‑149 (MQ=255)
                                                                                                                                              gcgAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTACGCGCt  >  4:279118/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CATGTCGAGGCGTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAAACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTACGCGCT  >  NC_002947/2052640‑2052930

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