Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,773,136 A→C 5.6% L2044R (CTC→CGC)  pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,773,1360AC 5.6% 219.2 / 3.1 74L2044R (CTC→CGC) pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base C (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

CAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAGCCCACCAGTTGCTTG  >  NC_002947/4772988‑4773282
                                                                                                                                                    |                                                                                                                                                  
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)
                                 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                                                                                                          <  7:26926/149‑1 (MQ=255)
                                          cAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCt                                                                                                          >  3:198156/1‑149 (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                                                                                            >  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                                                                    <  1:205188/118‑1 (MQ=25)
                                                                                                              tcgcCAACTTGCCGGTCACCGCCAACGGCAAGTTGGCGCGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTgg                                                                    >  2:205188/2‑119 (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                              >  6:249972/1‑149 (MQ=255)
                                                                                                                      ttGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCAt                              >  7:367526/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)
                                                                                                                                   tGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAg                 >  7:336464/1‑149 (MQ=255)
                                                                                                                                   tGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAg                 <  7:358980/149‑1 (MQ=255)
                                                                                                                                    gACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAGc                >  3:360470/1‑149 (MQ=255)
                                                                                                                                     aCCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAGc                <  3:368515/148‑1 (MQ=255)
                                                                                                                                      ccGGGAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAGcc               <  4:10717/148‑1 (MQ=255)
                                                                                                                                       cGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAGCcca             >  2:61942/1‑149 (MQ=255)
                                                                                                                                        ggCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAGCccac            <  4:198156/149‑1 (MQ=255)
                                                                                                                                          cAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAGCccacca          <  6:10636/149‑1 (MQ=255)
                                                                                                                                            aGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCgg                           >  6:68860/1‑130 (MQ=255)
                                                                                                                                            aGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCgg                           <  5:68860/130‑1 (MQ=255)
                                                                                                                                             gTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCg                                                                                                     <  2:787/55‑1 (MQ=255)
                                                                                                                                              ttGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAGCCCACCAGttg      >  5:357083/1‑149 (MQ=255)
                                                                                                                                                 gCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAGCCCACCAGttgct    <  5:300388/148‑1 (MQ=255)
                                                                                                                                                  cGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAGCCCACCAGttgcttg  <  5:334519/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACCACATAGCCCACCAGTTGCTTG  >  NC_002947/4772988‑4773282

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: ‑