Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,562,067 C→G 6.9% V155V (GTG→GTC) ‡ PP_3145 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,562,0670CG 6.9% 186.5 / 3.8 58V155V (GTG→GTC) ‡PP_3145hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/31);  new base G (2/2);  total (25/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 = 4.78e-01

TCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGA  >  NC_002947/3561937‑3562212
                                                                                                                                   |                                                                                                                                                 
tCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCgcgg                                                                                                                                 >  7:5055/1‑149 (MQ=255)
  acacCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaagg                                                                                                                                      <  6:2842/142‑4 (MQ=255)
  acacCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaagg                                                                                                                                      >  5:2842/1‑139 (MQ=255)
         ccGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGcc                                                                                                                        >  3:93573/1‑149 (MQ=255)
            cgATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCg                                                                                                                     >  2:159385/1‑149 (MQ=255)
               tAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGtt                                                                                                                   >  7:131758/1‑148 (MQ=255)
                 aGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                                <  6:162985/149‑1 (MQ=255)
                 aGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                                <  7:20793/149‑1 (MQ=255)
                 aGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGa                                                                                                                 <  8:132144/148‑1 (MQ=255)
                  gcgcGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAgg                                                                                                                           <  8:40139/137‑1 (MQ=255)
                  gcgcGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAgg                                                                                                                           >  7:40139/1‑137 (MQ=255)
                      gATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGC                                                                                                           <  1:159385‑M1/149‑1 (MQ=255)
                      gATAGGGATCACCTCGTAAGGCCGGAGGCCCTCGTTACTGGGGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGC                                                                                                           <  7:143815‑M1/149‑1 (MQ=255)
                           ggATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaaggg                                                                                                                                     <  8:170960/118‑5 (MQ=255)
                           ggATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAGCGTGATGCCGTTGGCCaaggg                                                                                                                                     >  7:170960/1‑114 (MQ=255)
                            gATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttat                                                                                                     >  5:1028‑M1/1‑145 (MQ=255)
                             aTCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGAt                                                                                                    <  4:93573/149‑1 (MQ=255)
                                aCCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCCCCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTc                                                                                                 >  2:318708/1‑149 (MQ=255)
                                 ccTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGt                                                                                                       >  3:271667/1‑142 (MQ=255)
                                 ccTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGt                                                                                                       <  4:271667/142‑1 (MQ=255)
                                  cTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCt                                                                                               >  1:216897/1‑149 (MQ=255)
                                  cTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCt                                                                                               <  7:11254/149‑1 (MQ=255)
                                      tAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCGACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTg                                                                                              >  1:241985/1‑147 (MQ=255)
                                       aaGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGc                                                                                          >  6:285762/1‑149 (MQ=255)
                                        aGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCGCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGc                                                                                          <  5:285762/148‑1 (MQ=255)
                                          gCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                                >  8:224920/1‑124 (MQ=255)
                                          gCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGaa                                                                                                                <  7:224920/124‑1 (MQ=255)
                                                 gCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCtgctgc                                                                                 >  7:28615/1‑148 (MQ=255)
                                                    cTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCG                                                                                                          >  4:252782‑M1/1‑120 (MQ=255)
                                                    cTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCG                                                                                                          <  3:252782‑M1/120‑1 (MQ=255)
                                                       gTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATa                                                                           <  8:5055/148‑1 (MQ=255)
                                                          aCTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAAccg                                                                       <  5:247132‑M2/149‑4 (MQ=255)
                                                            tggtggACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGAGAAccgcg                                                                     <  8:80284‑M2/149‑6 (MQ=255)
                                                              gtggACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGt                                                                   <  3:287060/149‑1 (MQ=255)
                                                               tggACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGtt                                                                  <  2:241985/149‑1 (MQ=255)
                                                                 gACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTgg                                                                <  3:153293/149‑1 (MQ=255)
                                                                  acacTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttatctcgcagcagacgctgcaggaaatcacccgctacttca                                                               <  6:1028‑M1/149‑43 (MQ=255)
                                                                      tGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttatctcgcagcagacgctgcaggaaatcac                                                                          >  1:92087‑M1/1‑103 (MQ=255)
                                                                      tGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGttatctcgcagcagacgctgcaggaaatcac                                                                          <  2:92087‑M1/134‑32 (MQ=255)
                                                                        ccGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCgg                                                          <  1:309204/148‑1 (MQ=255)
                                                                           tGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTg                                                      >  2:228375/1‑149 (MQ=255)
                                                                                gCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTa                                                 <  8:28615/149‑1 (MQ=255)
                                                                                 ccGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTgggg                                                   >  7:43944/1‑146 (MQ=255)
                                                                                 ccGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTgggg                                                   >  7:183866/1‑146 (MQ=255)
                                                                                   gggCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGaa                                              >  6:253517/1‑149 (MQ=255)
                                                                                    ggCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAAc                                             <  2:216897/149‑1 (MQ=255)
                                                                                       gcACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTg                                                      <  7:101645/137‑1 (MQ=255)
                                                                                       gcACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTg                                                      >  8:101645/1‑137 (MQ=255)
                                                                                                       ccGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGaa                            >  5:6853/1‑148 (MQ=255)
                                                                                                       ccGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAac                           <  8:131758/149‑1 (MQ=255)
                                                                                                         gggCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAacac                         <  8:183866/149‑1 (MQ=255)
                                                                                                             cAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTg                     <  1:228375/149‑1 (MQ=255)
                                                                                                                  cccTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTa                >  8:193783/1‑149 (MQ=255)
                                                                                                                   ccTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTAc               <  6:247132/149‑1 (MQ=255)
                                                                                                                         ccAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTg         <  8:43944/149‑1 (MQ=255)
                                                                                                                            aCGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGata      <  1:318708/149‑1 (MQ=255)
                                                                                                                            aCGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGata      <  6:6853/149‑1 (MQ=255)
                                                                                                                                cATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGa  >  7:196228/1‑149 (MQ=255)
                                                                                                                                   cACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTTATATTGa  >  2:217183/1‑146 (MQ=255)
                                                                                                                                   cACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGa  >  2:212564/1‑146 (MQ=255)
                                                                                                                                   cACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGa  <  1:217183/146‑1 (MQ=255)
                                                                                                                                   cACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGa  <  1:212564/146‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                                 
TCACACCTGCCGCGATAAGCGCGATAGGGATCACCTCGTAAGGCCTGAGGCCCTCGTTACTGGTGGACACTGCCGTGATCGCCGGGCGCACCACGATCACGCCCGGGCCAGGCCCCTTGGCCAACGGCATCACGCTGGCCATTTCGCGGCGCAGGGCCTCGTTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGA  >  NC_002947/3561937‑3562212

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