Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,120,858 G→T 8.7% A431D (GCC→GAC) ‡ PP_5371 ← rubredoxin/rubredoxin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,120,8580GT 8.7% 212.0 / 8.3 69A431D (GCC→GAC) ‡PP_5371rubredoxin/rubredoxin reductase
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base T (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 = 9.57e-01

CCGCAATTGCGGTCATCATGAGGCAGATTGACGATACTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCAATGGGTAGTGG  >  NC_002947/6120711‑6121001
                                                                                                                                                   |                                                                                                                                               
ccGCAATTGCGGTCATCATGAGGCAGATTGACGATACTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAggg                                                                                                                                                 >  6:50234/1‑148 (MQ=255)
         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                                                                                                                                        <  7:333949/105‑2 (MQ=255)
                                                    cTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCACGCAGGTGCCCTGCgc                                                                                                                                        >  8:333949/1‑104 (MQ=255)
                                                           cacaTTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCAtt                                                                                     <  2:90597/149‑1 (MQ=255)
                                                           cacaTTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCAtt                                                                                     >  8:167320/1‑149 (MQ=255)
                                                            acaTTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTg                                                                                    <  7:167320/149‑1 (MQ=255)
                                                               ttttGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGcc                                                                                 >  2:31261/1‑149 (MQ=255)
                                                                             cTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAgc                                                                                                                        >  3:215666/1‑96 (MQ=255)
                                                                             cTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAgc                                                                                                                        <  4:215666/96‑1 (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                                                           <  1:351771/147‑1 (MQ=255)
                                                                                       tGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTc                                                           >  2:351771/1‑147 (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              >  5:22760/1‑148 (MQ=255)
                                                                                                                                   caCGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGcc              >  6:30890/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)
                                                                                                                                                   |                                                                                                                                               
CCGCAATTGCGGTCATCATGAGGCAGATTGACGATACTGTTCATGGTTCAGCCTCAGGCCACATTTTGTTGGGTGTCCTGGCAGTTTTGCAGCCAGGTATTTCGTTGCATTTGTGCCTGACGGCCGAGCAGCACGAACCCGCGCAGGGCACCTGCGTGGTCATGAAACGCAGCGCTCAATCCATCGCTGGTCGCATCGCAGCGCCATTGGCCCTCGCAACCACCGGGAGGTGTCAGCAGGCACAGAGGGACGACCGGGGTTTTCACCGTCACGGGCGCCAATGGGTAGTGG  >  NC_002947/6120711‑6121001

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: