Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,773,117 G→C 9.7% T2050T (ACC→ACG) ‡ pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,773,1170GC 9.7% 215.8 / 11.0 72T2050T (ACC→ACG) ‡pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (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 = 9.28e-01

CTGCATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCAC  >  NC_002947/4772970‑4773260
                                                                                                                                                   |                                                                                                                                               
cTGCATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCgg                                                                                                                                                >  1:257010/1‑149 (MQ=255)
       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                                                                                                                    <  3:233081/149‑1 (MQ=255)
                            aCCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAg                                                                                                                    <  4:382996/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                                                                                                      <  5:343157/88‑1 (MQ=255)
                                                                                                       aaaTCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCa                                                                                                      >  6:343157/1‑88 (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                                                                      >  4:216690/1‑105 (MQ=255)
                                                                                                                      aaGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATgc                                                                      <  3:216690/105‑1 (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             >  5:152677/1‑146 (MQ=255)
                                                                                                                                      aCTTGCCGTTGGCGGTGACCGGCAAGTTGGCGCGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCa             <  6:152677/146‑1 (MQ=255)
                                                                                                                                       cTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAAGGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGc           >  8:344636/1‑147 (MQ=255)
                                                                                                                                        ttGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCAt        >  7:367526/1‑149 (MQ=255)
                                                                                                                                        ttGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCAt        >  6:249972/1‑149 (MQ=255)
                                                                                                                                          gCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCAt        <  4:251804/147‑1 (MQ=255)
                                                                                                                                          gCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCAt        <  3:13150/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)
                                                                                                                                                   |                                                                                                                                               
CTGCATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCAC  >  NC_002947/4772970‑4773260

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: