Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,052,786 A→C 6.0% intergenic (+18/+32) 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,7860AC 6.0% 250.1 / 4.1 82intergenic (+18/+32)prmB/PP_1828N5‑glutamine S‑adenosyl‑L‑methionine‑dependent methyltransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/39);  new base C (3/2);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGGCGTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTACGCGCTCA  >  NC_002947/2052644‑2052932
                                                                                                                                                 |                                                                                                                                                    
tCGAGGCGTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCTTCGCGGcc                                                                                                                                                >  1:425559/1‑149 (MQ=255)
gcGATGCGTTGTACCCCGAGGGGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTGATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGcc                                                                                                                                                <  2:115484/148‑1 (MQ=255)
   aGGCGTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCAATGGGGTGTGCTGGTTACCGGCGGCGAGTTCCGGCAGGATGATGGAATTGTGAAGGCCCGGCGTTTAGGGTTTTTGGGGCCGCAAGCCGGcct                                                                                                                                               >  1:336400/1‑147 (MQ=255)
      cGTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAAtt                                                                                                                                          <  2:267869/149‑1 (MQ=255)
       gTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCTTCGCGGCCCCaaata                                                                                                                                         <  2:425559/149‑4 (MQ=255)
       gTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAAttt                                                                                                                                         >  2:200486/1‑149 (MQ=255)
        ttGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTc                                                                                                                                        <  2:126647/149‑1 (MQ=255)
        ttGTACCCCGAGGTGGACTATGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTc                                                                                                                                        <  5:182234/149‑1 (MQ=255)
          gTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTAGGGCCGCTTCGCGGCCCCaaatagc                                                                                                                                       <  7:100317/148‑6 (MQ=255)
            aCCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTccccc                                                                                                                                    >  7:154902/1‑149 (MQ=255)
               ccGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGc                                                                                                                                 <  1:201756/149‑1 (MQ=255)
               ccGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGc                                                                                                                                 >  8:23457/1‑149 (MQ=255)
                cGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCt                                                                                                                                <  5:21744/149‑1 (MQ=255)
                    gTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATg                                                                                                                            <  5:358442/149‑1 (MQ=255)
                       gACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCg                                                                                                                          <  1:200486/148‑1 (MQ=255)
                         cTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTc                                                                                                                       >  6:83435/1‑149 (MQ=255)
                              ccTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGc                                                                                                                  <  1:239924/149‑1 (MQ=255)
                              ccTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGc                                                                                                                  <  8:154902/149‑1 (MQ=255)
                                tggctggAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTAGGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCgg                                                                                                                >  4:411124/1‑149 (MQ=255)
                                    tggAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAg                                                                                                            >  6:196110/1‑149 (MQ=255)
                                    tggAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCACTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAg                                                                                                            >  6:80202/1‑149 (MQ=255)
                                     ggAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAg                                                                                                            >  1:378288/1‑148 (MQ=255)
                                      gAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAg                                                                                                            <  1:219918/147‑1 (MQ=255)
                                          tCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATc                                                                                                       >  8:100367/1‑148 (MQ=255)
                                            gAGCGGGGCGGGCATGGGGTGTTCATGCTGACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCa                                                                                                     <  7:14213/148‑1 (MQ=255)
                                            gAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATc                                                                                                       >  8:14213/1‑146 (MQ=255)
                                              gcggggcgggCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGAt                                                                                                  <  5:398583/149‑1 (MQ=255)
                                                   gcgggcATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTAGGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACa                                                                                             <  3:411124/149‑1 (MQ=255)
                                                    cgggcATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAg                                                                                            <  5:196110/149‑1 (MQ=255)
                                                             gggTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCACATCAACAGCcctgcct                                                                                    >  1:394268/1‑149 (MQ=255)
                                                               gtgtTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTgg                                                                                  <  2:336400/149‑1 (MQ=255)
                                                                 gtTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCATAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGaa                                                                                >  6:135440/1‑149 (MQ=255)
                                                                    cATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAaca                                                                              <  6:185911/148‑1 (MQ=255)
                                                                       gCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACgg                                                                           >  8:227519/1‑148 (MQ=255)
                                                                          cACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCaa‑‑a                                                                     >  4:54613/1‑149 (MQ=255)
                                                                           aCCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCaa‑‑a                                                                     <  2:378288/148‑1 (MQ=255)
                                                                             cgccgAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACg                                                                            >  2:57935/1‑141 (MQ=255)
                                                                             cgccgAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACg                                                                            <  1:57935/141‑1 (MQ=255)
                                                                               ccgAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGcca                                                                <  2:394268/149‑1 (MQ=255)
                                                                                cgAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGccac                                                               >  3:68078/1‑149 (MQ=255)
                                                                                      gTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACa                                                                                             <  7:239397/115‑1 (MQ=255)
                                                                                      gTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACa                                                                                             >  8:239397/1‑115 (MQ=255)
                                                                                      gTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACcaggca                                                         <  2:208253/149‑1 (MQ=255)
                                                                                      gTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGCCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGCCTGATGTGGTCTAAGCGGGTAGCAATCCACATCAACACCCCTCCCTGCAACCCGGCAA‑‑ACGACAAcaggca                                                         >  2:125982/1‑149 (MQ=255)
                                                                                           gccAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATg                                                                                                                            >  7:74616/1‑79 (MQ=255)
                                                                                           gccAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATg                                                                                                                            <  8:74616/79‑1 (MQ=255)
                                                                                               gCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGt                                                >  5:211207/1‑149 (MQ=255)
                                                                                                  tCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGa                                              >  8:319431/1‑148 (MQ=255)
                                                                                                    aGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGaaa                                            <  6:211207/148‑1 (MQ=255)
                                                                                                    aGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAAc                                           >  3:213629/1‑149 (MQ=255)
                                                                                                      ggAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACg                                          >  8:95507/1‑149 (MQ=255)
                                                                                                      ggAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACg                                          >  6:390527/1‑149 (MQ=255)
                                                                                                        aaTTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGc                                         >  5:421572/1‑148 (MQ=255)
                                                                                                          ttGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAAt                                      >  4:118951/1‑149 (MQ=255)
                                                                                                            gTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATc                                     <  5:83435/148‑1 (MQ=255)
                                                                                                             ttCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATccg                                   >  4:333893/1‑149 (MQ=255)
                                                                                                              tCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATccgc                                  >  7:188516/1‑149 (MQ=255)
                                                                                                                aaGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATccgccg                                <  7:100367/149‑1 (MQ=255)
                                                                                                                aaGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATccgcc                                 <  7:23457/148‑1 (MQ=255)
                                                                                                                     ccTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGaa                                             >  2:210641/1‑131 (MQ=255)
                                                                                                                     ccTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGaa                                             <  1:210641/131‑1 (MQ=255)
                                                                                                                     ccTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTc                           >  7:1225/1‑149 (MQ=255)
                                                                                                                          gTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTcgc                         <  7:227519/147‑1 (MQ=255)
                                                                                                                          gTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCgcgc                       <  8:1225/149‑1 (MQ=255)
                                                                                                                             tAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTg                     <  5:390527/148‑1 (MQ=255)
                                                                                                                                gCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTAt                 >  1:17302/1‑149 (MQ=255)
                                                                                                                                  tATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTAttt               >  3:199305/1‑149 (MQ=255)
                                                                                                                                    tttGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACATCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTAttt               <  5:315839/147‑1 (MQ=255)
                                                                                                                                    tttGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACATCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTAttt               <  5:80202/147‑1 (MQ=255)
                                                                                                                                     ttagggCCGCTTCGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGc            >  2:275304/4‑149 (MQ=255)
                                                                                                                                      tagggCCGCTTCGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGc            >  1:208568/3‑148 (MQ=255)
                                                                                                                                       ggggCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCATGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTa          >  3:128266/1‑149 (MQ=255)
                                                                                                                                        gggCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTAc         <  7:95507/149‑1 (MQ=255)
                                                                                                                                           ccGCGAAGCGGCCCTAATTTCCCCCGGCTGATGTGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAACCACGCCACCAGGCAGGTATTGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTAcg        <  1:125982/149‑1 (MQ=255)
                                                                                                                                           ccGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTAcgcg      <  3:54613/149‑1 (MQ=255)
                                                                                                                                           ccGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTAcgcg      <  2:17302/149‑1 (MQ=255)
                                                                                                                                             gcgAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTACGCGCt    >  3:267100/1‑149 (MQ=255)
                                                                                                                                             gcgAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTACGCGCt    >  8:188999/1‑149 (MQ=255)
                                                                                                                                              cgAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTACGCGCt    >  6:127282/1‑148 (MQ=255)
                                                                                                                                               gAAGCGGCCCTAATTTCCCCCGTCTGATGCGGTCTTAGCGGGTAGCAATACAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTAcgcg      <  6:117349/145‑1 (MQ=255)
                                                                                                                                               gAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTACGCGCTCa  <  4:199305/149‑1 (MQ=255)
                                                                                                                                               gAAGCGGCCCTAATTTCACCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCATGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATACGCCGTCTTCGAGCTGGTATATGGCTAagcg      >  5:117349/1‑145 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TCGAGGCGTTGTACCCCGAGGTGGACTTTGCCTGGCTGGAGTTCGAGCGGGGCGGGCATGGGGTGTTCATGCTCACCGCCGAGCAGTGCCGCCAGCATCAGGAATTGTTCAAGGCCCTCGTTTAGGCTATTTGGGGCCGCGAAGCGGCCCTAATTTCCCCCGGCTGATGCGGTCTTAGCGGGTAGCAATCCAGATCAACAGCCCTGCCTGGAACACGGCAA‑‑ACGCCACCAGGCAGGTAATGGTGAAACGCAATCCGCCGTCTTCGCGCTGGTATTTGGCTACGCGCTCA  >  NC_002947/2052644‑2052932

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: