Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,168,191 G→C 5.5% intergenic (+91/+597) PP_1921 → / ← PP_1922 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,168,1910GC 5.5% 228.9 / 3.6 73intergenic (+91/+597)PP_1921/PP_1922hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/39);  new base C (2/2);  total (32/41)
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.17e-01

GGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCGTCATCT  >  NC_002947/2168043‑2168336
                                                                                                                                                    |                                                                                                                                                 
ggATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCCACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCag                                                                                                                                                   <  1:250327/149‑1 (MQ=255)
        agagGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCa                                                                                                                                           <  2:98034/149‑1 (MQ=255)
           gtgACAAAACGACTAGCATCTATTTTTGTGTGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCTCTGCTAGCAcc                                                                                                                                         <  8:157210/146‑1 (MQ=255)
           gggACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCg                                                                                                                                        >  5:250463/1‑149 (MQ=255)
                aaaaCGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCtttt                                                                                                                                   <  3:152671/149‑1 (MQ=255)
                aaaaCGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCtttt                                                                                                                                   <  4:134157/149‑1 (MQ=255)
                aaaaCGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCtttt                                                                                                                                   <  5:367101/149‑1 (MQ=255)
                 aaaCGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTg                                                                                                                                  <  3:266653/149‑1 (MQ=255)
                  aaCGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTTGCACCGCTTTTg                                                                                                                                  >  1:360000/1‑148 (MQ=255)
                    cGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCa                                                                                                                               <  2:116972/149‑1 (MQ=255)
                        tAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATtgt                                                                                                                           <  7:142840/149‑1 (MQ=255)
                            aTCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGtct                                                                                                                       <  1:55278/149‑1 (MQ=255)
                             tCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCTCTGCTAGCACCGcatc                                                                                                                                    >  6:71327/1‑132 (MQ=255)
                             tCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCTCTGCTAGCACCGcatc                                                                                                                                    <  5:71327/135‑4 (MQ=255)
                                 tttttGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCgg                                                                                                                    <  5:76152/147‑1 (MQ=255)
                                 tttttGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCg                                                                                                                     <  8:211690/146‑1 (MQ=255)
                                 tttttGGGGGCTCTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCggg                                                                                                                   >  4:188086/1‑148 (MQ=255)
                                     tGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCggg                                                                                                                   >  4:340757/1‑144 (MQ=255)
                                     tGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCggg                                                                                                                   <  3:340757/144‑1 (MQ=255)
                                         ggCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGc                                                                                                                                <  7:17500/127‑1 (MQ=255)
                                         ggCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGc                                                                                                                                >  8:17500/1‑127 (MQ=255)
                                         ggCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAgctgct                                                                                                          <  5:96277/149‑1 (MQ=255)
                                            tttATGCTGCGTGACTACGCGGGATAGAGCTACCGAAAAGCGCGAATTAAGTAATATGTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCCCCGCTTTTGGAATtgt                                                                                                                           >  7:236033/1‑129 (MQ=255)
                                            tttATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGt                                                                                                       >  4:38059/1‑149 (MQ=255)
                                              tATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAgct                                                                                                             <  7:238471/141‑1 (MQ=255)
                                              tATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAgct                                                                                                             >  8:238471/1‑141 (MQ=255)
                                                tgctgcGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTg                                                                                                      <  3:188086/146‑1 (MQ=255)
                                                 gctgcGTGACTACGCGCCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGAGGTGCTAGCAGAGCTAGCACCGCTTTTGCCATTGTGTCTCGGGCAAGCTGCAGGAGAcc                                                                                                   >  4:144031/1‑148 (MQ=255)
                                                     cGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCtata                                                                                               <  7:198074/148‑1 (MQ=255)
                                                      gTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAgc                                                                                             <  8:49291/149‑1 (MQ=255)
                                                        gACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAgcgc                                                                                           >  2:123062/1‑149 (MQ=255)
                                                           tACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGctgctg                                                                                                         <  6:362570/132‑1 (MQ=255)
                                                           tACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGctgctg                                                                                                         >  5:362570/1‑132 (MQ=255)
                                                             cgcgGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGaaa                                                                                      <  3:91047/149‑1 (MQ=255)
                                                             cgcgGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGaa                                                                                       <  3:223768/148‑1 (MQ=255)
                                                                  cATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTgg                                                                                 <  6:121002/149‑1 (MQ=255)
                                                                   aTAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCTTATGTGCTTTGAAAGTCAGATGAGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCTGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTggg                                                                                <  7:322553/149‑1 (MQ=255)
                                                                   aTAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTggg                                                                                >  1:335996/1‑149 (MQ=255)
                                                                     agagCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTg                                                                              >  2:255436/1‑149 (MQ=255)
                                                                      gagCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGc                                                                             >  3:7061/1‑149 (MQ=255)
                                                                              gAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGt                                                                     <  7:84821/149‑1 (MQ=255)
                                                                               aaaaGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGTACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTc                                                                    >  8:299230/1‑149 (MQ=255)
                                                                                         aTTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAAt                                                          <  3:88358/149‑1 (MQ=255)
                                                                                         aTTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAAt                                                          <  1:123062/149‑1 (MQ=255)
                                                                                         aTTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAAt                                                          <  6:250463/149‑1 (MQ=255)
                                                                                         aTTAAGTAATATATTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAAt                                                          <  3:83105/149‑1 (MQ=255)
                                                                                            aaGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGt                                                        >  5:46457/1‑148 (MQ=255)
                                                                                            aaGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGt                                                        >  5:239608/1‑148 (MQ=255)
                                                                                             aGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGt                                                                                                       <  8:164038/100‑1 (MQ=255)
                                                                                             aGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGt                                                                                                       >  7:164038/1‑100 (MQ=255)
                                                                                                 atatCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTacaa                                                    >  2:279460/1‑147 (MQ=255)
                                                                                                   atCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTg                                                                         <  1:374491/124‑1 (MQ=255)
                                                                                                   atCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTg                                                                         >  2:374491/1‑124 (MQ=255)
                                                                                                    tCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGtt                                               >  3:248669/1‑149 (MQ=255)
                                                                                                    tCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGt                                                <  3:38059/148‑1 (MQ=255)
                                                                                                     cTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTg                                              >  3:303464/1‑149 (MQ=255)
                                                                                                         tttGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGGTCTCGAATGTACACCAGTTGGCAc                                          >  4:74697/1‑149 (MQ=255)
                                                                                                            ggggCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCg                                        <  6:46457/148‑1 (MQ=255)
                                                                                                             gggCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTc                                                                    <  2:316109/119‑1 (MQ=255)
                                                                                                                        gCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCtctc                                                              <  3:330367/114‑1 (MQ=255)
                                                                                                                        gCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCtctc                                                              >  4:330367/1‑114 (MQ=255)
                                                                                                                        gCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGc                           <  6:239608/149‑1 (MQ=255)
                                                                                                                         cTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCt                          <  5:213206/149‑1 (MQ=255)
                                                                                                                          tttGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCt                          >  8:221419/1‑148 (MQ=255)
                                                                                                                             gAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACt                      >  6:55785/1‑149 (MQ=255)
                                                                                                                              aaaGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACtt                     >  4:279668/1‑149 (MQ=255)
                                                                                                                                aGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCCCCGCTGCTCACAAGGCTGACtt                     >  7:293869/1‑147 (MQ=255)
                                                                                                                                  tCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACTTCAGt                 >  4:154185/1‑149 (MQ=255)
                                                                                                                                          ggTGCTAGCTCTGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACTTCAGTCCTTTTcc         >  7:119334/1‑149 (MQ=255)
                                                                                                                                            tGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGTCCTATAGCGCCGAAAGCTTGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCg        <  5:211049/148‑1 (MQ=255)
                                                                                                                                                 gCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACtt                     <  5:103812/130‑1 (MQ=255)
                                                                                                                                                 gCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACtt                     >  6:103812/1‑130 (MQ=255)
                                                                                                                                                 gCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCGTCAtct  <  1:255436/149‑1 (MQ=255)
                                                                                                                                                 gCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCGTCAt    <  2:335996/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGATCTCCAGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCGTCATCT  >  NC_002947/2168043‑2168336

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