Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,773,119 T→C 9.7% T2050A (ACC→GCC) ‡ pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,773,1190TC 9.7% 192.1 / 9.5 72T2050A (ACC→GCC) ‡pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base C (3/4);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

CCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAG  >  NC_002947/4772977‑4773267
                                                                                                                                              |                                                                                                                                                    
ccGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCgg                                                                                                                                                <  6:308027/149‑1 (MQ=255)
       ggggCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTg                                                                                                                                         >  6:355358/1‑149 (MQ=255)
        gggCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTgg                                                                                                                                        <  1:35640/149‑1 (MQ=255)
        gggCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGGTCACCGCCAACGGCAAGTTg                                                                                                                                         <  3:220595/148‑1 (MQ=39)
         ggCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGc                                                                                                                                       >  4:298375/1‑149 (MQ=255)
          gCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCg                                                                                                                                      >  4:233081/1‑149 (MQ=255)
          gCAGGTTCAACACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCg                                                                                                                                      >  3:38818/1‑149 (MQ=255)
           cAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGa                                                                                                                                     <  5:355358/149‑1 (MQ=255)
                 cAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTTCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTGGCGCGCCATCTTGCCGTTGGCGGTGACCGGCAATTTGGCGagcagc                                                                                                                                <  4:358497/148‑1 (MQ=255)
                    cACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCAt                                                                                                                              >  4:186849/1‑147 (MQ=255)
                     aCCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAg                                                                                                                           <  4:382996/149‑1 (MQ=255)
                     aCCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAg                                                                                                                           <  3:233081/149‑1 (MQ=255)
                      ccGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAgg                                                                                                                          >  4:368515/1‑149 (MQ=255)
                       cGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGt                                                                                                                         <  7:95536/149‑1 (MQ=255)
                         tctcCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTgc                                                                                                                       <  4:370750/149‑1 (MQ=255)
                           tcCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTgcgc                                                                                                                     >  4:115818/1‑149 (MQ=255)
                             ccAGCTCGACGCCAGTTGCTGCTCCCGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGCGGTTGCGCGCCAACTTGCCGGTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCg                                                                                                                   <  4:88859/149‑1 (MQ=255)
                             ccAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCg                                                                                                                   >  1:200986/1‑149 (MQ=255)
                                 aTCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCACGTGCGCCGGCAc                                                                                                               >  4:143039/1‑149 (MQ=255)
                                    gACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCa                                                                                                             <  3:298375/148‑1 (MQ=255)
                                     aCGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGc                                                                                                                                       >  8:119672/1‑121 (MQ=255)
                                      cGCCAGTTGCTGCTCCTGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGGTCACCGCCAACGGCAAGTTGGCg                                                                                                                                      <  6:191802/121‑1 (MQ=37)
                                      cGCCAGTTGCTGCTCCCGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGGTCACCGCCAACGGCAAGTTGGCg                                                                                                                                      >  5:191802/1‑121 (MQ=37)
                                            ttGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCgg                                                                                                    <  2:257010/149‑1 (MQ=255)
                                              gctgctCCTGCGACTTGCGCTGATCGACATACACCTGCCTGGCGGCGTTGTAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCggg                                                                                                   <  7:269221/148‑1 (MQ=255)
                                                tgctCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGc                                                                                                  >  4:13150/1‑147 (MQ=255)
                                                 gctCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGGCGGCGCCATGTAATCGGGCagg                                                                                               >  4:288634/1‑147 (MQ=255)
                                                 gctCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAg                                                                                               >  3:251804/1‑149 (MQ=255)
                                                   tCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGcc                                                                                             <  6:360304/149‑1 (MQ=255)
                                                     cAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCt                                                                                           >  3:198156/1‑149 (MQ=255)
                                                     cAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCt                                                                                           <  7:26926/149‑1 (MQ=255)
                                                      aGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACATGGCGTTGGCGGTGGCCGGCAACGTGGGGAGCAGCATAACGTGCGCCGGCCCCATGTTATACGGCCAGcccg                                                                                           >  1:199485/1‑147 (MQ=255)
                                                                gcgGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGaaa                                                                                <  2:200986/149‑1 (MQ=255)
                                                                   gAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATgc                                                                             <  3:115818/149‑1 (MQ=255)
                                                                      cGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTTCCGGTGGCGGTGACCGGCAAGGTTGCCGGCAGCATCAGGTGCGCCGGCACCATGTTATAGGGCAACCCCGGCACGAGGCAAttcac                                                                           >  8:126950/1‑144 (MQ=255)
                                                                      cGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATgcgc                                                                           >  5:10636/1‑148 (MQ=255)
                                                                         cATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGt                                                                       >  4:341752/1‑149 (MQ=255)
                                                                          ctACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGtt                                                                      <  4:199436/148‑1 (MQ=255)
                                                                                ctgcctgGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCAc                                                                >  6:300388/1‑149 (MQ=255)
                                                                                        cggcgTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGt                                                                       >  1:264639/1‑134 (MQ=255)
                                                                                        cggcgTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGt                                                                       <  2:264639/134‑1 (MQ=255)
                                                                                            gTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGc                                                    >  4:338412/1‑149 (MQ=255)
                                                                                              tGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGGTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCAACATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGccc                                                  >  3:10717/1‑149 (MQ=255)
                                                                                               gAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCg                                                 >  7:216860/1‑149 (MQ=255)
                                                                                                aaaTCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCa                                                                                                             >  6:343157/1‑88 (MQ=255)
                                                                                                aaaTCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCa                                                                                                             <  5:343157/88‑1 (MQ=255)
                                                                                                  aTCCGGTGCCGGCAAGAGGTGGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGGTGGCGAGCAGCATCAGGTGCGCCGTCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCt                                              <  3:143039/149‑1 (MQ=255)
                                                                                                     cGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTg                                                                                                                        <  5:243156/72‑1 (MQ=255)
                                                                                                     cGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTg                                                                                                                        >  6:243156/1‑72 (MQ=255)
                                                                                                          ccGGCAAGAGGTTGCGCGCCAACTTGCCGTCGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGcc                                       <  8:204869/148‑1 (MQ=255)
                                                                                                           cGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCg                                     >  2:47893/1‑149 (MQ=255)
                                                                                                            ggCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCgg                                    >  6:334519/1‑149 (MQ=255)
                                                                                                               aaGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATgc                                                                             <  3:216690/105‑1 (MQ=255)
                                                                                                               aaGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATgc                                                                             >  4:216690/1‑105 (MQ=255)
                                                                                                               aaGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTc                                 <  2:199485/149‑1 (MQ=255)
                                                                                                                agagGGTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCACTGGCCCGCCTCGCTGCCCGGCTCg                                <  7:126950/149‑1 (MQ=255)
                                                                                                                         tcgcCAACTTGCCGGTCACCGCCAACGGCAAGTTGGCGCGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTgg                                                     >  2:205188/2‑119 (MQ=25)
                                                                                                                         tcgcCAACTTGCCGGTCACCGCCAACGGCAAGTTGGCGCGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTgg                                                     <  1:205188/118‑1 (MQ=25)
                                                                                                                         tcgcCAACTTGCCGGTCACCGCCAACGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGcc                             <  4:247699/142‑1 (MQ=38)
                                                                                                                         tcgcCAACTTGCCGGTCACCGCCAACGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCtagatt                       >  3:247699/2‑144 (MQ=38)
                                                                                                                           gcCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGc                     <  2:136477/149‑1 (MQ=255)
                                                                                                                            ccAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCa                    <  4:38818/149‑1 (MQ=255)
                                                                                                                               aCTTGCCGTTGGCGGTGACCGGCAAGTTGGCGCGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCa                    <  6:152677/146‑1 (MQ=255)
                                                                                                                               aCTTGCCGTTGGCGGTGACCGGCAAGTTGGCGCGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCa                    >  5:152677/1‑146 (MQ=255)
                                                                                                                                cTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAAGGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGc                  >  8:344636/1‑147 (MQ=255)
                                                                                                                                 ttGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCAt               >  6:249972/1‑149 (MQ=255)
                                                                                                                                 ttGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCAt               >  7:367526/1‑149 (MQ=255)
                                                                                                                                   gCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCAt               <  3:13150/147‑1 (MQ=255)
                                                                                                                                   gCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCAt               <  4:251804/147‑1 (MQ=255)
                                                                                                                                    ccGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCgg            <  3:288634/149‑1 (MQ=255)
                                                                                                                                     cGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCAGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTAGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGc           <  1:283729/149‑1 (MQ=255)
                                                                                                                                       ttGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGcac         >  8:358980/1‑149 (MQ=255)
                                                                                                                                              tGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAg  >  7:336464/1‑149 (MQ=255)
                                                                                                                                              tGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAg  <  7:358980/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAG  >  NC_002947/4772977‑4773267

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