Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,430,452 Δ1 bp 6.3% intergenic (‑86/+75) 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,4520T. 6.3% 185.6 / 7.0 63intergenic (‑86/+75)PP_4769/PP_4770hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base . (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

GGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCCT  >  NC_002947/5430307‑5430600
                                                                                                                                                   |                                                                                                                                                    
ggCCTGGTGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTatac                                                                                                                                                  <  1:1285252/150‑1 (MQ=255)
ggCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGCCCTTCACGGACAGGCCGGCTCCCACAGGTATCGCGCGTacac                                                                                                                                                  >  1:925874/1‑150 (MQ=255)
    gggCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGt                                                                                                                                              <  2:200796/149‑1 (MQ=255)
          aCGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGGCACGCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGc                                                                                                                                        <  1:968657/150‑1 (MQ=25)
          aCGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGc                                                                                                                                        >  1:872956/1‑150 (MQ=37)
             gtgtAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGCggg                                                                                                                                     <  2:872956/150‑1 (MQ=37)
             gtgtAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGc                                                                                                                                    <  2:1088042/150‑1 (MQ=255)
                tagtagAGGTTATGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGt                                                                                                                                  <  2:1222963/150‑1 (MQ=255)
                          tCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAg                                                                                                                        <  2:533949/150‑1 (MQ=255)
                           cTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTTGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAgg                                                                                                                       <  2:22881/150‑1 (MQ=255)
                              cacaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGgc                                                                                                                    <  1:278750/150‑1 (MQ=255)
                                 aTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCa                                                                                                                 >  1:510328/1‑150 (MQ=255)
                                 aTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCa                                                                                                                 >  1:1067797/1‑150 (MQ=255)
                                     gCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTacac                                                                                                                                                  <  1:930903/113‑1 (MQ=255)
                                     gCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTacac                                                                                                                                                  >  2:930903/1‑113 (MQ=255)
                                       gAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                                      >  1:339280/1‑139 (MQ=255)
                                       gAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                                      <  2:339280/139‑1 (MQ=255)
                                            gcgcGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCg                                                                                                      >  2:1279596/1‑150 (MQ=255)
                                            gcgcGATCTTTTTCACCTGGGAACCTCATTTTTGGCATGGAGTGAAGGAGGGTTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAATGGGCGCAAAGCGCCCCCg                                                                                                      >  1:863671/1‑150 (MQ=39)
                                              gcgATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGa                                                                                                    >  1:1174772/1‑150 (MQ=255)
                                                gATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACCGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGaaa                                                                                                  <  2:863671/150‑1 (MQ=255)
                                                 aTCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACAcc                                                                                                                                                 >  2:783190/1‑102 (MQ=255)
                                                 aTCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACAcc                                                                                                                                                 <  1:783190/102‑1 (MQ=255)
                                                  tCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTg                                                                                                                                             >  2:118081/1‑105 (MQ=255)
                                                  tCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTg                                                                                                                                             <  1:118081/105‑1 (MQ=255)
                                                    tttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGcc                                                                                              >  1:575580/1‑150 (MQ=255)
                                                    tttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGcc                                                                                              >  2:807508/1‑150 (MQ=255)
                                                    tttttCACCTGGGAACCCCATTTTTTTTTTTTTTTTCCGTGGTGCTGCTGCGCAGACCTTCTTTTTCATGCCCGCTCCCACAGGTATCGCGCTTACACCCGTGGTATCGTGCGTGCCCGCGAAGGGGCTCAAAGCGCCCCCGGAAATGcc                                                                                              >  2:487742/1‑150 (MQ=255)
                                                     ttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCca                                                                                             <  2:723926/150‑1 (MQ=255)
                                                      tttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCcag                                                                                            <  2:962226/150‑1 (MQ=255)
                                                           ccTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTa                                                                                       <  2:1067797/150‑1 (MQ=255)
                                                                gAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCa                                                                                  <  2:510328/150‑1 (MQ=255)
                                                                gAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCa                                                                                  <  1:501002/150‑1 (MQ=255)
                                                                  aCCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGa                                                                                <  1:363304/150‑1 (MQ=255)
                                                                   ccTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAc                                                                               <  2:575580/150‑1 (MQ=255)
                                                                    cTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGGGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAcc                                                                              >  1:785656/1‑150 (MQ=255)
                                                                    cTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAcc                                                                              <  2:94759/150‑1 (MQ=255)
                                                                     tCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCt                                                                             >  2:74390/1‑150 (MQ=255)
                                                                       aTTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCAGAAATGCCAGCAGTATGCCAGACCTGt                                                                           >  2:263291/1‑150 (MQ=255)
                                                                        tttttGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAgg                                                                                                                       >  2:438908/1‑105 (MQ=255)
                                                                        tttttGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAgg                                                                                                                       <  1:438908/105‑1 (MQ=255)
                                                                            tGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGAc                                                                      <  2:925874/150‑1 (MQ=255)
                                                                             ggCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACg                                                                     <  1:807508/150‑1 (MQ=255)
                                                                              gCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGt                                                                    >  1:1170953/1‑150 (MQ=255)
                                                                                cTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTaa                                                                  >  1:1002555/1‑150 (MQ=255)
                                                                                cTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTaa                                                                  >  2:372628/1‑150 (MQ=255)
                                                                                cTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTaa                                                                  >  2:1085997/1‑150 (MQ=255)
                                                                                 tGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAgc                                                                                                             <  2:812129/106‑1 (MQ=255)
                                                                                 tGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAgc                                                                                                             >  1:812129/1‑106 (MQ=255)
                                                                                        ccGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                                          <  1:1279596/150‑1 (MQ=255)
                                                                                              ggCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCga                                                    >  1:693691/1‑150 (MQ=255)
                                                                                                    tgcGCAGACCTTCGCGGACACGCCCGCTCCCACAGGAGTCAGCGCTTCGCCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAgaagaa                                              >  2:887702/1‑150 (MQ=38)
                                                                                                     gcgcAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCGGACGTAAGAATCGTCGTCCGAGAAGAAc                                             <  1:74390/150‑1 (MQ=255)
                                                                                                     cccacctccTTCGCGGACACGCCCCCTCCCCCAGCTCTCGCGCGTCCACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAAc                                             <  2:1174652/143‑1 (MQ=255)
                                                                                                           aCCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCg                                       >  2:748126/1‑150 (MQ=255)
                                                                                                             cTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAAccagatcgg                                     >  1:1174652/1‑143 (MQ=255)
                                                                                                                           ccGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGCGGGCGTGTCCGCGAAGGGGAGCAGGGAGCCCCCGGAGTATTCGGGAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAac                       >  1:156801/1‑150 (MQ=12)
                                                                                                                           ccGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAac                       >  2:1029671/1‑150 (MQ=255)
                                                                                                                                 ccACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTgtc                                                                                                                               >  2:553923/1‑38 (MQ=255)
                                                                                                                                 ccACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTgtc                                                                                                                               <  1:553923/40‑3 (MQ=255)
                                                                                                                                              gcgcgTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGgc                                                                                                                    <  1:219329/40‑1 (MQ=255)
                                                                                                                                              gcgcgTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGgc                                                                                                                    >  2:219329/1‑40 (MQ=255)
                                                                                                                                                  gTACACCAGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCCt  >  1:644977/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCCT  >  NC_002947/5430307‑5430600

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