Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,120,859 C→G 8.7% A431P (GCC→CCC) ‡ PP_5371 ← rubredoxin/rubredoxin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,120,8590CG 8.7% 217.4 / 9.6 69A431P (GCC→CCC) ‡PP_5371rubredoxin/rubredoxin reductase
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base G (2/4);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

CGGTCATCATGAGGCAGATTGACGATACTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCAATGGGTAGTGGGCGGCG  >  NC_002947/6120720‑6121007
                                                                                                                                           |                                                                                                                                                    
cGGTCATCATGAGGCAGATTGACGATACTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGTGCCCTGCgcg                                                                                                                                             <  6:228000/149‑3 (MQ=255)
cGGTCATCATGAGGCAGATTGACGATACTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTg                                                                                                                                             <  5:50234/149‑1 (MQ=255)
 ggTCATCATGAGGCAGATTGACGATACTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGTGCCCTGCGCgg                                                                                                                                            <  3:234332/149‑1 (MQ=255)
                       gATACTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATc                                                                                                                      <  8:111454/149‑1 (MQ=255)
                         tACTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCa                                                                                                                    <  8:242652/149‑1 (MQ=255)
                          aCTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCa                                                                                                                    >  3:182624/1‑148 (MQ=255)
                          aCTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTCTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCAt                                                                                                                   >  8:176736/1‑149 (MQ=255)
                           cTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCAt                                                                                                                   <  4:182624/148‑1 (MQ=255)
                            tGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGt                                                                                                                                           <  6:125859/123‑1 (MQ=255)
                            tGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGt                                                                                                                                           >  5:125859/1‑123 (MQ=255)
                            tGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGAAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCg                                                                                                                 >  6:286723/1‑149 (MQ=255)
                               tCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCa                                                                                                                                                     <  6:106009/110‑1 (MQ=255)
                               tCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCa                                                                                                                                                     >  5:106009/1‑110 (MQ=255)
                                 aTGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGt                                                                                                            >  5:278416/1‑149 (MQ=255)
                                  tGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGtc                                                                                                           >  3:105366/1‑149 (MQ=255)
                                     ttCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGtcgca                                                                                                        >  1:168553/1‑149 (MQ=255)
                                      tCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGTGCCCTGCGCgg                                                                                                                                            <  2:136835/112‑1 (MQ=255)
                                      tCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGTGCCCTGCGCgg                                                                                                                                            >  1:136835/1‑112 (MQ=255)
                                       cAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGtcgcat                                                                                                       <  2:168553/148‑1 (MQ=255)
                                          ccTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGtcgcatcgca                                                                                                   <  4:73940/149‑1 (MQ=255)
                                           cTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCACGCAGGTGCCCTGCgc                                                                                                                                              >  8:333949/1‑104 (MQ=255)
                                           cTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCACGCAGGTGCCCTGCgc                                                                                                                                              <  7:333949/105‑2 (MQ=255)
                                                  cacaTTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCAtt                                                                                           >  8:167320/1‑149 (MQ=255)
                                                  cacaTTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCAtt                                                                                           <  2:90597/149‑1 (MQ=255)
                                                   acaTTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTg                                                                                          <  7:167320/149‑1 (MQ=255)
                                                      ttttGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGcc                                                                                       >  2:31261/1‑149 (MQ=255)
                                                                    cTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAgc                                                                                                                              <  4:215666/96‑1 (MQ=255)
                                                                    cTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAgc                                                                                                                              >  3:215666/1‑96 (MQ=255)
                                                                    cTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAAt                                                                                                                       <  3:39946/103‑1 (MQ=255)
                                                                    cTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAAt                                                                                                                       >  4:39946/1‑103 (MQ=255)
                                                                      ggCAGTTTTGCAGCGAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGa                                                                       <  2:157138/149‑1 (MQ=255)
                                                                       gCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAg                                                                      <  3:145729/149‑1 (MQ=255)
                                                                         aGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGgt                                                                    <  7:172681/149‑1 (MQ=255)
                                                                           ttttGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGgtgt                                                                  <  5:286723/149‑1 (MQ=255)
                                                                              tGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTc                                                                 >  2:351771/1‑147 (MQ=255)
                                                                              tGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTc                                                                 <  1:351771/147‑1 (MQ=255)
                                                                                 aGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTcagcag                                                            <  1:31261/149‑1 (MQ=255)
                                                                                       gTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCTCCGGGAGGTGTcagcag                                                            >  4:186641/1‑143 (MQ=255)
                                                                                       gTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCTCCGGGAGGTGTcagcag                                                            <  3:186641/143‑1 (MQ=255)
                                                                                             cGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGga                                                 <  6:342206/148‑1 (MQ=255)
                                                                                                   aTTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAg                                                                      <  4:300702/121‑1 (MQ=255)
                                                                                                   aTTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAg                                                                      >  3:300702/1‑121 (MQ=255)
                                                                                                   aTTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCgg                                          >  7:251419/1‑149 (MQ=255)
                                                                                                    gttGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCggg                                         <  6:278416/148‑1 (MQ=255)
                                                                                                          ccTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGtttt                                    <  8:251419/148‑1 (MQ=255)
                                                                                                           cTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTgg                                                                                                             <  3:207975/74‑1 (MQ=255)
                                                                                                           cTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTgg                                                                                                             >  4:207975/1‑74 (MQ=255)
                                                                                                            agACGGCCGAGCAGCACGAACCCGCGCACCGCACCAGCGTGCTCATGAAACGCACCGCTCAAACCATCGCTCGTCGCATCGCAGCGCCAATGGCCCTCGCAACCACCGCGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCAc                                 >  8:130909/2‑149 (MQ=255)
                                                                                                                ggCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTc                             >  7:29805/1‑149 (MQ=255)
                                                                                                                  ccGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCAc                           >  3:340587/1‑149 (MQ=255)
                                                                                                                   cGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACg                          >  5:25363/1‑149 (MQ=255)
                                                                                                                         gcaCGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGgcgc                     >  3:223735/1‑148 (MQ=255)
                                                                                                                          caCGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGcc                    >  6:30890/1‑148 (MQ=255)
                                                                                                                          caCGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGcc                    >  5:22760/1‑148 (MQ=255)
                                                                                                                          caCGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGcc                    >  4:272546/1‑148 (MQ=255)
                                                                                                                           aCGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCaa                  <  4:105366/149‑1 (MQ=255)
                                                                                                                             gAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCAATg                >  3:205642/1‑149 (MQ=255)
                                                                                                                              aaCCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCAATgg               <  6:22760/149‑1 (MQ=255)
                                                                                                                              aaCCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCAATgg               <  4:205642/149‑1 (MQ=255)
                                                                                                                                 ccGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCAATggg              >  7:236396/1‑147 (MQ=255)
                                                                                                                                 ccGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCAATggg              >  4:235804/1‑147 (MQ=255)
                                                                                                                                 ccGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCAATGGGTa            <  4:223735/149‑1 (MQ=255)
                                                                                                                                    cgcAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCg                                                                          >  7:19295/1‑84 (MQ=255)
                                                                                                                                    cgcAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCg                                                                          <  8:19295/84‑1 (MQ=255)
                                                                                                                                    cgcAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCggg                                         >  2:108453/1‑117 (MQ=255)
                                                                                                                                    cgcAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCggg                                         <  1:108453/117‑1 (MQ=255)
                                                                                                                                     gcAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCAATGGGTAGTgg        <  3:272546/149‑1 (MQ=255)
                                                                                                                                       aGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGgcgc                     >  7:18359/1‑134 (MQ=255)
                                                                                                                                       aGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGgcgc                     <  8:18359/134‑1 (MQ=255)
                                                                                                                                           cACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCAATGGGTAGTGGgcggcg  <  8:29805/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CGGTCATCATGAGGCAGATTGACGATACTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCAATGGGTAGTGGGCGGCG  >  NC_002947/6120720‑6121007

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: