Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,898,939 G→A 9.6% A67V (GCC→GTC)  PP_3441 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,898,9390GA 9.6% 172.0 / 6.4 62A67V (GCC→GTC) PP_3441hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base A (3/3);  total (29/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 = 9.98e-01

CACGGCGCTCAGGTAGTGCGGCACTGACAGCGACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACTTTGCTGGCCAGCAGCCTGATTTCG  >  NC_002947/3898794‑3899082
                                                                                                                                                 |                                                                                                                                               
cacGGCGCTCAGGTAGTGCGGCACTGACAGCGACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGACCa                                                                                                                                              >  8:222093/1‑149 (MQ=255)
        tCAGGTAGTGCGGCACTGACAGCGACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGACCATTTCCCGt                                                                                                                                      >  5:23598/1‑149 (MQ=255)
        tCAGGTAGTGCGGCACTGACAGCGACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGACCATTTCCCGt                                                                                                                                      >  3:294057/1‑149 (MQ=255)
          aGGTAGTGCGGCACTGACAGCGACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTtg                                                                                                                                    >  7:221784/1‑149 (MQ=255)
          aGGTAGTGCGGCACTGACAGCGACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTtg                                                                                                                                    <  4:353176/149‑1 (MQ=255)
            gTAGTGCGGCACTGACAGCGACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTtgct                                                                                                                                  >  4:364874/1‑149 (MQ=255)
                    gCACTGACAGCGACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGGGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCgg                                                                                                                          <  8:327790/149‑1 (MQ=255)
                          aCAGCGACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTCGCTGCTTGCGGCCGa                                                                                                                      <  5:118772/147‑1 (MQ=255)
                            aGCGACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATg                                                                                                                    >  8:124297/1‑147 (MQ=255)
                               gACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTagag                                                                                                               >  7:312712/1‑149 (MQ=255)
                               gACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTagag                                                                                                               >  8:279774/1‑149 (MQ=255)
                               gACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGACCATTTCCCGTTGCTGCTTGCGGCCGATGTagag                                                                                                               <  7:222093/149‑1 (MQ=255)
                                   ccGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGAtt                                                                                                           <  3:102294/149‑1 (MQ=255)
                                     gACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTc                                                                                                         <  3:364874/149‑1 (MQ=255)
                                            gCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCTGATGTAGAGGATTTCGTACggg                                                                                                  <  6:186590/149‑1 (MQ=255)
                                              aTTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTc                                                                                                         >  1:26351/1‑140 (MQ=255)
                                              aTTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTc                                                                                                         <  2:26351/140‑1 (MQ=255)
                                               tttCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCCCTTGCACTGTCGACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGctc                                                                                               >  2:237941/1‑149 (MQ=255)
                                                   gCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCt                                                                                           >  4:46938/1‑149 (MQ=255)
                                                       atgaCCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTg                                                                                       >  3:347941/1‑149 (MQ=255)
                                                        tgaCCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTgc                                                                                      <  2:268105/149‑1 (MQ=255)
                                                          aCCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTgcg                                                                                     <  4:280458/148‑1 (MQ=255)
                                                          aCCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTgcg                                                                                     <  7:133284/148‑1 (MQ=255)
                                                            cACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTgcgggc                                                                                  >  5:339582/1‑149 (MQ=255)
                                                            cACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGAGTTCGTACGGGCTCTCCTGTTgcgggc                                                                                  <  3:46938/149‑1 (MQ=255)
                                                                    ggctggcAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTg                                                                          >  5:211896/1‑149 (MQ=255)
                                                                     gctggcAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTgg                                                                         >  4:174786/1‑149 (MQ=255)
                                                                      ctggcAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGc                                                                        <  5:189247/149‑1 (MQ=255)
                                                                        ggcAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGa                                                                      <  8:312712/149‑1 (MQ=255)
                                                                           aGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACtt                                                                   <  2:262611/149‑1 (MQ=255)
                                                                           aGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACtt                                                                   <  1:350976/149‑1 (MQ=255)
                                                                              cgccgTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTccc                                                                <  6:211896/149‑1 (MQ=255)
                                                                                 cgTGGCACTTGCACTGTCCACGATCGAAGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAAc                                                             <  7:176350/149‑1 (MQ=255)
                                                                                  gTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACtt                                                                   >  5:76010/1‑142 (MQ=255)
                                                                                  gTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACtt                                                                   <  6:76010/142‑1 (MQ=255)
                                                                                      cACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCag                                                        <  7:279774/149‑1 (MQ=255)
                                                                                        cTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCagag                                                      <  3:174786/149‑1 (MQ=255)
                                                                                        cTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCa                                                         >  7:281470/1‑146 (MQ=255)
                                                                                         ttGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAgaga                                                     >  2:246733/1‑149 (MQ=255)
                                                                                             aCTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAgg                                                 >  7:27584/1‑149 (MQ=255)
                                                                                               tGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCa                                               >  5:338843/1‑149 (MQ=255)
                                                                                                   cACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCCGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCAcgcg                                           >  5:306242/1‑149 (MQ=255)
                                                                                                    aCGCTCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGACCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACgcgc                                          <  4:294057/149‑1 (MQ=255)
                                                                                                    aCGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGACCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACgcgc                                          <  6:23598/149‑1 (MQ=255)
                                                                                                     cGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACgcgc                                          >  3:138296/1‑148 (MQ=255)
                                                                                                          gAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTCGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCGTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACAt                                    >  2:310183/1‑149 (MQ=255)
                                                                                                              cgcTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTtct                                <  8:281470/149‑1 (MQ=255)
                                                                                                               gcTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTtctc                               <  6:339582/149‑1 (MQ=255)
                                                                                                                cTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGTAACCTGGCGACTTCCCTACGCCAGAGACAGGCACGCGCCAACATTTCTCg                              <  1:237941/149‑1 (MQ=255)
                                                                                                                     cTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACtt                         >  5:384227/1‑149 (MQ=255)
                                                                                                                          ccTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACTTTGCTg                    <  4:347941/149‑1 (MQ=255)
                                                                                                                            ttGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTtctc                               <  6:268906/136‑1 (MQ=255)
                                                                                                                            ttGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTtctc                               >  5:268906/1‑136 (MQ=255)
                                                                                                                               ccACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACTTTGCTg                    >  8:41254/1‑144 (MQ=255)
                                                                                                                               ccACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACTTTGCTg                    <  7:41254/144‑1 (MQ=255)
                                                                                                                                cACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACTTTGCTGGCcagc              >  3:32734/1‑149 (MQ=255)
                                                                                                                                 aCCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACTTTGCTGGCcagca             <  6:338843/149‑1 (MQ=255)
                                                                                                                                 aCCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACTTTGCTGGCcag               >  2:186953/1‑147 (MQ=255)
                                                                                                                                   cGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACTTTGCTGGCCagcagc           <  4:138296/149‑1 (MQ=255)
                                                                                                                                    gATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACTTTGCTGGCCAGCAGcc          <  7:124297/149‑1 (MQ=255)
                                                                                                                                         gcagcTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACTTTGCTGGCCAGCAGCCTGAtt     <  1:246733/149‑1 (MQ=255)
                                                                                                                                             cTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACTTTGCTGGCCAGCAGCCTGATTTCg  >  7:231844/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CACGGCGCTCAGGTAGTGCGGCACTGACAGCGACCCCGACGTGGGCATTTCGCTGATGACCACGACGTGGCTGGCAGCCGCCGTGGCACTTGCACTGTCCACGATCGAGGCGCTGCCCTGGCCCTTGCCACCGATCAGCAGCTTGGCCATTTCCCGTTGCTGCTTGCGGCCGATGTAGAGGATTTCGTACGGGCTCTCCTGTTGCGGGCGGAACCTGGCGACTTCCCAACGCCAGAGACAGGCACGCGCCAACATTTCTCGCACTTTGCTGGCCAGCAGCCTGATTTCG  >  NC_002947/3898794‑3899082

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