Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,289,002:1 +G 5.7% coding (165/822 nt) PP_3760 ← methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,289,0021.G 5.7% 281.3 / 9.3 70coding (165/822 nt)PP_3760methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (33/33);  new base G (2/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

ACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTT  >  NC_002947/4288856‑4289151
                                                                                                                                                   |                                                                                                                                                      
aCTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTcc                                                                                                                                                     >  1:70799/1‑150 (MQ=255)
 cTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCt                                                                                                                                                    <  1:258054/150‑1 (MQ=255)
   ttGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGc                                                                                                                                                  <  2:481155/150‑1 (MQ=255)
       tCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACg                                                                                                                                              >  1:997931/1‑150 (MQ=255)
       tCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACg                                                                                                                                              >  1:1376562/1‑150 (MQ=255)
           cGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGa                                                                                                                                          <  1:736905/150‑1 (MQ=255)
             gCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGaca                                                                                                                                        >  1:1280596/1‑150 (MQ=255)
               cAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACAcc                                                                                                                                      >  2:228400/1‑150 (MQ=255)
                 ggTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGt                                                                                                                                    <  2:1411142/150‑1 (MQ=255)
                    gcgcAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAg                                                                                                                                 <  2:921122/150‑1 (MQ=255)
                          cagcGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATc                                                                                                                           <  2:864512/150‑1 (MQ=255)
                          cagcGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATc                                                                                                                           <  1:67781/150‑1 (MQ=255)
                            gcGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGa                                                                                                                         <  1:482228/150‑1 (MQ=255)
                               gCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACt                                                                                                                      <  2:70799/150‑1 (MQ=255)
                                 accacTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGa                                                                                                                    >  2:535909/1‑150 (MQ=255)
                                   cacTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGAACCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACg                                                                                                                  >  1:13545/1‑150 (MQ=255)
                                    acTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGc                                                                                                                 >  2:24231/1‑150 (MQ=255)
                                              gcACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGga                                                                                                       <  1:1166182/150‑1 (MQ=255)
                                              gcACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGga                                                                                                       >  2:1166182/1‑150 (MQ=255)
                                                  ggCAAGGAAGTGGTCAGGGTCACGGAACATCACGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAAt                                                                                                   >  2:872896/1‑150 (MQ=255)
                                                      acgAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCgg                                                                                               <  2:1101061/148‑1 (MQ=255)
                                                                   ggTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAACAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTg                                                                                                         >  1:1073527/1‑127 (MQ=255)
                                                                   ggTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAACAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTg                                                                                                         <  2:1073527/127‑1 (MQ=255)
                                                                      cACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGTAACGCAGACACCGTCAGAGAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAgcgc                                                                              >  1:564297/1‑150 (MQ=255)
                                                                         ggAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCgg                                                                           >  2:1176865/1‑150 (MQ=255)
                                                                         ggAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCgg                                                                           >  2:255001/1‑150 (MQ=255)
                                                                          gAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATCAAGCGCCGGa                                                                          <  2:1399621/150‑1 (MQ=255)
                                                                               tctcGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCg                                                                                                                          >  1:1144554/1‑98 (MQ=255)
                                                                               tctcGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCg                                                                                                                          <  2:1144554/98‑1 (MQ=255)
                                                                               tctcGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGaataat                                                                     >  2:1315485/1‑150 (MQ=255)
                                                                                 cgatctaCCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTg                                                                   <  1:156254/144‑1 (MQ=255)
                                                                                        ccGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGAGTGGAGAATCCGGTGCTGGATCGAAGGCCCGGGATAATTGcagatcg                                                            >  2:156254/1‑144 (MQ=255)
                                                                                         cGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATc                                                           >  1:1150302/1‑150 (MQ=255)
                                                                                            gATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAAt                                                                                                                            >  2:15158/1‑83 (MQ=255)
                                                                                            gATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAAt                                                                                                                            <  1:15158/83‑1 (MQ=255)
                                                                                            gATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTa                                                        <  1:872896/150‑1 (MQ=255)
                                                                                                 tCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGcc                                                                             >  2:405616/1‑124 (MQ=255)
                                                                                                 tCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGcc                                                                             <  1:405616/124‑1 (MQ=255)
                                                                                                    ggTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACt                                                >  1:1055864/1‑150 (MQ=255)
                                                                                                          gcagcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGt                                          >  1:400396/1‑150 (MQ=255)
                                                                                                             gcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGc                                                                  <  2:279803/123‑1 (MQ=255)
                                                                                                             gcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGc                                                                  >  1:279803/1‑123 (MQ=255)
                                                                                                             gcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTaa                                                       >  1:850903/1‑133 (MQ=255)
                                                                                                             gcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAg                                                       <  2:850903/134‑1 (MQ=255)
                                                                                                              caACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGATGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAg                                            >  1:1259981/1‑144 (MQ=255)
                                                                                                              caACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGATGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAg                                            <  2:1259981/144‑1 (MQ=255)
                                                                                                                  tgcatgcaCATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATc                                     <  1:929144/147‑1 (MQ=255)
                                                                                                                  tGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATc                                     >  2:929144/1‑147 (MQ=255)
                                                                                                                  tGCATGAAAATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGcc                                  <  1:292414/150‑1 (MQ=255)
                                                                                                                   gCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAg                                                                                                               >  1:279747/1‑73 (MQ=255)
                                                                                                                   gCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAg                                                                                                               <  2:279747/73‑1 (MQ=255)
                                                                                                                    cATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTc                                <  2:1280596/150‑1 (MQ=255)
                                                                                                                    cATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTc                                <  2:1376562/150‑1 (MQ=255)
                                                                                                                       gAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGAt                             >  1:1197045/1‑150 (MQ=255)
                                                                                                                       gAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGAt                             <  1:24231/150‑1 (MQ=255)
                                                                                                                             gCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGATAATCCGGCGCTTGATCGaa                                                                                  <  2:728987/91‑1 (MQ=255)
                                                                                                                             gCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGaa                                                                                  >  1:728987/1‑91 (MQ=255)
                                                                                                                               ccGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTc                     >  1:467479/1‑150 (MQ=255)
                                                                                                                               ccGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTc                     <  2:400396/150‑1 (MQ=255)
                                                                                                                                cgtgTCGTGCAGGATGCTCTTCTGTAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTCAGGTAGATCGCCTCGATCAGCAGTCg                    <  2:564297/147‑1 (MQ=255)
                                                                                                                                atcccgttGCAGGAGCGCGTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCg                    <  1:1285783/143‑1 (MQ=255)
                                                                                                                                     gttGCAGGAGCGCGTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGAtct               >  1:695470/3‑150 (MQ=255)
                                                                                                                                     gttGCAGGAGCGCGTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGCGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAt                                      >  1:55946/3‑127 (MQ=255)
                                                                                                                                     gttGCAGGAGCGCGTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGCGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAt                                      <  2:55946/125‑1 (MQ=255)
                                                                                                                                       tGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCg             <  1:535909/150‑1 (MQ=255)
                                                                                                                                       tGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCg             >  1:1408991/1‑150 (MQ=255)
                                                                                                                                          aGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATg          >  1:385055/1‑150 (MQ=255)
                                                                                                                                            gACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTc        <  2:385055/150‑1 (MQ=255)
                                                                                                                                             aCGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGAt                                                                                      >  1:1310684/1‑71 (MQ=255)
                                                                                                                                             aCGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGAt                                                                                      <  2:1310684/71‑1 (MQ=255)
                                                                                                                                               gcgcTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGgt     <  1:228400/150‑1 (MQ=255)
                                                                                                                                                  cTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGtt  <  1:1176865/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTT  >  NC_002947/4288856‑4289151

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