Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 924,320 C→G 6.7% A291A (GCG→GCC) ‡ PP_0804 ← protein secretion ABC transporter permease/ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947924,3200CG 6.7% 209.1 / 5.0 71A291A (GCG→GCC) ‡PP_0804protein secretion ABC transporter permease/ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base G (2/3);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

CCGTGCCAGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCACCCGCAGGTCGGC  >  NC_002947/924179‑924460
                                                                                                                                             |                                                                                                                                            
ccGTGCCAGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGc                                                                                                                                       <  3:2878/149‑1 (MQ=255)
ccGTGCCAGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGc                                                                                                                                       >  2:275858/1‑149 (MQ=255)
     ccAGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGTCCAGCACCGCGACGGTGCTggcc                                                                                                                                  <  2:179656/149‑3 (MQ=255)
       aGGCGCCTTTGTGCAAGCAACCTGGGCAGCACCATCGCCAGCAGCGCGACCAGAGGGATGAACACCAGTACGCCGCCCATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGACGGCCAGCACCCTCCCAGTGCGGGACGt                                                                                                                                >  6:203793/1‑149 (MQ=255)
        ggCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTg                                                                                                                               <  3:117725/149‑1 (MQ=255)
              tttGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCACCACCGTCGCGGTGCTGGACGGGATCAgg                                                                                                                         >  5:195235/1‑149 (MQ=255)
                tgtgCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGt                                                                                                                        <  8:21998/148‑1 (MQ=255)
                    caagcaaCCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTcgcg                                                                                                                    >  2:56946/1‑148 (MQ=255)
                     aagcaaCCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGa                                                                                                                  >  6:153566/1‑149 (MQ=255)
                       gcaaCCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGa                                                                                                                  >  7:116133/1‑147 (MQ=255)
                       gcaaCCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGa                                                                                                                  <  8:116133/147‑1 (MQ=255)
                       gcaaCCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATc                                                                                                                <  5:153566/149‑1 (MQ=255)
                            ccGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGGCGTGATCAGGTCGCGGATCGACTc                                                                                                           <  2:349951/149‑1 (MQ=255)
                             cGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTcc                                                                                                          <  6:328219/149‑1 (MQ=255)
                             cGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTcc                                                                                                          <  1:56946/149‑1 (MQ=255)
                              gggCAGCACCATGGCCAGGAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCa                                                                                                         <  8:135361/149‑1 (MQ=255)
                                     aCCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTcgcg                                                                                                  >  6:150051/1‑149 (MQ=255)
                                      ccATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCgcgc                                                                                                 <  1:109090/149‑1 (MQ=255)
                                       cATGGCCAGCATCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGa                                                                                                              <  7:325224/135‑1 (MQ=255)
                                         tGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTcgcagt                                                                                                >  1:131770/1‑144 (MQ=255)
                                            ccAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCCCCGTCGCGGTGGTGGACGTGATCAGGGCGCGGATCGACTCCAGCTCGCGCAACTgg                                                                                           >  2:230091/1‑148 (MQ=255)
                                            ccAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGc                                                                                           <  5:150051/149‑1 (MQ=255)
                                             cagcagCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGCCTCCAGCTCGCGCAACTGc                                                                                           >  7:292612/1‑148 (MQ=255)
                                                cagcGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCgaga                                                                                       <  2:40846/149‑1 (MQ=255)
                                                         cAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGc                                                                              >  8:285941/1‑149 (MQ=255)
                                                          aGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGaa                                                                                   <  2:261562/143‑1 (MQ=255)
                                                          aGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGaa                                                                                   >  1:261562/1‑143 (MQ=255)
                                                           gCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCGGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCg                                                                            >  8:19378/1‑149 (MQ=255)
                                                                   aaCACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACtt                                                                    >  5:103277/1‑149 (MQ=255)
                                                                     cacCAGTCCGCCGCGGAGCAGCCCGAACACGAACAGGAACAGCAGGAAGAACGGCAGGGCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGTTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTgg                                                                  <  5:97947/149‑1 (MQ=255)
                                                                      acCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGACCGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTggg                                                                 <  7:285941/149‑1 (MQ=255)
                                                                        cAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGcc                                                               <  1:275858/149‑1 (MQ=255)
                                                                        cAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGcc                                                               >  1:58668/1‑149 (MQ=255)
                                                                               ccgccgATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCACCACCGTCGCGGTGGAGGACGTGATCAGGGGCCCGAAAGCCAACAGCAACCCCAACGGACAGAAGAACGAGCCCGAGGACAAGGGCCGCAcc                                                          >  4:155072/1‑147 (MQ=255)
                                                                                  ccgATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGtt                                                     <  6:103277/149‑1 (MQ=255)
                                                                                     aTCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTgcg                                                  >  3:337469/1‑149 (MQ=255)
                                                                                         gCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTgcgcagg                                              >  4:201323/1‑149 (MQ=255)
                                                                                         gCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTgcg                                                  <  4:136903/145‑1 (MQ=255)
                                                                                         gCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTgcg                                                  >  3:136903/1‑145 (MQ=255)
                                                                                           cAGAACACGAGCAGGAACAGCAGGAAGAACGGCAGGTCGGCCACCCCCGTCGCGGTGGGGGGCGTGATCAGGTCGCGGATGGACGCCAGCACGCGCCACAGCGAGATGAACGAGCCGGGGGCCGTTGGCCGCCCCGCGTTgcgccggcg                                            >  4:187485/1‑149 (MQ=255)
                                                                                             gAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTgcgcaggcgca                                          <  2:58668/149‑1 (MQ=255)
                                                                                               acacGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCGGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGcgcaggcgcagg                                        <  7:19378/149‑1 (MQ=255)
                                                                                                  cGAAGAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACGGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCa                                     <  3:155072/149‑1 (MQ=255)
                                                                                                   gaacaggaacagCAGGAAGAACGGCAGGTCGGCCAGCACCCTCGCGGTGCTGGACGTGATCAGGTCGCGGATAGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTTGGCCGCACCGAGTTGCGCAGGCGCAGGGCAt                                    >  4:70291/1‑149 (MQ=255)
                                                                                                     acatgaacagCAGTAAGAACGGCTGGTCGTCCAGCACCGTGGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCTCAACTTCGAGATTAACGAGCCGGTGGACTTTTGCCGCACCGAGTTGCGCAGGCGCAGGGCATgg                                  <  5:203793/149‑1 (MQ=255)
                                                                                                        ggaacagCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCg                               >  8:169436/1‑149 (MQ=255)
                                                                                                         gaacagCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATg                                                                                     >  7:67010/1‑94 (MQ=255)
                                                                                                         gaacagCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATg                                                                                     <  8:67010/94‑1 (MQ=255)
                                                                                                          aacagCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTa                             <  7:169436/149‑1 (MQ=255)
                                                                                                           ccatcagGAAGAACGGCAGGGCGGCCAGCACCGTCTCGGTGCTGGACGTGATCAGGTCGCGGATCGAGTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGTTGGACTTGGGCCTCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTa                             <  6:195235/144‑1 (MQ=255)
                                                                                                               cagGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGAc                          >  6:184121/1‑147 (MQ=255)
                                                                                                                agGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATc                                                                                                                <  6:330399/60‑1 (MQ=255)
                                                                                                                agGAAGAACGGCAGGTCGGCCACCACCGTCGCGGTGCATGACGTGATCAGGTCGCGGAta                                                                                                                >  5:330399/1‑59 (MQ=255)
                                                                                                                 ggAAGAACGGCAGGTCGGCCAGCACCGGCGCGGGGCTGGACGTGATCAGGTCGCGGATCGCCTACAGCGCGCGCAACTTCGAGATGAACGAGCCGGTGGACGTGGGCCGCCCCCAGTTGCGCAGGCGCCGGGCATGGCCGTAGACCAgg                      >  4:212230/1‑149 (MQ=255)
                                                                                                                  gaagaaCGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATg                                   <  1:186665/135‑1 (MQ=255)
                                                                                                                       aCGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTa                             >  5:217978/1‑136 (MQ=255)
                                                                                                                       aCGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTa                             <  6:217978/136‑1 (MQ=255)
                                                                                                                       aCGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTc                <  3:187485/149‑1 (MQ=255)
                                                                                                                         ggCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGTGCAGGGCATGGCCGTAGACCTGGTCGCTCAc              >  4:113983/1‑149 (MQ=255)
                                                                                                                         ggCAGGTCGGCCAGCACCGTCGCGGTGCGGGACGTGATCAGGTCGCGGGTGGACTCCACCTCGCGCAACAGGGAGATGAAAGAGCCCGTGGGCATGGGGCGCACCGAGTTGCGCAGGCGCAGGGCCAGGCCGGAGACCAGGTCGCGCAc              >  4:49009/1‑149 (MQ=255)
                                                                                                                           cAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCAccc            <  1:230091/149‑1 (MQ=255)
                                                                                                                              gTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCACCCGCa         <  3:201323/149‑1 (MQ=255)
                                                                                                                               acgtccAGCACCGCGACGGTGCTGGCCGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGcgcaggcgcagg                                        <  3:75114/113‑1 (MQ=255)
                                                                                                                               acgtccAGCACCGCGACGGTGCTGGCCGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGcgcaggcgcagg                                        >  4:75114/5‑117 (MQ=255)
                                                                                                                                  gCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCg                                                         >  4:57587/1‑97 (MQ=255)
                                                                                                                                  gCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCg                                                         <  3:57587/97‑1 (MQ=255)
                                                                                                                                  gCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCAcc             <  2:275877/141‑1 (MQ=255)
                                                                                                                                  gCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCAcc             >  1:275877/1‑141 (MQ=255)
                                                                                                                                   ccAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCACCCGCAGGTCg    >  3:270128/1‑149 (MQ=255)
                                                                                                                                    cAGCACCGCGACGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCACCCGCAGGTCgg   >  2:288652/1‑149 (MQ=255)
                                                                                                                                     aGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCACCCGCAGGTCGgc  <  8:351222/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                            
CCGTGCCAGGCGCCTTTGTGCAAGCAACCCGGGCAGCACCATGGCCAGCAGCGCGACCAGCGGGATGAACACCAGTACGCCGCCGATCAGCCAGAACACGAACAGGAACAGCAGGAAGAACGGCAGGTCGGCCAGCACCGTCGCGGTGCTGGACGTGATCAGGTCGCGGATCGACTCCAGCTCGCGCAACTGCGAGATGAACGAGCCGGTGGACTTGGGCCGCACCGAGTTGCGCAGGCGCAGGGCATGGCCGTAGACCAGGTCGCTCACCCGCAGGTCGGC  >  NC_002947/924179‑924460

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