Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,821,134 G→C 5.8% L3352L (CTC→CTG pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,821,1340GC 5.8% 261.9 / 4.1 80L3352L (CTC→CTGpvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (34/41);  new base C (3/2);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCGCCAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCATCGCCTCGCCGC  >  NC_002947/4820986‑4821279
                                                                                                                                                      |                                                                                                                                                 
cTTCGCCAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTg                                                                                                                                                   >  4:56523/1‑149 (MQ=255)
     ccAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGc                                                                                                                                                >  7:116688/1‑147 (MQ=255)
      cAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAg                                                                                                                                             >  2:250190/1‑149 (MQ=255)
      cAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAg                                                                                                                                             >  4:208598/1‑149 (MQ=255)
         gggCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGa                                                                                                                                                  >  2:297639/1‑141 (MQ=255)
         gggCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGa                                                                                                                                                  <  1:297639/141‑1 (MQ=255)
                gTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGCCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCgg                                                                                                                                   <  6:59286/149‑1 (MQ=255)
                   gACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCt                                                                                                                                <  5:100127/149‑1 (MQ=255)
                        cGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCCGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTAGCTCAGCAAGGCAAAGCCGGCCTCccgc                                                                                                                           >  4:237827/1‑149 (MQ=255)
                         gaagaaATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGTGCAAGGCAAAGCCGGCCTcgcg                                                                                                                            <  7:55838/147‑1 (MQ=255)
                           agaaATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTcc                                                                                                                        <  3:233822/149‑1 (MQ=255)
                            gaaATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAGGGCAAAGCCGGCCTCGCGCTCCt                                                                                                                       >  2:236729/1‑149 (MQ=255)
                              aaTCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGc                                                                                                                     <  1:250190/149‑1 (MQ=255)
                                  gTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTtgt                                                                                                                 >  4:126611/1‑149 (MQ=255)
                                   ttCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTTGCTCAGCGAGGCaccgt                                                                                                                                      <  4:7116/127‑5 (MQ=255)
                                   ttCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTTGCTCAGCGAGGCaccgt                                                                                                                                      >  3:7116/1‑123 (MQ=255)
                                      atcaACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCGAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAg                                                                                                             >  5:112348/1‑149 (MQ=255)
                                      atcaACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAg                                                                                                             <  3:40643/149‑1 (MQ=255)
                                       tcaACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAgg                                                                                                            <  1:327739/149‑1 (MQ=255)
                                                   ggCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGATCCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGgcgc                                                                                                <  1:315115/149‑1 (MQ=255)
                                                     cACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACTGGCGCACCCGGCGCAGGTGGAACGGTCCCTCGCTGAGAAAGGCAAACCCGCCCTCGCGATCCTGATTGTGCAGGGCCTGCAGgcgc                                                                                                >  1:302806/1‑149 (MQ=255)
                                                     cACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTcgc                                                                                                                             <  8:67612/118‑1 (MQ=255)
                                                     cACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTcgc                                                                                                                             >  7:67612/1‑118 (MQ=255)
                                                      aCCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGCGCCTGCAGGCcctgc                                                                                             <  3:208598/149‑1 (MQ=255)
                                                       ccAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAGGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCgctgct                                                                                            <  1:236729/149‑1 (MQ=255)
                                                        cAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTc                                                                                           >  5:14085/1‑149 (MQ=255)
                                                        cAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTc                                                                                           >  3:347257/1‑149 (MQ=255)
                                                          ggCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCt                                                                                         <  5:135146/149‑1 (MQ=255)
                                                            cATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTgg                                                                                       <  3:126611/149‑1 (MQ=255)
                                                              tCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTggg                                                                                      >  5:52077/1‑148 (MQ=255)
                                                                    aGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCAtcgt                                                                                <  7:273190/148‑1 (MQ=255)
                                                                      gATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAg                                                                                                                                             >  7:264448/1‑85 (MQ=255)
                                                                      gATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAg                                                                                                                                             <  8:264448/85‑1 (MQ=255)
                                                                       aTGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCa                                                                                                     <  4:152175/124‑1 (MQ=255)
                                                                       aTGTGGTGGTTGCTCATCATGAACCAGTAGAGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCa                                                                                                     >  3:152175/1‑124 (MQ=255)
                                                                            gtggtTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAgg                                                                        <  3:56523/148‑1 (MQ=255)
                                                                                ttGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGc                                                                                                               <  6:42780/105‑1 (MQ=255)
                                                                                ttGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGc                                                                                                               >  5:42780/1‑105 (MQ=255)
                                                                                    tcatcaTGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCAGCGCCAg                                                               >  4:114728/1‑149 (MQ=255)
                                                                                      atcatGATCCAGTTGCGTTCCTCGGCCAC‑‑GCGCTCCAGGCGCAGGTGGAACGGGGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGgcgc                                                                                                <  6:112401/114‑1 (MQ=255)
                                                                                      atcatGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGt                                                              >  4:137179/1‑148 (MQ=255)
                                                                                      atcatGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGt                                                              <  3:137179/148‑1 (MQ=255)
                                                                                      atcatGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGCTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGgcgc                                                                                                >  5:112401/1‑114 (MQ=255)
                                                                                          tGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTgg                                                         >  1:169771/1‑149 (MQ=255)
                                                                                          tGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTgg                                                         <  5:308277/149‑1 (MQ=255)
                                                                                          tGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTgg                                                         <  2:302806/149‑1 (MQ=255)
                                                                                           gAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTgg                                                         >  7:252667/1‑148 (MQ=255)
                                                                                            aaCCAGTAGCGTTCCTGGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTa                                                       <  3:237827/149‑1 (MQ=255)
                                                                                            aaCCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTgg                                                         >  7:220087/1‑147 (MQ=255)
                                                                                            aaCCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTgg                                                         >  5:287986/1‑147 (MQ=255)
                                                                                               cAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGt                                                     >  2:164890/1‑148 (MQ=255)
                                                                                               cAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTc                                                    >  3:162817/1‑149 (MQ=255)
                                                                                               cAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTc                                                    <  8:116688/149‑1 (MQ=255)
                                                                                                    gCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCg                                               >  2:256851/1‑149 (MQ=255)
                                                                                                       ttCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGc                                              >  2:219284/1‑147 (MQ=255)
                                                                                                         ccTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAg                                           <  1:219284/148‑1 (MQ=255)
                                                                                                         ccTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAg                                           <  8:252667/148‑1 (MQ=255)
                                                                                                          cTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAgg                                          <  4:347257/148‑1 (MQ=255)
                                                                                                                     gcgcACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCt                                >  3:74215/1‑149 (MQ=255)
                                                                                                                     gcgcACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCt                                >  5:345504/1‑149 (MQ=255)
                                                                                                                       gcACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTtg                              <  4:162817/149‑1 (MQ=255)
                                                                                                                       gcACCAGGCGCAGGTGCAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTtg                              >  4:253695/1‑149 (MQ=255)
                                                                                                                        cACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTtgt                             <  6:52077/149‑1 (MQ=255)
                                                                                                                           cAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGga                          <  1:256851/149‑1 (MQ=255)
                                                                                                                            aGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGgat                         <  2:169771/149‑1 (MQ=255)
                                                                                                                                  aGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCt                                <  4:183585/136‑1 (MQ=255)
                                                                                                                                  aGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCt                                >  3:183585/1‑136 (MQ=255)
                                                                                                                                  aGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTg                   <  6:345504/149‑1 (MQ=255)
                                                                                                                                  aGGTGCAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTg                   <  3:253695/149‑1 (MQ=255)
                                                                                                                                      ggAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCAtc                                                                                  <  4:303871/82‑1 (MQ=255)
                                                                                                                                      ggAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCAtc                                                                                  >  3:303871/1‑82 (MQ=255)
                                                                                                                                           tttGCCTTGCTCAGCGAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATGGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGGCAGGTTGCCCGGCTTGTGga                          <  5:209051/131‑1 (MQ=255)
                                                                                                                                           tttGCCTTGCTCAGCGAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTCCCCGGCATGTGga                          >  6:209051/3‑133 (MQ=255)
                                                                                                                                             tGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCa                                            <  6:377433/113‑1 (MQ=255)
                                                                                                                                             tGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCa                                            >  5:377433/1‑113 (MQ=255)
                                                                                                                                               ccTCGCTGAGCGAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCAtcgcctc       <  6:112348/148‑1 (MQ=255)
                                                                                                                                                   gCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGTTTGACCGGATTGTGGATGATCTGCAGCATCGCCTcgccgc  <  1:164890/149‑1 (MQ=255)
                                                                                                                                                   gCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCAtc            >  1:124393/1‑139 (MQ=255)
                                                                                                                                                   gCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCAtc            <  2:124393/139‑1 (MQ=255)
                                                                                                                                                    cTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGt                                                              >  6:98030/1‑88 (MQ=255)
                                                                                                                                                    cTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGt                                                              <  5:98030/88‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CTTCGCCAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCAC‑‑GCGCACCAGGCGCAGGTGGAACGGTGCCTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCATCGCCTCGCCGC  >  NC_002947/4820986‑4821279

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: