Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,430,392 A→G 6.5% intergenic (‑74/+87) PP_4769 ← / ← PP_4770 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,430,3920AG 6.5% 176.9 / 2.6 61intergenic (‑74/+87)PP_4769/PP_4770hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base G (2/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

GCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCAC  >  NC_002947/5430250‑5430534
                                                                                                                                                |                                                                                                                                               
gCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACATCATTTTTGGCATGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACATTCGCGGACACGCCCCATCCCACAGGTATc                                                                                                                                          >  7:6523/1‑149 (MQ=255)
    ccGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATcgcgc                                                                                                                                      <  5:228100/149‑1 (MQ=255)
            cTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACAccc                                                                                                                              <  3:261127/149‑1 (MQ=255)
              ggCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATGTTTTTCACCTGGGAACCGAATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGAGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGt                                                                                                                            >  5:157118/1‑149 (MQ=255)
               gcgcACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTg                                                                                                                           <  5:117947/149‑1 (MQ=255)
                  cACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGCGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTgg                                                                                                                          <  1:112925/149‑1 (MQ=255)
                    cGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGc                                                                                                                      <  5:337135/149‑1 (MQ=37)
                     ggTGTAGTAGAGGTTCTGCACATCGGCGTACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCg                                                                                                                     >  1:211774/1‑149 (MQ=255)
                      gtgtAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCgg                                                                                                                    <  2:210042/149‑1 (MQ=255)
                         tagtagACGGTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGc                                                                                                                  <  7:219095/148‑1 (MQ=255)
                            tagAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGc                                                                                                              <  5:143817/149‑1 (MQ=255)
                            tagAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGCACGCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGCGGGCGTgtc                                                                                                             >  5:207109/1‑147 (MQ=25)
                              gagGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGc                                                                                                                  >  6:358503/1‑143 (MQ=255)
                              gagGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGc                                                                                                                  <  5:358503/143‑1 (MQ=255)
                              gagGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGCACGCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGCGGGCGTGTCcg                                                                                                           <  6:207109/149‑1 (MQ=25)
                               agGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCcg                                                                                                           >  6:358358/1‑149 (MQ=255)
                                 gTTCTGCACTTCGCCGAAAATCGCTATATTTTTCACCTGGGAACCTCATTTTTGGCCTGGTGTGCTG‑‑CGGTGATGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTa                                                                                                                                             <  1:273107/114‑1 (MQ=37)
                                  ttCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCcgcg                                                                                                         >  8:235999/1‑148 (MQ=255)
                                     agCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGTGAAg                                                                                                      <  5:166989/147‑1 (MQ=255)
                                      gCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                    >  2:89595/1‑149 (MQ=255)
                                        acaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTa                                                                                                                                   <  7:187512/116‑1 (MQ=255)
                                        acaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTa                                                                                                                                   >  8:187512/1‑116 (MQ=255)
                                          aTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGgcg                                                                                                 >  4:40103/1‑148 (MQ=255)
                                          aTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGgcg                                                                                                 >  4:312777/1‑148 (MQ=255)
                                             ggCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCaaa                                                                                             <  2:211774/149‑1 (MQ=255)
                                             ggCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCaaa                                                                                             >  4:133235/1‑149 (MQ=255)
                                               cGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCCCCCCCCCACGTTTAGCGCGGGAACCCCGGGGGGCGGGCGGGGCCCGGCAGGGGCGCcaacg                                                                                           >  2:354319/1‑147 (MQ=255)
                                                gAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAgcg                                                                                          >  8:114088/1‑149 (MQ=255)
                                                  aCAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGcc                                                                                        <  6:165996/149‑1 (MQ=255)
                                                  aCAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGcc                                                                                        <  6:107874/149‑1 (MQ=255)
                                                    aGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAgcgc                                                                                         <  8:342642/146‑1 (MQ=255)
                                                    aGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAgcgc                                                                                         >  7:342642/1‑146 (MQ=255)
                                                       gcgATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCg                                                                                    >  3:17256/1‑148 (MQ=255)
                                                               tttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGaaa                                                                                <  3:175223/144‑1 (MQ=255)
                                                               tttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGAGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAAt                                                                               <  3:368777/145‑1 (MQ=255)
                                                                ttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGTCCCAAAGCGCCCCCGGAAATGCca                                                                           <  7:235999/148‑1 (MQ=255)
                                                                   aCCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCcagcag                                                                       >  5:66081/1‑149 (MQ=255)
                                                                             cTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAc                                                             >  7:95660/1‑149 (MQ=255)
                                                                                 tttttGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCAACAGGTATCGCGCGTACACCCGTGGTAGGGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGt                                                         <  1:89595/149‑1 (MQ=255)
                                                                                       gCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTGCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGAc                                                    <  7:128191/148‑1 (MQ=255)
                                                                                                  cG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                        <  8:6523/149‑1 (MQ=255)
                                                                                                      gggggCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATcgtcgt                                      <  3:333180/149‑1 (MQ=255)
                                                                                                        gggCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTcc                                    <  5:305184/149‑1 (MQ=255)
                                                                                                        gggCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTcc                                    <  4:17256/149‑1 (MQ=255)
                                                                                                          gctgctGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTcc                                    >  8:331102/1‑147 (MQ=255)
                                                                                                          gctgctGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCga                                  <  3:133235/149‑1 (MQ=255)
                                                                                                              ctgcGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGagaa                               <  6:66081/148‑1 (MQ=255)
                                                                                                                 cgcAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGGCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAAc                           >  5:293308/1‑149 (MQ=255)
                                                                                                                      aCCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATcgtcgt                                      <  7:254545/133‑1 (MQ=255)
                                                                                                                      aCCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATcgtcgt                                      >  8:254545/1‑133 (MQ=255)
                                                                                                                           cgcgGACACGCCGGCTCCCACAGGTATCGCGCGTACCCCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTc                   <  1:371641/147‑1 (MQ=255)
                                                                                                                            gcgGACACGCCCCCTCCCACAGGTATCGCGCGTAAACCCGTGGGAGCGGGCCTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTAAGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATc                >  6:19161/1‑149 (MQ=255)
                                                                                                                             cgGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCg               >  1:30132/1‑149 (MQ=255)
                                                                                                                                 cacGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGc           <  3:40103/149‑1 (MQ=255)
                                                                                                                                  acGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGcc          <  4:27791/149‑1 (MQ=255)
                                                                                                                                  acGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGc           >  6:199382/1‑148 (MQ=255)
                                                                                                                                  acGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGc           >  2:45839/1‑148 (MQ=255)
                                                                                                                                    gCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGCGGGCGTGCGCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCa         >  2:34534/1‑148 (MQ=37)
                                                                                                                                    gCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGCACCCACCGTCATCGCTGCCa         >  1:53600/1‑148 (MQ=255)
                                                                                                                                    gCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCaa        <  6:233812/149‑1 (MQ=255)
                                                                                                                                          tCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGACAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCAACGCTGCCCAAaccac  >  8:165304/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCG‑‑GGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCAC  >  NC_002947/5430250‑5430534

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: