Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,781,591 G→C 5.8% P1908R (CCC→CGC)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,781,5910GC 5.8% 259.3 / 4.3 86P1908R (CCC→CGC) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (42/39);  new base C (2/3);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

GAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGTGACTTGAGCCACGCCTTGGCCGT  >  NC_002947/4781446‑4781739
                                                                                                                                                 |                                                                                                                                                      
gAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGgcg                                                                                                                                                     >  1:241990/1‑149 (MQ=255)
gAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGgcg                                                                                                                                                     >  4:265611/1‑149 (MQ=255)
  aCAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCTAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGcgcg                                                                                                                                                   <  5:203458/149‑1 (MQ=255)
    aGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGtt                                                                                                                                                 >  6:129289/1‑149 (MQ=255)
            gTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaac                                                                                                                                          <  4:124717/148‑1 (MQ=255)
              cATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaacacg                                                                                                                                       <  8:97057/149‑1 (MQ=255)
               aTACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaacacg                                                                                                                                       <  2:32318/148‑1 (MQ=255)
                tACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaacacgaa                                                                                                                                     <  6:93613/149‑1 (MQ=255)
                   gACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAg                                                                                                                                  <  8:199161/149‑1 (MQ=255)
                    aCAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAgg                                                                                                                                 >  2:131448/1‑149 (MQ=255)
                           cAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCg                                                                                                                           >  8:192575/1‑148 (MQ=255)
                            aGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑aa                                                                                                                       >  5:38051/1‑149 (MQ=255)
                                ccAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaa                                                                                                                                           <  6:162984/127‑1 (MQ=255)
                                ccAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaa                                                                                                                                           >  5:162984/1‑127 (MQ=255)
                                 cAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAggg                                                                                                                  <  2:241990/149‑1 (MQ=255)
                                      ccTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGa                                                                                                             >  1:713/1‑149 (MQ=255)
                                       cTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGCGCGCGTtctc                                                                                                                                              >  2:188895/1‑117 (MQ=255)
                                       cTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGCGCGCGTtctc                                                                                                                                              >  2:183017/1‑117 (MQ=255)
                                       cTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGCGCGCGTtctc                                                                                                                                              <  1:183017/117‑1 (MQ=255)
                                       cTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGCGCGCGTtctc                                                                                                                                              <  1:188895/117‑1 (MQ=255)
                                          aCCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGAcc                                                                                                           <  4:153585/147‑1 (MQ=255)
                                          aCCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGCGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGc                                                                                                         <  1:270971/149‑1 (MQ=255)
                                              gggTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGAATCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTAcc                                                                                                     >  2:316339/1‑149 (MQ=255)
                                                   caccaccacGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCt                                                                                                <  6:132775/149‑1 (MQ=255)
                                                      caccacGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGc                                                                                             >  7:322086/1‑149 (MQ=255)
                                                      aaccacGGTGAGCGTGTAGTTGGTGTGCGTGCGCGAACTCACTGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGc                                                                                             <  8:334310/148‑1 (MQ=255)
                                                         cacGGTGAGCGGGTAGTTGGTGTGCGTGCGCTAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTcgaacacgaac                                                                                                                                    <  7:156195/109‑1 (MQ=255)
                                                         cacGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCg                                                                                                    <  3:38350/139‑1 (MQ=255)
                                                         cacGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCg                                                                                                    >  4:38350/1‑139 (MQ=255)
                                                         cacGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCAAGGGCGCGTTCTccaacacgaaa                                                                                                                                    >  8:156195/1‑108 (MQ=255)
                                                                aGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTg                                                                                   <  7:192575/149‑1 (MQ=255)
                                                                     gTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGc                                                                                             <  5:204388/134‑1 (MQ=255)
                                                                     gTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGc                                                                                             >  6:204388/1‑134 (MQ=255)
                                                                                   tgcgCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAg                                                                <  4:14908/149‑1 (MQ=255)
                                                                                       cgAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAgtgct                                                            >  5:65659/1‑149 (MQ=255)
                                                                                        gAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATccac                                                                          >  4:194782/1‑134 (MQ=255)
                                                                                        gAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATccac                                                                          <  3:194782/134‑1 (MQ=255)
                                                                                        gAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAgtgctc                                                           >  5:159975/1‑149 (MQ=255)
                                                                                            tCCCATCCTTGATCTTCATGGTCGAAGCCTGGGCGCGAACGCTGGACTCCACGGGCGCTTTCTCGAACACGAACAGGCTGTAG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCAc                                                                                        <  4:49865/116‑1 (MQ=39)
                                                                                            tCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCAc                                                                                        >  3:49865/1‑116 (MQ=255)
                                                                                            tCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCAc                                                                                        >  3:124504/1‑116 (MQ=255)
                                                                                            tCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCAc                                                                                        <  4:124504/116‑1 (MQ=255)
                                                                                            tCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAgtgctcgtgc                                                       >  4:355028/1‑149 (MQ=255)
                                                                                              ccAGCCTTGAGCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAgtgctcgtgc                                                       <  1:131448/147‑1 (MQ=255)
                                                                                              ccAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTc                                                                                                                            <  6:66046/80‑1 (MQ=255)
                                                                                              ccAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTc                                                                                                                            >  5:66046/1‑80 (MQ=255)
                                                                                              ccAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGc                                                                                                                 <  6:92522/89‑1 (MQ=255)
                                                                                              ccAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCt                                                                                    >  3:184046/1‑118 (MQ=255)
                                                                                              ccAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCt                                                                                    <  4:184046/118‑1 (MQ=255)
                                                                                              ccAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAgtgctcgtgctc                                                     >  1:296312/1‑149 (MQ=255)
                                                                                              ccAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAgtgctcgtgctc                                                     >  2:371409/1‑149 (MQ=255)
                                                                                              ccAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAgtgctcgtgctc                                                     >  4:311855/1‑149 (MQ=255)
                                                                                              ccAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAgtgctcgtgctc                                                     >  6:234353/1‑149 (MQ=255)
                                                                                              ccAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAgtgctcgtgctc                                                     >  5:306687/1‑149 (MQ=255)
                                                                                              ccAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAgtgctcgtgctc                                                     >  3:45553/1‑149 (MQ=255)
                                                                                                      gACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCagg                                              <  3:311855/148‑1 (MQ=255)
                                                                                                       aCCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCaggtt                                            >  6:215863/1‑149 (MQ=255)
                                                                                                       aCCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCaggtt                                            <  6:159975/149‑1 (MQ=255)
                                                                                                         cTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAgtgc                                                             >  8:186422/1‑130 (MQ=255)
                                                                                                         cTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAgtgc                                                             <  7:186422/130‑1 (MQ=255)
                                                                                                         cTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCaggtt                                            >  1:355855/1‑147 (MQ=255)
                                                                                                          ttCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCaggttgag                                         >  5:126361/1‑149 (MQ=255)
                                                                                                           tCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCaggttgagg                                        >  5:328727/1‑149 (MQ=255)
                                                                                                            cATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCaggttgaggt                                       >  5:16273/1‑149 (MQ=255)
                                                                                                              tGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCaggttgaggttg                                     <  3:265611/149‑1 (MQ=255)
                                                                                                                gTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTtgtg                                   <  3:355028/149‑1 (MQ=255)
                                                                                                                gTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTtgt                                    <  5:129289/148‑1 (MQ=255)
                                                                                                                  cGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGCCCGCTACCGACCTCGCTGCACGCCTGGATATCCCCCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGt                                 >  2:304746/1‑149 (MQ=255)
                                                                                                                  cGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGt                                 <  4:45553/149‑1 (MQ=255)
                                                                                                                     aGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCaa                              <  5:234353/149‑1 (MQ=255)
                                                                                                                         tGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGTCCTGACCGCTACCGACCTCGCTGCACGCATGGATATCCACCAAAGGCAACTGCTCGTGCTCCCGCACGTTGAAGTTGTGGTCaaa                             >  5:292158/1‑145 (MQ=255)
                                                                                                                           ggCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCaaa                                                                      <  6:114637/103‑1 (MQ=255)
                                                                                                                           ggCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCaaa                                                                      >  5:114637/1‑103 (MQ=255)
                                                                                                                               ccAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGTGACtt                    <  6:16273/149‑1 (MQ=255)
                                                                                                                                  aCGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGTGACTTGAg                 <  1:371409/149‑1 (MQ=255)
                                                                                                                                   cGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGTGACTTGAGc                >  3:309374/1‑149 (MQ=255)
                                                                                                                                    gCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGt                                 <  7:17971/131‑1 (MQ=255)
                                                                                                                                    gCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGt                                 >  8:17971/1‑131 (MQ=255)
                                                                                                                                        gACTCCACGGGCGCGTTCTCGAACAAGAACAGGATGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGTGACTTGAGCCAc             <  6:292158/147‑1 (MQ=255)
                                                                                                                                             cACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCaa                              <  7:332898/125‑1 (MQ=255)
                                                                                                                                             cACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCaa                              >  8:332898/1‑125 (MQ=255)
                                                                                                                                             cACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGTGACTTGAGCCACGCCTTgg      >  1:163790/1‑149 (MQ=255)
                                                                                                                                             cACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGTGACTTGAGCCACGCCTTgg      >  5:354998/1‑149 (MQ=255)
                                                                                                                                              aCGGGCGCGTTCTCGAACAAGAACATGCTGTCGAAAACAGGGTCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGTGACTTGAGCCACGCCTTg       <  1:316339/149‑1 (MQ=255)
                                                                                                                                               cGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGt                         <  6:126744/128‑1 (MQ=255)
                                                                                                                                               cGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGt                         >  5:126744/1‑128 (MQ=255)
                                                                                                                                                 ggCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGTGACTTGAGCCACGCCTTGGCCGt  <  8:322086/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
GAACAGGCGTTGGTCATACGACAGGTGCAGGCCCAGCACCTCACCCGGGTACACCACCACGGTGAGCGGGTAGTTGGTGTGCGTGCGCGAACTCCCAGCCTTGACCTTCATGGTCGAAGCCTGGGCGCCAACGCTGGACTCCACGGGCGCGTTCTCGAACACGAACAGGCTGTCG‑‑AACAGGGCCTGACCGCTACCGACCTCGCTGCACGCCTGGATATCCACCAAAGGCAAGTGCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGTGACTTGAGCCACGCCTTGGCCGT  >  NC_002947/4781446‑4781739

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: