Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,681,711 C→G 7.3% A63G (GCC→GGC)  recJ → ssDNA‑specific exonuclease RecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,681,7110CG 7.3% 201.2 / 5.5 64A63G (GCC→GGC) recJssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base C (33/26);  new base G (3/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTACG  >  NC_002947/1681568‑1681857
                                                                                                                                               |                                                                                                                                                  
gggTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGa                                                                                                                                              >  1:991163/1‑150 (MQ=255)
   tAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAcca                                                                                                                                           >  2:537820/1‑150 (MQ=255)
     aCCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagc                                                                                                                                         <  1:7236/150‑1 (MQ=255)
     aCCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGATCGAAGCCATCGAccagc                                                                                                                                         >  2:20421/1‑150 (MQ=255)
      ccTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcg                                                                                                                                        >  2:1180530/1‑150 (MQ=255)
           cACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgccagc                                                                                                                                   <  1:572544/150‑1 (MQ=255)
           cACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgccagc                                                                                                                                   <  1:537820/150‑1 (MQ=255)
               cTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCAt                                                                                                                               >  1:944825/1‑150 (MQ=255)
                ttGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATc                                                                                                                                                <  2:697774/132‑1 (MQ=255)
                ttGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATc                                                                                                                                                >  1:697774/1‑132 (MQ=255)
                      aCCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtc                                                                                                                        >  1:737598/1‑150 (MQ=255)
                        ccGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtcgt                                                                                                                      >  2:1169472/1‑150 (MQ=255)
                        ccGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtcgt                                                                                                                      >  2:636018/1‑150 (MQ=255)
                          gCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtcgtcg                                                                                                                    >  2:948721/1‑150 (MQ=255)
                            cTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGc                                                                                                                  <  2:987050/150‑1 (MQ=255)
                             tgcaCGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCg                                                                                                                 >  2:692229/4‑150 (MQ=255)
                                  gccgcGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTc                                                                                                            >  1:565771/1‑150 (MQ=255)
                                  gccgcGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTAACAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTc                                                                                                            >  1:752406/1‑150 (MQ=255)
                                     gggcgtcgcgTGCAGAGGGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCTGGTGTGGACCTTCTGGTCGAAGCCATCGATCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGAt                                                                                                         <  2:328999/143‑1 (MQ=255)
                                     gcgcgcgGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGAt                                                                                                         <  1:692229/150‑1 (MQ=255)
                                     gcgcgcgGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGAt                                                                                                         <  1:948721/150‑1 (MQ=255)
                                       gcgcgGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGc                                                                                                       >  1:795821/1‑150 (MQ=255)
                                                 cAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGAt                                                                                                         >  2:1181460/1‑138 (MQ=255)
                                                 cAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGAt                                                                                                         <  1:1181460/138‑1 (MQ=255)
                                                  agagCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCa                                                                                                                                                  >  1:396879/1‑96 (MQ=255)
                                                  agagCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCa                                                                                                                                                  <  2:396879/96‑1 (MQ=255)
                                                  agagCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCa                                                                                            <  2:106321/150‑1 (MQ=255)
                                                     gCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGATGGCTTCGAccagcagg                                                                                                                                      <  1:819963/105‑4 (MQ=255)
                                                     gCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGATGGCTTCGAccagcagg                                                                                                                                      >  2:819963/1‑102 (MQ=255)
                                                                 tGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCg                                                                             <  1:1180530/150‑1 (MQ=255)
                                                                  ggACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCgg                                                                            <  2:991163/150‑1 (MQ=255)
                                                                   gACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCggg                                                                           <  1:636018/150‑1 (MQ=255)
                                                                        gagCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCt                                                                      >  1:1064210/1‑150 (MQ=255)
                                                                         agCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTg                                                                     >  2:408562/1‑150 (MQ=255)
                                                                               gCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGgcctg                                                               <  2:1064210/150‑1 (MQ=255)
                                                                                 aCGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGgcctgcg                                                             >  1:844774/1‑150 (MQ=255)
                                                                                   gCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGgcctgcgcc                                                           >  1:281012/1‑150 (MQ=255)
                                                                                     ctgctgCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGgcctgcgcctg                                                         >  2:906853/1‑150 (MQ=255)
                                                                                       gctgcCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCg                                                                                         <  1:1000091/116‑1 (MQ=255)
                                                                                       gctgcCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCg                                                                                         >  2:1000091/1‑116 (MQ=255)
                                                                                           cccTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGt                                                   >  2:1153753/1‑150 (MQ=255)
                                                                                              tACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCg                                                >  2:257182/1‑150 (MQ=255)
                                                                                                  agcagcTCAAGGGCATCGAGGCAGGTGTGGAACTCCTATTAGCCGGCATGCAACAGCGCCATCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCCCCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGccc                                            <  2:281012/150‑1 (MQ=255)
                                                                                                    cagcTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGATCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCAt                                          <  1:20421/150‑1 (MQ=255)
                                                                                                       cTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTc                                       >  1:27521/1‑150 (MQ=255)
                                                                                                       cTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTc                                       >  1:1006102/1‑150 (MQ=255)
                                                                                                        tCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGt                                                   <  1:943772/137‑1 (MQ=255)
                                                                                                        tCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGAGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGt                                                   >  2:943772/1‑137 (MQ=255)
                                                                                                          aaGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGAt                                    <  1:257182/150‑1 (MQ=255)
                                                                                                                aTCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTg                              <  2:737598/150‑1 (MQ=255)
                                                                                                                   gAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTg                                                                     >  2:20650/1‑108 (MQ=255)
                                                                                                                   gAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTg                                                                     <  1:20650/108‑1 (MQ=255)
                                                                                                                    aGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGc                          <  2:844774/150‑1 (MQ=255)
                                                                                                                             tgGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGt                                                   >  1:855224/1‑116 (MQ=255)
                                                                                                                             tgGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGt                                                   <  2:855224/116‑1 (MQ=255)
                                                                                                                                    gctgGTCGATGGCTTCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTa          >  2:728067/1‑150 (MQ=255)
                                                                                                                                    gctgGTCGATGGCTTCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTa          >  2:226426/1‑150 (MQ=255)
                                                                                                                                    gctgGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTa          >  1:153955/1‑150 (MQ=255)
                                                                                                                                    gctgGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTa          <  1:906853/150‑1 (MQ=255)
                                                                                                                                     ctgGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTAt         >  2:1030438/1‑150 (MQ=255)
                                                                                                                                      tgGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATg        <  2:944825/150‑1 (MQ=255)
                                                                                                                                        gTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGc      >  1:871905/1‑150 (MQ=255)
                                                                                                                                           gAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTAc   >  1:70202/1‑150 (MQ=255)
                                                                                                                                            aaGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTACg  <  1:1169472/150‑1 (MQ=255)
                                                                                                                                               ccATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTa    <  2:1051222/145‑1 (MQ=255)
                                                                                                                                               ccATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTa    >  1:1051222/1‑145 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTACG  >  NC_002947/1681568‑1681857

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