Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,288,998 Δ1 bp 5.7% coding (169/822 nt) PP_3760 ← methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,288,9980C. 5.7% 218.8 / 7.2 70coding (169/822 nt)PP_3760methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base . (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 = 9.14e-01

ACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTC  >  NC_002947/4288856‑4289145
                                                                                                                                              |                                                                                                                                                     
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                                                                                                 >  2:1166182/1‑150 (MQ=255)
                                              gcACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGga                                                                                                 <  1:1166182/150‑1 (MQ=255)
                                                  ggCAAGGAAGTGGTCAGGGTCACGGAACATCACGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAAt                                                                                             >  2:872896/1‑150 (MQ=255)
                                                      acgAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCgg                                                                                         <  2:1101061/148‑1 (MQ=255)
                                                                   ggTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAACAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTg                                                                                                   <  2:1073527/127‑1 (MQ=255)
                                                                   ggTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAACAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTg                                                                                                   >  1:1073527/1‑127 (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                                                                       <  1:405616/124‑1 (MQ=255)
                                                                                                 tCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGcc                                                                       >  2:405616/1‑124 (MQ=255)
                                                                                                    ggTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACt                                          >  1:1055864/1‑150 (MQ=255)
                                                                                                          gcagcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGt                                    >  1:400396/1‑150 (MQ=255)
                                                                                                             gcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGc                                                            >  1:279803/1‑123 (MQ=255)
                                                                                                             gcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGc                                                            <  2:279803/123‑1 (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:24231/150‑1 (MQ=255)
                                                                                                                       gAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGAt                       >  1:1197045/1‑150 (MQ=255)
                                                                                                                             gCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGATAATCCGGCGCTTGATCGaa                                                                            <  2:728987/91‑1 (MQ=255)
                                                                                                                             gCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGaa                                                                            >  1:728987/1‑91 (MQ=255)
                                                                                                                               ccGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTc               <  2:400396/150‑1 (MQ=255)
                                                                                                                               ccGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTc               >  1:467479/1‑150 (MQ=255)
                                                                                                                                cgtgTCGTGCAGGATGCTCTTCTGTAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTCAGGTAGATCGCCTCGATCAGCAGTCg              <  2:564297/147‑1 (MQ=255)
                                                                                                                                atcccgttGCAGGAGCGCGTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCg              <  1:1285783/143‑1 (MQ=255)
                                                                                                                                     gttGCAGGAGCGCGTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGAtct         >  1:695470/3‑150 (MQ=255)
                                                                                                                                     gttGCAGGAGCGCGTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGCGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAt                                <  2:55946/125‑1 (MQ=255)
                                                                                                                                     gttGCAGGAGCGCGTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGCGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAt                                >  1:55946/3‑127 (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)
                                                                                                                                              |                                                                                                                                                     
ACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTGATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTC  >  NC_002947/4288856‑4289145

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: