Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,792,250 A→G 5.5% S530S (AGT→AGC pvdI ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,792,2500AG 5.5% 214.9 / 2.3 73S530S (AGT→AGCpvdInon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (34/35);  new base G (2/2);  total (36/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 = 9.25e-01

CGCCCCCCCGGCCTTGAGAATCGCCAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGCAACAGGTGCGCCCGCCAGCGGCCAGTCACCCAGC  >  NC_002947/4792103‑4792397
                                                                                                                                                    |                                                                                                                                                   
cGCCCCCCCGGCCTTGAGAATCGCCAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAAc                                                                                                                                                    >  3:24068/1‑149 (MQ=255)
cGCCCCCCCGGCCTTGAGAATCGCCAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAAc                                                                                                                                                    <  2:87684/149‑1 (MQ=255)
   ccccccGGCCTTGAGAATCGCCAGCAGGCCAACCAACAGGTTCACCGAGGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACt                                                                                                                                                   >  1:257693/1‑148 (MQ=255)
       ccGGCCTTGAGAATCGCCAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTgagcga                                                                                                                                             <  4:24068/149‑1 (MQ=255)
       ccGGCCTTGAGAAGCGCCAGCAGGCCACCCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTgagcga                                                                                                                                             <  7:141074/149‑1 (MQ=255)
        cGGCCTTGAGAATCGCCAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTgagcgag                                                                                                                                            <  5:236792/149‑1 (MQ=255)
            cTTGAGAATCGCCAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCAt                                                                                                                                        >  6:335428/1‑149 (MQ=255)
                  aaTCGCCAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCa                                                                                                                                   >  8:289956/1‑148 (MQ=255)
                  aaTCGCCAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCAc                                                                                                                                  <  2:257693/149‑1 (MQ=255)
                    tCGCCAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCACGCCGTAACTGAGCGAGCCATCGGCACTg                                                                                                                                >  3:201666/1‑149 (MQ=255)
                      gCCAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGa                                                                                                                               <  4:201666/148‑1 (MQ=255)
                       ccAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTgagcgag                                                                                                                                            >  2:75874/1‑134 (MQ=255)
                       ccAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTgagcgag                                                                                                                                            <  1:75874/134‑1 (MQ=255)
                        cagcagGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACaa                                                                                                                            <  5:335428/149‑1 (MQ=255)
                           cagGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACgg                                                                                                                         <  2:175215/149‑1 (MQ=255)
                              gCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGa                                                                                                                      >  8:269479/1‑149 (MQ=255)
                               ccaaccaaCAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGAt                                                                                                                     >  3:300578/1‑149 (MQ=255)
                               ccaaccaaCAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGAt                                                                                                                     <  4:300578/149‑1 (MQ=255)
                                caaccaacAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATg                                                                                                                    >  3:111467/1‑149 (MQ=255)
                                caaccaacAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATg                                                                                                                    >  6:109328/1‑149 (MQ=255)
                                caaccaacAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATg                                                                                                                    >  8:264135/1‑149 (MQ=255)
                                     aacAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCgctg                                                                                                               <  4:23215/149‑1 (MQ=255)
                                          gTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCgctgcggct                                                                                                          >  2:321038/1‑149 (MQ=255)
                                              aCCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGCGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTgcgc                                                                                                      >  6:276523/1‑149 (MQ=255)
                                               ccGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTAAc                                                                                                                              >  2:281556/1‑124 (MQ=255)
                                               ccGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTAAc                                                                                                                              <  1:281556/124‑1 (MQ=255)
                                                 gAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCtt                                                                                                   <  7:269479/149‑1 (MQ=255)
                                                  aGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTg                                                                                                  <  5:109328/149‑1 (MQ=255)
                                                  aGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTg                                                                                                  <  7:264135/149‑1 (MQ=255)
                                                    gcgcTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTgagcga                                                                                                                                             >  7:49702/1‑105 (MQ=255)
                                                    gcgcTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTgagcga                                                                                                                                             <  8:49702/105‑1 (MQ=255)
                                                       cTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCaa                                                                                              <  7:289956/149‑1 (MQ=255)
                                                        tCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCa                                                                                               >  2:284393/1‑147 (MQ=255)
                                                          ggCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAGCTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAgg                                                                                           >  8:227118/1‑149 (MQ=255)
                                                          ggCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAGCTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAgg                                                                                           >  8:227119/1‑149 (MQ=255)
                                                                 aGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAGCTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCa                                                                                     <  7:227118/148‑1 (MQ=255)
                                                                 aGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAGATGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGGCa                                                                                     <  7:227119/148‑1 (MQ=255)
                                                                   gCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCAcgc                                                                                  <  4:111467/149‑1 (MQ=255)
                                                                      cACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCtt                                                                                                   <  8:19607/129‑1 (MQ=255)
                                                                      cACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCtt                                                                                                   >  7:19607/1‑129 (MQ=255)
                                                                               aTCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTgcagc                                                                       >  5:332693/1‑148 (MQ=255)
                                                                                  ggCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGTCAAAGGCCGTCACGCGCTGGTGCAGCAAc                                                                    <  1:321038/148‑1 (MQ=255)
                                                                                      aCCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGcc                                                               >  5:114738/1‑149 (MQ=255)
                                                                                          cacCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCgg                                                             >  8:37306/1‑147 (MQ=255)
                                                                                          cacCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCgg                                                             <  7:37306/147‑1 (MQ=255)
                                                                                           acCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCa                                                           >  5:321827/1‑148 (MQ=255)
                                                                                             cGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCAt                                                          >  8:144040/1‑147 (MQ=255)
                                                                                             cGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATaa                                                        <  6:332693/149‑1 (MQ=255)
                                                                                             cGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGAAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATaa                                                        <  1:284393/149‑1 (MQ=255)
                                                                                                    ttGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTgg                                                                                                 >  6:205710/1‑101 (MQ=255)
                                                                                                    ttGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTgg                                                                                                 <  5:205710/101‑1 (MQ=255)
                                                                                                        aGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCtt                                                                                                   >  7:328963/1‑95 (MQ=255)
                                                                                                        aGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCtt                                                                                                   <  8:328963/95‑1 (MQ=255)
                                                                                                          cTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCTTAACTGAGCGAGCGATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTTGGAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGTCATAATGCAGCACCtgcg                                           <  6:361886/149‑1 (MQ=255)
                                                                                                           tGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGgca                                                    <  1:277668/139‑1 (MQ=255)
                                                                                                           tGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGgca                                                    >  2:277668/1‑139 (MQ=255)
                                                                                                               cacacaGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCtgcggtgcg                                      <  6:114738/149‑1 (MQ=255)
                                                                                                                acacaGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCtgcggtgcgg                                     >  4:196980/1‑149 (MQ=255)
                                                                                                                 cacaGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGc                                    <  7:144040/149‑1 (MQ=255)
                                                                                                                  acaGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGc                                    >  2:240257/1‑148 (MQ=255)
                                                                                                                   caGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCtgcggtgc                                       <  5:299240/144‑1 (MQ=255)
                                                                                                                   caGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCtgcggtgc                                       >  6:299240/1‑144 (MQ=255)
                                                                                                                     gCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGCaa                                  >  7:198316/1‑147 (MQ=255)
                                                                                                                     gCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGCaa                                  >  5:159524/1‑147 (MQ=255)
                                                                                                                        gATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACgg                                                                  <  1:232948/112‑1 (MQ=255)
                                                                                                                        gATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACgg                                                                  >  2:232948/1‑112 (MQ=255)
                                                                                                                         aTCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGCAACAGGt                             <  6:321827/148‑1 (MQ=255)
                                                                                                                          tCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGCAACAGGTgc                           >  3:347036/1‑149 (MQ=255)
                                                                                                                                     gTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGCAACAGGTGCGCCCGCCAGCg                <  4:347036/149‑1 (MQ=255)
                                                                                                                                          gCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGCAACAGGTGCGCCCGCCAGCGGc              >  3:236407/1‑146 (MQ=255)
                                                                                                                                          gCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGCAACAGGTGCGCCCGCCAGCGGc              <  4:236407/146‑1 (MQ=255)
                                                                                                                                                 gTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCGCGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGCAACAGGTGCGCCCGCCAGCGGCCAGTCAcc      >  3:4265/1‑147 (MQ=255)
                                                                                                                                                 gTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTGGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGCAACAGGTGCGCCCGCCAGCGGCCAGTCACCca    <  5:276523/149‑1 (MQ=255)
                                                                                                                                                   aaCTGAGCGAGCCATCGGCACTGAAAACGGCGATGCGCTGCGGCTGCGCCTAGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCCGTGCGGCAACAGGTGCGCACGCCAGCGGCAAGTCACCcagc  >  7:181927/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCCCCCCCGGCCTTGAGAATCGCCAGCAGGCCAACCAACAGGTTCACCGAGCGCTCGGCGATCAGGCCCACCGGCTGATCGGCGACCACACCGTGCGCTTGCAGCTGGGCACACAGCTGATCGGCGCGCGCGTCGAGCTCGCCGTAACTGAGCGAGCCATCGGCACTGACAACGGCGATGCGCTGCGGCTGCGCCTTGGCAAAGGCCGTCACGCGCTGGTGCAGCAACGGGCCGGCATAAGGCAGCACCTGCGGTGCGGCAACAGGTGCGCCCGCCAGCGGCCAGTCACCCAGC  >  NC_002947/4792103‑4792397

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: