Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,778,460 C→A 9.5% E269D (GAG→GAT)  pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,778,4600CA 9.5% 225.9 / 11.6 74E269D (GAG→GAT) pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (39/28);  new base A (5/2);  total (44/30)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

GTTGGCGATCGGTACGCCCACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCTTCCATCCAGTTG  >  NC_002947/4778313‑4778604
                                                                                                                                                   |                                                                                                                                                
gTTGGCGATCGGTACGCCCACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGAtct                                                                                                                                                 >  4:66860/1‑149 (MQ=255)
    gCGATCGGTACGCCCACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGAt                                                                                                                                             >  5:62374/1‑149 (MQ=255)
      gATCGGTACGCCCACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGc                                                                                                                                           >  4:263797/1‑149 (MQ=255)
       aTCGGTACGCCCACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGAtct                                                                                                                                                 >  2:334481/1‑142 (MQ=255)
       aTCGGTACGCCCACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGAtct                                                                                                                                                 <  1:334481/142‑1 (MQ=255)
           gTACGCCCACGCGCACCTGGGGCTGGCCGCTGTAGCCGTGCAGCAAGGACACGAACGATGCCAGCAACACCATCAACAACCTTACACCCTCGACAATCGCCAACCCCTTCAGCGACAGAAGAAACAAGGCAACCAACAAAATGCCCAcc                                                                                                                                      >  6:280980/1‑149 (MQ=255)
            tACGCCCACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCAccc                                                                                                                                     >  5:28551/1‑149 (MQ=255)
               gCCCACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCcgcg                                                                                                                                  >  6:332829/1‑149 (MQ=255)
               gCCCACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGAGATCGATGccgtgtt                                                                                                                                     <  4:39068/146‑6 (MQ=255)
               gCCCACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGAGATCGATGccgtgtt                                                                                                                                     >  3:39068/1‑141 (MQ=255)
                cccACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCcgcgc                                                                                                                                 <  5:332829/149‑1 (MQ=255)
                 ccACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCcgcgcg                                                                                                                                >  3:143801/1‑149 (MQ=255)
                 ccACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTCACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCcgcgcg                                                                                                                                >  3:167560/1‑149 (MQ=255)
                  cACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACCCGGCATCGATCTCGATGCCCACCCgcgcgc                                                                                                                               <  1:205532/149‑1 (MQ=255)
                  cACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCgcgcgc                                                                                                                               <  6:42862/149‑1 (MQ=255)
                          ccTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCg                                                                                                                       <  6:28551/149‑1 (MQ=255)
                                 cTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGgcg                                                                                                                >  6:212074/1‑149 (MQ=255)
                                   ggCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGcc                                                                                                              >  6:56384/1‑149 (MQ=255)
                                      cGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGAtctc                                                                                                                                                >  2:163918/1‑112 (MQ=255)
                                      cGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGAtctc                                                                                                                                                <  1:163918/112‑1 (MQ=255)
                                      cGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCAcc                                                                                                           >  7:37412/1‑149 (MQ=255)
                                                     gcaAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGcaaccg                                                                                            >  7:225587/1‑144 (MQ=255)
                                                     gcaAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGAt                                                                                              <  2:173223/147‑1 (MQ=255)
                                                     gcaAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGAt                                                                                              >  1:173223/1‑147 (MQ=255)
                                                     gcaAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTCACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTg                                                                                            <  4:167560/149‑1 (MQ=255)
                                                      caAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTgg                                                                                           >  1:33838/1‑149 (MQ=255)
                                                       aaGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGt                                                                                          >  4:297908/1‑149 (MQ=255)
                                                        aGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTc                                                                                         <  5:56384/149‑1 (MQ=255)
                                                          gTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCgg                                                                                       <  6:66950/149‑1 (MQ=255)
                                                           tCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGc                                                                                      <  5:280980/149‑1 (MQ=255)
                                                            cTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCgg                                                                                       >  8:137435/1‑147 (MQ=255)
                                                               gAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGc                                                                                  <  6:62374/149‑1 (MQ=255)
                                                               gAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATAGGTCGGCAGGc                                                                                  >  1:173555/1‑149 (MQ=255)
                                                                   gATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAg                                                                              >  2:15536/1‑149 (MQ=255)
                                                                    aTGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCAccggccgg                                                                                                     <  3:298100/125‑1 (MQ=255)
                                                                    aTGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCAccggccgg                                                                                                     >  4:298100/1‑125 (MQ=255)
                                                                       ccAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCAc                                                                           <  3:66860/148‑1 (MQ=255)
                                                                           cAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGtt                                                                      >  1:141450/1‑149 (MQ=255)
                                                                                gCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCAccggcc                                                                                                       >  6:337615/1‑111 (MQ=255)
                                                                                gCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCAccggcc                                                                                                       <  5:337615/111‑1 (MQ=255)
                                                                                  aTGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCgg                                                                                       >  1:143487/1‑125 (MQ=255)
                                                                                  aTGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCgg                                                                                       <  2:143487/125‑1 (MQ=255)
                                                                                   tGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCccg                                                              <  4:143801/149‑1 (MQ=255)
                                                                                   tGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCccg                                                              >  1:299180/1‑149 (MQ=255)
                                                                                         aGGGTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAgg                                                        <  7:137435/149‑1 (MQ=255)
                                                                                          ggTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAggg                                                       >  5:342521/1‑149 (MQ=255)
                                                                                                acCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCt                                                 >  2:257354/1‑149 (MQ=255)
                                                                                                  ccTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTc                                               >  3:27098/1‑149 (MQ=255)
                                                                                                   cTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTcc                                              <  1:257354/149‑1 (MQ=255)
                                                                                                    tGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCt                                                                   >  4:258207/1‑127 (MQ=255)
                                                                                                    tGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCt                                                                   <  3:258207/127‑1 (MQ=255)
                                                                                                     gggCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTc                                               <  5:141372/146‑1 (MQ=255)
                                                                                                     gggCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTc                                               >  6:141372/1‑146 (MQ=255)
                                                                                                      ggCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAgg                                                        <  2:99764/136‑1 (MQ=255)
                                                                                                      ggCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAgg                                                        >  1:99764/1‑136 (MQ=255)
                                                                                                      ggCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGt                                           <  3:297908/149‑1 (MQ=255)
                                                                                                        cATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTaa                                         <  4:27098/149‑1 (MQ=255)
                                                                                                          ttGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGc                                       <  3:263797/149‑1 (MQ=255)
                                                                                                              ccAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGc                                   >  2:260626/1‑149 (MQ=255)
                                                                                                              ccAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGc                                   <  6:342521/149‑1 (MQ=255)
                                                                                                               cAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCt                                  >  6:226071/1‑149 (MQ=255)
                                                                                                                aaCGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTACGCCAGCTg                                 >  8:63259/1‑149 (MQ=255)
                                                                                                                   gCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTTCAGTAAGCCAGCTGCcg                              >  7:1797/1‑149 (MQ=255)
                                                                                                                     cTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCcgcg                            >  7:271906/1‑149 (MQ=255)
                                                                                                                            gACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTc                     <  5:212074/149‑1 (MQ=255)
                                                                                                                                   gcaACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCt              >  8:236944/1‑149 (MQ=255)
                                                                                                                                   gcaACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCt              <  2:299180/149‑1 (MQ=255)
                                                                                                                                    caACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCtt             >  7:74614/1‑149 (MQ=255)
                                                                                                                                    caACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCtt             >  8:159970/1‑149 (MQ=255)
                                                                                                                                       gtgggCATCGAGATCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCt              <  6:96450/142‑1 (MQ=255)
                                                                                                                                       gtgggCATCGAGATCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCt              >  5:96450/4‑145 (MQ=255)
                                                                                                                                       gtgggCATCGAGATCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCTtcc           >  4:262060/4‑148 (MQ=255)
                                                                                                                                       cacGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCTtcca          >  7:158307/1‑149 (MQ=255)
                                                                                                                                        tgggCATCGAGATCGATGCCGACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCTtccat         >  8:237136/3‑149 (MQ=255)
                                                                                                                                        tgggCATCGAGATCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCTtccat         >  3:298713/3‑149 (MQ=255)
                                                                                                                                        acGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCT‑CCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCTtccatc        >  4:299359/1‑149 (MQ=255)
                                                                                                                                         cGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCTCTtccatc        <  8:37412/149‑1 (MQ=255)
                                                                                                                                            cATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCTTCCATCCAg     <  7:159970/149‑1 (MQ=255)
                                                                                                                                               cGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCTTCCATCCAGTTg  >  2:73462/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTTGGCGATCGGTACGCCCACGCGCACCTGGGGCTGGCCGCTGTAGCGGTGCAGCAAGGTCTGGAACGATGCCAGCAACAGCATGAACAAGGTTACACCCTGGGCATTGGCCAACGCCTTGAGCGACTGCAGCAACACGGCATCGATCTCGATGCCCACCCGCGCGCCGCGCCCGCTGGGCGCCACCGGCCGGGAGCGATTGGTCGGCAGGCTGAGCACCGGTTGCTCCCCGCCGAGGGTGTCCTTCCAGTAAGCCAGCTGCCGCGCCTGCTCACCCGCTTCCATCCAGTTG  >  NC_002947/4778313‑4778604

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