Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,127,636 Δ1 bp 6.4% coding (1189/2187 nt) PP_4514 → alpha/beta hydrolase superfamily esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,127,6360T. 6.4% 196.9 / 7.0 62coding (1189/2187 nt)PP_4514alpha/beta hydrolase superfamily esterase
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base . (2/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCTTGGCGCCGAATGGC  >  NC_002947/5127493‑5127783
                                                                                                                                               |                                                                                                                                                   
atgTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGt                                                                                                                                                 <  3:322828/148‑1 (MQ=255)
   tgATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGAt                                                                                                                                              <  1:329274/148‑1 (MQ=255)
    gATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTgg                                                                                                                                                  <  5:38597/143‑1 (MQ=255)
    gATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTgg                                                                                                                                                  >  6:38597/1‑143 (MQ=255)
     aTCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCa                                                                                                                                            <  7:258146/148‑1 (MQ=255)
          cTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGccc                                                                                                                                      <  1:208516/149‑1 (MQ=255)
                  tCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGa                                                                                                                              <  2:86489/149‑1 (MQ=255)
                     gCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGAcgc                                                                                                                           <  6:164606/149‑1 (MQ=255)
                     gCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGAcgc                                                                                                                           >  2:200492/1‑149 (MQ=255)
                          ccGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGcggcgg                                                                                                                      <  1:200492/149‑1 (MQ=255)
                             aTCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCAc                                                                                                                   >  6:336527/1‑149 (MQ=255)
                             aTCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCAc                                                                                                                   <  5:151052/149‑1 (MQ=255)
                                    cACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTAc                                                                                                            <  3:17375/149‑1 (MQ=255)
                                                tGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGc                                                                                                >  5:310048/1‑149 (MQ=255)
                                                tGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGc                                                                                                >  4:109177/1‑149 (MQ=255)
                                                    ccGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGa                                                                                            <  3:72353/149‑1 (MQ=255)
                                                      gCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAAc                                                                                          <  5:336527/149‑1 (MQ=255)
                                                          gCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGatgcaag                                                                                      >  1:95830/1‑143 (MQ=255)
                                                          gCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGt                                                                                      >  2:334530/1‑149 (MQ=255)
                                                            aGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTcc                                                                                    >  6:152541/1‑149 (MQ=255)
                                                                   cATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGCCCTGACCATGTCCGATGAcc                                                                             >  7:28385/1‑149 (MQ=255)
                                                                      gggTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAAc                                                                                          <  7:125443/133‑1 (MQ=255)
                                                                      gggTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAAc                                                                                          >  8:125443/1‑133 (MQ=255)
                                                                         tACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGcgcg                                                                       >  8:54892/1‑149 (MQ=255)
                                                                                  tGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAgc                                                               <  3:109177/148‑1 (MQ=255)
                                                                                    gCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAgcgc                                                             <  5:149563/148‑1 (MQ=255)
                                                                                      cTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCAt                                                           >  8:306003/1‑148 (MQ=255)
                                                                                      cTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCAt                                                           <  7:306003/148‑1 (MQ=255)
                                                                                         gACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCaa                                                       <  2:78944/149‑1 (MQ=255)
                                                                                            tactTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCt                                                                                        >  2:246601/1‑113 (MQ=255)
                                                                                              ctTCGACAGGGTGCAGTACCGCGTGTTGCACAAGGTCAACGCCCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCGGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTgg                                                  <  8:186303/149‑1 (MQ=255)
                                                                                              ctTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTg                                                   <  6:310048/148‑1 (MQ=255)
                                                                                               ttCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTggg                                                 >  4:312471/1‑149 (MQ=255)
                                                                                                tCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGgc                                                >  5:120021/1‑149 (MQ=255)
                                                                                                 cGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTggg                                                 <  7:54892/147‑1 (MQ=255)
                                                                                                  gACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGgcgc                                              <  3:312471/149‑1 (MQ=255)
                                                                                                  gACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGgcgc                                              <  6:120021/149‑1 (MQ=255)
                                                                                                            aGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCgg                                    >  6:237522/1‑149 (MQ=255)
                                                                                                            aGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCgg                                    <  5:237522/149‑1 (MQ=255)
                                                                                                                ccGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTc                                 >  8:227306/1‑148 (MQ=255)
                                                                                                                   cgTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACg                             <  1:334530/149‑1 (MQ=255)
                                                                                                                     tggtgCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGc                           <  5:23060/149‑1 (MQ=255)
                                                                                                                            caAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCg                                     >  5:177197/1‑132 (MQ=255)
                                                                                                                            caAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCg                                     <  6:177197/132‑1 (MQ=255)
                                                                                                                             aaGGGCAACGACCACAGCTTGGTGATCAAGGCCCGCGGCCGACGCGGCGGCACCGTCTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAg                   <  4:27777/149‑1 (MQ=255)
                                                                                                                             aaGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAg                   >  6:325363/1‑149 (MQ=255)
                                                                                                                             aaGGGCAACGACCACACCTTGGTGATCAACGCACGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCaa                      >  3:27777/1‑146 (MQ=255)
                                                                                                                             aaGGGCAACGACCACACC‑TGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGc                  >  5:13788/1‑149 (MQ=255)
                                                                                                                             aaGGGCAACGACCACACC‑TGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGc                  >  5:8877/1‑149 (MQ=255)
                                                                                                                               gggCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGcc                 >  6:273933/1‑149 (MQ=255)
                                                                                                                               gggCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGcc                 >  1:293575/1‑149 (MQ=255)
                                                                                                                                ggCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCt                <  7:227306/149‑1 (MQ=255)
                                                                                                                                 gCAACGACCACACC‑TGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCTTg              <  6:13788/149‑1 (MQ=255)
                                                                                                                                 gCAACGACCACACC‑TGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCTTg              <  6:8877/149‑1 (MQ=255)
                                                                                                                                  cAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCTTg              >  7:127166/1‑149 (MQ=255)
                                                                                                                                         cacaCCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCTTgg             >  1:238649/1‑143 (MQ=255)
                                                                                                                                         cacaCCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCTTgg             <  2:238649/143‑1 (MQ=255)
                                                                                                                                          acacCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAATCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCTTGGCGCCGaa      >  2:283120/1‑149 (MQ=255)
                                                                                                                                            acCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCTTg              >  5:226654/1‑139 (MQ=255)
                                                                                                                                            acCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCTTg              <  6:226654/139‑1 (MQ=255)
                                                                                                                                            acCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCTTGGCGCCGAATg    <  2:293575/149‑1 (MQ=255)
                                                                                                                                              cTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAATCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCTTGGCGCCGAATGGc  <  1:283120/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
ATGTGATCCGCTCGTATATCCGCCAGCCGATCGACGCACCGCTGGAACTGGACCGCCTGCAGACCGACATGGGTACGCTGTATGGCCTGGACTACTTCGACAGGGTGCAGTACCGCGTGGTGCACAAGGGCAACGACCACACCTTGGTGATCAACGCCCGCGGCCGACGCGGCGGCACCGACTACCTGCGCCTGGGCCTGAACCTGTCCGATGACCTGCGCGGCGACAGCGCATTCAACCTGGGCGCCAGCTACCGGGTCAACGGCATCAACAGCCTTGGCGCCGAATGGC  >  NC_002947/5127493‑5127783

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