Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,291,548 C→T 6.1% D41D (GAC→GAT)  PP_4665 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,291,5480CT 6.1% 192.5 / 2.4 65D41D (GAC→GAT) PP_4665hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/33);  new base T (2/2);  total (30/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.49e-01

GCCAGTCACCAGGAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTGAAAGCAGGTGATGTGT  >  NC_002947/5291404‑5291696
                                                                                                                                                |                                                                                                                                                    
gCCAGTCACCAGGAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGgca                                                                                                                                                  <  3:321317/149‑1 (MQ=255)
  cAGTCACCAGGAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGgcag                                                                                                                                                 >  3:339834/1‑148 (MQ=255)
   aGTCACCAGGAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGgcagc                                                                                                                                                >  8:280099/1‑148 (MQ=255)
   aGTCACCAGGAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGgcagc                                                                                                                                                >  6:164948/1‑148 (MQ=255)
    gTCACCAGGAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCaa                                                                                                                                              <  2:121133/149‑1 (MQ=255)
    gTCACCAGGAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCaa                                                                                                                                              <  7:207508/149‑1 (MQ=255)
         cAGGAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCa                                                                                                                                         >  5:189852/1‑149 (MQ=255)
          aGGAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCAt                                                                                                                                        <  8:272826/149‑1 (MQ=255)
           ggAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGTCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATg                                                                                                                                       >  6:313416/1‑149 (MQ=255)
            gAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATgg                                                                                                                                      <  5:272914/149‑1 (MQ=255)
            gAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATgg                                                                                                                                      <  7:280099/149‑1 (MQ=255)
                  tgCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGATGGCAGCAAGATCATGGGCACCt                                                                                                                                >  6:195481/1‑149 (MQ=255)
                  tgCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCt                                                                                                                                >  6:169853/1‑149 (MQ=255)
                    ccATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTgg                                                                                                                              <  1:329263/149‑1 (MQ=255)
                          ccGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTg                                                                                                                         >  5:376945/1‑148 (MQ=255)
                           cGATTTGATCACCGTCCTTGGCGAAACCTGCGCGACGCCGGTCGTTGAATCCACAAAGTGTAATCGCATCGGCGGCGTTTCGCACAGAGTCTACCTCAATGCCTATCTGTAGGCTGACTGTAGCAAGATCATGGGCACCTGTATCTGCa                                                                                                                       <  4:66136/149‑1 (MQ=255)
                           cGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGc                                                                                                                        <  6:40351/148‑1 (MQ=255)
                            gATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCAc                                                                                                                      >  7:10076/1‑149 (MQ=255)
                                  atcaCCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCGCCTGGATCTGCACCCCgg                                                                                                                 <  6:189852/148‑1 (MQ=255)
                                     aCCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCaa                                                                                                             <  5:164948/149‑1 (MQ=255)
                                              cgcgAAACCTGCCCGACGTCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAgg                                                                                                     <  5:313416/148‑1 (MQ=255)
                                              cgcgAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGATGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAg                                                                                                      <  5:195481/147‑1 (MQ=255)
                                               gcgAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGt                                                                                                    >  8:224891/1‑148 (MQ=255)
                                               gcgAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTc                                                                                                   <  5:169853/149‑1 (MQ=255)
                                                              cGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGAGCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATgg                                                                                    <  5:330318/149‑1 (MQ=255)
                                                                ccGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGa                                                                                  <  7:224891/149‑1 (MQ=255)
                                                                      gTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCt                                                                             <  5:5733/148‑1 (MQ=255)
                                                                      gTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCt                                                                             <  6:376945/148‑1 (MQ=255)
                                                                      gTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCt                                                                             >  6:5733/1‑148 (MQ=255)
                                                                        ggccGCCACCTAGTCGAAACGCATCCGCGGCGAACCGCACACCGTCACCCACAATGCCAATCTGTCGTCTGACGGCAGCAAGATCAACTGCACCAGGATCTCCACCCCTGGCTAAATCGAGGTCAACTACGATAAATGGGAGTTCTGc                                                                           >  4:116367/4‑148 (MQ=255)
                                                                         gACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCa                                                                         >  6:27558/1‑149 (MQ=255)
                                                                              caccaAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTcc                                                                    >  4:348243/1‑149 (MQ=255)
                                                                              caccaAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTcc                                                                    >  3:256945/1‑149 (MQ=255)
                                                                                   aGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGTGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGAt                                                               <  6:255142/149‑1 (MQ=255)
                                                                                   aGTCGAAGAGGATCGGAGGCGATCCGCACACGGTCAACCTCAATGCCTATCTGTCTGATGACGGCAGCAAGATGATTGGAACCTGGTTCTGCACCCCGGTAAATTTCGATGTAAACTACTATAAATTGGAGTTCTGCCACTACCTCGAt                                                               <  4:339834/149‑1 (MQ=255)
                                                                                    gTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTc                                                                  <  4:260029/145‑1 (MQ=255)
                                                                                    gTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTc                                                                  >  3:260029/1‑145 (MQ=255)
                                                                                       gAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGc                                                            >  8:334505/1‑148 (MQ=255)
                                                                                         aGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTAc                                                         >  8:271826/1‑149 (MQ=255)
                                                                                             cATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACt                                                        >  4:100061/1‑146 (MQ=255)
                                                                                             cATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACt                                                        <  3:100061/146‑1 (MQ=255)
                                                                                              aTCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGca                                                     >  4:285531/1‑148 (MQ=255)
                                                                                                         tCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCg                                          >  6:163486/1‑148 (MQ=255)
                                                                                                             cacaCCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAgg                                     <  4:256945/149‑1 (MQ=255)
                                                                                                                  cGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCg                                   >  1:198973/1‑146 (MQ=255)
                                                                                                                  cGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGAACTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTACCTCGATGGCTACTGCATAATCACCCCGGAAGGCg                                   >  2:257627/1‑146 (MQ=255)
                                                                                                                      aaCCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATgg                                                                                    <  4:202030/93‑1 (MQ=255)
                                                                                                                      aaCCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATgg                                                                                    >  3:202030/1‑93 (MQ=255)
                                                                                                                      aaCCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGTGAAGAGc                             <  5:27558/148‑1 (MQ=255)
                                                                                                                        cgtgactGCCTCTATGTCGTCTTACGGCAGCAAGATCATGCGCACCTGGATCTGCACCCGGGCGAAGTTCTAGGGCAACTACGATAAATGGTAGTTCTGCCACTTGCTCGATGGCTTCTGCATCATCACCCCGGAAGGCGAAGAGCCGa                          <  3:116367/143‑1 (MQ=255)
                                                                                                                                 cTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTg                  >  1:214729/1‑148 (MQ=255)
                                                                                                                                 cTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTg                  >  6:372702/1‑148 (MQ=255)
                                                                                                                                  tATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTGa                 >  2:38200/1‑148 (MQ=255)
                                                                                                                                   aTCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGaa                         <  7:119124/139‑1 (MQ=255)
                                                                                                                                   aTCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGaa                         >  8:119124/1‑139 (MQ=255)
                                                                                                                                    tCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTGaaa               <  8:10076/148‑1 (MQ=255)
                                                                                                                                        tCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTGAAAGCAgg          <  2:198973/149‑1 (MQ=255)
                                                                                                                                         aggttGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTGAAAGCAGGt         <  2:362752/145‑1 (MQ=255)
                                                                                                                                           gCTGATGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAaga                               <  2:10362/125‑1 (MQ=255)
                                                                                                                                           gCTGATGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAaga                               >  1:10362/1‑125 (MQ=255)
                                                                                                                                           gCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTGAAAGCAGGTGa       <  7:334505/149‑1 (MQ=255)
                                                                                                                                           gCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTGAAAGCAGGTGa       <  2:83706/149‑1 (MQ=255)
                                                                                                                                           gCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTGAAAGCAGGTGa       >  6:18799/1‑149 (MQ=255)
                                                                                                                                            cTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGGCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTGAAAGCAGGTGAt      <  7:156298/149‑1 (MQ=255)
                                                                                                                                             tGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCt                                                                   <  6:108348/87‑1 (MQ=255)
                                                                                                                                                cGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTGAAAGCAGGTGATgtgt  <  5:372702/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCCAGTCACCAGGAGTTGTGCCATGTCCGATTTCATCACCGTCCTGCGCGAAACCTGCCCGACGCCGGTCGTGGACGCCACCAAGTGGAAGCGCATCGGCGGCGATCCGCACACCGTCAACCTCAATGCCTATCTGTCGGCTGACGGCAGCAAGATCATGGGCACCTGGATCTGCACCCCGGGCAAGTTCGAGGTCAACTACGATAAATGGGAGTTCTGCCACTTCCTCGATGGCTACTGCATCATCACCCCGGAAGGCGAAGAGCCGAAGCACCTGAAAGCAGGTGATGTGT  >  NC_002947/5291404‑5291696

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: ‑