Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,773,114 G→C 9.5% A2051A (GCC→GCG) ‡ pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,773,1140GC 9.5% 216.2 / 12.3 74A2051A (GCC→GCG) ‡pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base C (3/4);  total (37/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 = 3.83e-01

TCCTGCATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCAC  >  NC_002947/4772968‑4773260
                                                                                                                                                  |                                                                                                                                                  
tCCTGCATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTgg                                                                                                                                                   <  7:50232/148‑1 (MQ=255)
tCCTGCATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGc                                                                                                                                                  <  4:40999/149‑1 (MQ=255)
  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                                                                                                                    <  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                                                                <  2:264639/134‑1 (MQ=255)
                                                                                                 cggcgTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGt                                                                >  1:264639/1‑134 (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                                                                                                                 >  6:243156/1‑72 (MQ=255)
                                                                                                              cGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTg                                                                                                                 <  5:243156/72‑1 (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             >  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)
                                                                                                                                                  |                                                                                                                                                  
TCCTGCATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCAC  >  NC_002947/4772968‑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: