Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,168,192 A→T 5.5% intergenic (+92/+596) 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,1920AT 5.5% 226.6 / 2.5 73intergenic (+92/+596)PP_1921/PP_1922hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/39);  new base T (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 = 4.32e-01

AGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCGTCATCT  >  NC_002947/2168051‑2168336
                                                                                                                                             |                                                                                                                                                
agagGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCa                                                                                                                                           <  2:98034/149‑1 (MQ=255)
   gtgACAAAACGACTAGCATCTATTTTTGTGTGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCTCTGCTAGCAcc                                                                                                                                         <  8:157210/146‑1 (MQ=255)
   gggACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCg                                                                                                                                        >  5:250463/1‑149 (MQ=255)
        aaaaCGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCtttt                                                                                                                                   <  4:134157/149‑1 (MQ=255)
        aaaaCGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCtttt                                                                                                                                   <  3:152671/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                                                                                                                                    <  5:71327/135‑4 (MQ=255)
                     tCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCTCTGCTAGCACCGcatc                                                                                                                                    >  6:71327/1‑132 (MQ=255)
                         tttttGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCgg                                                                                                                    <  5:76152/147‑1 (MQ=255)
                         tttttGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCg                                                                                                                     <  8:211690/146‑1 (MQ=255)
                         tttttGGGGGCTCTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCggg                                                                                                                   >  4:188086/1‑148 (MQ=255)
                             tGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCggg                                                                                                                   <  3:340757/144‑1 (MQ=255)
                             tGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCggg                                                                                                                   >  4:340757/1‑144 (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                                                                                                             >  8:238471/1‑141 (MQ=255)
                                      tATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAgct                                                                                                             <  7:238471/141‑1 (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                                                                                                         >  5:362570/1‑132 (MQ=255)
                                                   tACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGctgctg                                                                                                         <  6:362570/132‑1 (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                                                          <  1:123062/149‑1 (MQ=255)
                                                                                 aTTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAAt                                                          <  3:88358/149‑1 (MQ=255)
                                                                                 aTTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAAt                                                          <  6:250463/149‑1 (MQ=255)
                                                                                 aTTAAGTAATATATTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAAt                                                          <  3:83105/149‑1 (MQ=255)
                                                                                    aaGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGt                                                        >  5:239608/1‑148 (MQ=255)
                                                                                    aaGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGt                                                        >  5:46457/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                                                                         >  2:374491/1‑124 (MQ=255)
                                                                                           atCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTg                                                                         <  1:374491/124‑1 (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                                                              >  4:330367/1‑114 (MQ=255)
                                                                                                                gCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCtctc                                                              <  3:330367/114‑1 (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)
                                                                                                                                             |                                                                                                                                                
AGAGGGACAAAACGAATAGCATCTATTTTTGGGGGCTTTATGCTGCGTGACTACGCGGCATAGAGCTACCGAAAAGCGCGAATTAAGTAATATCTTATTTGGGGCCGTATGTGCTTTGAAAGTCAGATGCGGTGCTAGCAGAGCTAGCACCGCTTTTGGCATTGTGTCTCGGGCAAGCTGCTGGTGACCTATAGCGCCGAAAGCTGGGTGCTTTGGAGTCGCTCTCGAATGTACAACAGTTGGCACCGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCGTCATCT  >  NC_002947/2168051‑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: ‑