Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 857,639 G→C 5.3% D395E (GAC→GAG) ‡ phr ← deoxyribodipyrimidine photolyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947857,6390GC 5.3% 249.4 / 3.0 75D395E (GAC→GAG) ‡phrdeoxyribodipyrimidine photolyase
Reads supporting (aligned to +/- strand):  ref base G (31/40);  new base C (2/2);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

AAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGC  >  NC_002947/857493‑857781
                                                                                                                                                   |                                                                                                                                                  
aaaaGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTACTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTCAAACCGGGTTGAAGAAACGGAAATAGGGCACagcgac                                                                                                                                                  >  2:153791/1‑147 (MQ=255)
  aaGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCgg                                                                                                                                                >  5:165311/1‑149 (MQ=255)
   agttcagTCGATTTAATTGGTTGATGGATGGCTTTTTCATCTAGCTCCTGTAGTTCCGGTACCCAGTTGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGTCACCGCGTCGGt                                                                                                                                               <  3:308885/143‑1 (MQ=255)
   aGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGt                                                                                                                                               <  7:45352/149‑1 (MQ=255)
        cgtctttttATTTGGTTGATGGATGGCTTTTTCATGGATCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCgg                                                                                                                                          <  8:331664/143‑1 (MQ=255)
          tCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTg                                                                                                                                        <  3:67270/149‑1 (MQ=255)
            gATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTgg                                                                                                                                       <  4:322930/148‑1 (MQ=255)
             aTTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGaa                                                                                                                                     <  3:149136/149‑1 (MQ=255)
             aTTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGaa                                                                                                                                     <  5:405991/149‑1 (MQ=255)
               ttAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGc                                                                                                                                   >  2:254654/1‑149 (MQ=255)
                 aaTGGGTTGATGGATTGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCa                                                                                                                                  <  5:262579/148‑1 (MQ=255)
                     ggctGATGCATGGCTGTATCATCGAGCGCCGCGAGGACCGGTACCCAGTGGCGAAAGCACCGACCCTGTGGATCGAAGCGCGGCGACTGTGAAACCGGGGTGAAGATACGGAAATAGCGCCCCGCGTCGGTGCCGGTGGAAGCACTCCa                                                                                                                            <  4:262195/146‑1 (MQ=255)
                       ttGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACt                                                                                                                          <  5:216936/149‑1 (MQ=255)
                       ttGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACt                                                                                                                          >  5:157970/1‑149 (MQ=255)
                          atggatggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGc                                                                                                                        <  6:157970/148‑1 (MQ=255)
                              atggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGc                                                                                                                                             <  7:175961/124‑1 (MQ=255)
                              atggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTgg                                                                                                                                       >  8:237259/1‑130 (MQ=255)
                              atggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTgg                                                                                                                                       <  7:237259/130‑1 (MQ=255)
                              atggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGATCAAGATACGGAAATAGGGCACCGCGTCgcagc                                                                                                                                             >  8:175961/1‑120 (MQ=255)
                                 gCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAAccccc                                                                                                                >  2:134375/1‑149 (MQ=255)
                                  cTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCg                                                                                                               <  7:296098/149‑1 (MQ=255)
                                   tttttCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCgt                                                                                                              >  3:392739/1‑149 (MQ=255)
                                     tttCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCgtt                                                                                                             <  4:147038/148‑1 (MQ=255)
                                        cATCGAGCTCCTGTAGTTCCGGTACCCAGTTGCGAATGAACCGTCCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCgttgtt                                                                                                          <  8:340564/148‑1 (MQ=255)
                                         aTCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGttgttg                                                                                                         <  8:86884/148‑1 (MQ=255)
                                            gAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTggcgg                                                                                                     <  2:94892/149‑1 (MQ=255)
                                                  cTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGt                                                                                               >  8:14212/1‑149 (MQ=255)
                                                  cTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGACGCGGTGccctat                                                                                                                                        <  7:129936/109‑5 (MQ=255)
                                                  cTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGACGCGGTGccctat                                                                                                                                        >  8:129936/1‑105 (MQ=255)
                                                   tGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGACCCGCGTCGGTGCCGGTGGAAGCCCTCCACTGCCAACCCCCGTTGTTGGCGGCCCGGTc                                                                                              >  5:131966/1‑149 (MQ=255)
                                                     tAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCAGTGCCAACCCCCGTTGTTGGCGGCCAGGTc                                                                                              <  1:153791/147‑1 (MQ=255)
                                                            ggTACCCAGTGGCGAATGAACCGAGCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTAt                                                                                     >  2:323124/1‑149 (MQ=255)
                                                                 ccAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAgg                                                                                >  1:313740/1‑149 (MQ=255)
                                                                 ccAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCCACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGt                                                                               >  8:11415/1‑149 (MQ=255)
                                                                  cAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACt                                                                                                                          <  1:300056/106‑1 (MQ=255)
                                                                  cAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACt                                                                                                                          >  2:300056/1‑106 (MQ=255)
                                                                  cAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTg                                                                              >  5:271866/1‑149 (MQ=255)
                                                                  cAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTg                                                                              <  5:229876/149‑1 (MQ=255)
                                                                   aGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCCACTGCCAACCCCCGTTGTTGGCGGTCAGGTCGCCGTCTATCAGGTg                                                                              <  1:90653/149‑1 (MQ=255)
                                                                         gAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCaag                                                                       >  1:237571/1‑147 (MQ=255)
                                                                          aaTGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATg                                                                       >  5:177145/1‑148 (MQ=255)
                                                                           aTGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGaa                                                                     <  4:392739/149‑1 (MQ=255)
                                                                               aCCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATc                                                                                   <  5:81073/131‑1 (MQ=255)
                                                                               aCCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATc                                                                                   >  6:81073/1‑131 (MQ=255)
                                                                                ccGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGcc                                                              >  5:284957/1‑149 (MQ=255)
                                                                                     ccTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCg                                                         >  7:384820/1‑149 (MQ=255)
                                                                                      cTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGc                                                        >  2:90939/1‑149 (MQ=255)
                                                                                            aTCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAACTAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCTCTTCCAACGCCCGGTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTTCCGCTCGCCCTTgc                                                  <  5:288770/149‑1 (MQ=255)
                                                                                                       gCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAAt                                        >  3:252615/1‑148 (MQ=255)
                                                                                                        cGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAAt                                        <  7:14212/147‑1 (MQ=255)
                                                                                                        cGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATc                                       <  6:271866/148‑1 (MQ=255)
                                                                                                         gACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATcag                                     >  1:320014/1‑149 (MQ=255)
                                                                                                           cTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATcag                                     >  3:30646/1‑147 (MQ=255)
                                                                                                            tgtgAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATcagca                                   <  1:90939/149‑1 (MQ=255)
                                                                                                              tgAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATcagcag                                  <  6:165311/147‑1 (MQ=255)
                                                                                                              tgAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGt                                >  4:159168/1‑149 (MQ=255)
                                                                                                                 aaCCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCGGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTGGCCGTCTATCAGGTGGCGCATGAA‑‑GTGACGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCt                             <  6:131966/149‑1 (MQ=255)
                                                                                                                        gTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGa                       <  6:177145/149‑1 (MQ=255)
                                                                                                                           gAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAgg                     <  6:284957/148‑1 (MQ=255)
                                                                                                                            aaGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGaa                   >  8:321419/1‑149 (MQ=255)
                                                                                                                            aaGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGaa                   >  8:188616/1‑149 (MQ=255)
                                                                                                                            aaGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGaa                   >  2:28378/1‑149 (MQ=255)
                                                                                                                               aTACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGa                       >  5:256203/1‑142 (MQ=255)
                                                                                                                               aTACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGa                       <  6:256203/142‑1 (MQ=255)
                                                                                                                                  cGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATgg              <  2:313740/148‑1 (MQ=255)
                                                                                                                                     aaaTAGGGCACCGCGTCGGTGGCGGTGGAAGGAGTCGACTTCCAACCAGAGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGTA‑‑GTGCCGCTCGGCCTTGCGCCAGTGAATCAGCAGGTTCTTGCTGAGGAACATGGCCAc          <  4:30646/149‑1 (MQ=255)
                                                                                                                                        tAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGAc       <  1:28378/149‑1 (MQ=255)
                                                                                                                                         ccggCACCGACGCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGaa                   <  1:335342/134‑1 (MQ=255)
                                                                                                                                         ccggCACCGACGCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGaa                   >  2:335342/3‑136 (MQ=255)
                                                                                                                                          gggCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAATCGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGa        <  1:323124/148‑1 (MQ=255)
                                                                                                                                           ggCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCAt    >  2:29858/1‑149 (MQ=255)
                                                                                                                                            gCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGTCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTgcgc                                                >  5:90400/1‑104 (MQ=255)
                                                                                                                                            gCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGTCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTgcgc                                                <  6:90400/104‑1 (MQ=255)
                                                                                                                                            gCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGa                                                                      >  7:221817/1‑84 (MQ=255)
                                                                                                                                            gCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGa                                                                      <  8:221817/84‑1 (MQ=255)
                                                                                                                                             cACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATgc  <  1:134375/149‑1 (MQ=255)
                                                                                                                                                  cgTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCtt                                                    <  8:21418/94‑1 (MQ=255)
                                                                                                                                                  cgTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCtt                                                    >  7:21418/1‑94 (MQ=255)
                                                                                                                                                   gTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCAt    >  4:169646/1‑141 (MQ=255)
                                                                                                                                                   gTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCAt    <  3:169646/141‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAAGCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAA‑‑GTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGC  >  NC_002947/857493‑857781

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: