Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,129,330 G→T 7.3% Q692K (CAG→AAG) ‡ recQ ← ATP‑dependent DNA 3'‑5' helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,129,3300GT 7.3% 215.6 / 4.6 68Q692K (CAG→AAG) ‡recQATP‑dependent DNA 3'‑5' helicase
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base T (3/2);  total (35/33)
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.63e-01

TACGCGTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCCAGCAGGC  >  NC_002947/5129182‑5129470
                                                                                                                                                    |                                                                                                                                             
tACGCGTTCGATAATCGTCAAAAACCATGACAAAGGTCAGCACATGGCGCTTCTGCCTTCGATCATGCATGCGTTTACAATTCGAATTCCGAAGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCCCGCCCCGAACTg                                                                                                                                               >  5:323397/1‑149 (MQ=255)
  cgcgTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCg                                                                                                                                             <  3:34795/149‑1 (MQ=255)
   gcgTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCgg                                                                                                                                            <  8:177105/149‑1 (MQ=255)
    cgTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCggg                                                                                                                                           <  5:202342/149‑1 (MQ=255)
    agTTCGATAATCTTCAAAAAGCATGACAAAGGTCAGCTGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTAGAATTAGAATTCCGCCGCGAGTGCTGCGCGCACGCAATTGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCggg                                                                                                                                           <  6:323397/148‑1 (MQ=255)
     gTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCgggg                                                                                                                                          <  5:196295/149‑1 (MQ=255)
        cGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCaa                                                                                                                                       >  8:5039/1‑149 (MQ=255)
         gATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAAt                                                                                                                                      >  1:6983/1‑149 (MQ=255)
              tCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCggg                                                                                                                                           >  3:21935/1‑139 (MQ=255)
              tCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCggg                                                                                                                                           <  4:21935/139‑1 (MQ=255)
                   aaaaaGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACAcc                                                                                                                            <  5:245681/149‑1 (MQ=255)
                    aaaaGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGCAGTTCGGGGcgcgg                                                                                                                                     <  4:129728/139‑5 (MQ=255)
                    aaaaGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGCAGTTCGGGGcgcgg                                                                                                                                     >  3:129728/1‑135 (MQ=255)
                      aaGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGGCGACACCGGt                                                                                                                         <  6:19123/149‑1 (MQ=255)
                          aTGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCa                                                                                                                     <  4:280414/149‑1 (MQ=255)
                             aCAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCt                                                                                                                 >  5:28078/1‑149 (MQ=255)
                             aCAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCt                                                                                                                 >  8:134415/1‑149 (MQ=255)
                                  ggTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCAccc                                                                                                            >  1:190388/1‑149 (MQ=255)
                                     cagcagATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCg                                                                                                         >  3:172453/1‑149 (MQ=255)
                                      agcagATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCGCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGc                                                                                                        >  2:319444/1‑149 (MQ=255)
                                           aTGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCCCGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGCCGCCGGTGGCAGCTCACCCTCGCCATc                                                                                                    >  4:59568/1‑148 (MQ=255)
                                                gcTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTgg                                                                                                                       >  4:115694/1‑125 (MQ=255)
                                                gcTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTgg                                                                                                                       >  4:120228/1‑125 (MQ=255)
                                                gcTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTgg                                                                                                                       <  3:115694/125‑1 (MQ=255)
                                                gcTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTgg                                                                                                                       <  3:120228/125‑1 (MQ=255)
                                                 cTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAg                                                                                             >  1:288396/1‑149 (MQ=255)
                                                    cTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGc                                                                                           <  1:125604/148‑1 (MQ=255)
                                                     tGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGt                                                                                         >  3:292519/1‑149 (MQ=255)
                                                     tGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGt                                                                                         >  6:69154/1‑149 (MQ=255)
                                                      gCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaaa                                                                                              >  1:254423/1‑143 (MQ=255)
                                                      gCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaaa                                                                                              <  2:254423/143‑1 (MQ=255)
                                                      gCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTc                                                                                        >  8:233323/1‑149 (MQ=255)
                                                       ccTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGACTCACCCTCGCCATCAAGAAAGGc                                                                                           <  7:47044/146‑1 (MQ=255)
                                                         ttCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGCGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaaa                                                                                              <  7:85784/140‑1 (MQ=255)
                                                         ttCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGCGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaaa                                                                                              >  8:85784/1‑140 (MQ=255)
                                                           cGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGc                                                                                    <  7:5039/148‑1 (MQ=255)
                                                            gATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCAc                                                                                  <  7:134415/149‑1 (MQ=255)
                                                              tCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTcc                                                                                       >  3:117683/1‑142 (MQ=255)
                                                              tCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTcc                                                                                       <  4:117683/142‑1 (MQ=255)
                                                               ccTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGCCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTc                                                                               <  2:190388/149‑1 (MQ=255)
                                                                                ttCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACTTCCATCATCAGGTCTTCCg                                                              >  5:96119/1‑149 (MQ=255)
                                                                                ttCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTACATCATCAGGTCTTcag                                                              >  4:215400/1‑147 (MQ=255)
                                                                                ttCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAAATGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCg                                                              >  3:59651/1‑149 (MQ=255)
                                                                                 tagAATACCGCCGACAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACAGCGGGGCAATGGACGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTACTGCACCTCCCTCATCAGGTCTTCCgg                                                             >  3:36637/3‑149 (MQ=255)
                                                                                   gAATTCCGCCTCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGGAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCa                                                           <  4:36637/149‑1 (MQ=255)
                                                                                           ccgccAGTGCTGCGCGCACGCAATAGAGTACGACCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGTCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGTTCTTCCGGCAGGTCGATg                                                   <  2:288396/149‑1 (MQ=255)
                                                                                             gccAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGc                                                 <  4:292519/149‑1 (MQ=255)
                                                                                             gccAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGc                                                 >  4:42133/1‑149 (MQ=255)
                                                                                              ccAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCa                                                                                                                     >  5:167448/1‑81 (MQ=255)
                                                                                              ccAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCa                                                                                                                     <  6:167448/81‑1 (MQ=255)
                                                                                              ccAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGcc                                                >  4:205535/1‑149 (MQ=255)
                                                                                                    ctgcGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAAATGCGGTGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTcc                                          <  4:59651/149‑1 (MQ=255)
                                                                                                               caATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCTCCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTc                               <  2:6983/149‑1 (MQ=255)
                                                                                                               caATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTc                               <  4:172453/149‑1 (MQ=255)
                                                                                                                aaTAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTcc                              >  3:117385/1‑149 (MQ=255)
                                                                                                                 aTAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCt                             <  7:233323/149‑1 (MQ=255)
                                                                                                                            ccTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAAc                  >  1:191859/1‑149 (MQ=255)
                                                                                                                            ccTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAAc                  >  7:296274/1‑149 (MQ=255)
                                                                                                                              tCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATAGCCTGGTCCAGACTCAGCTCCTggcgg                        >  3:229978/1‑141 (MQ=255)
                                                                                                                              tCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATAGCCTGGTCCAGACTCAGCTCCTggcgg                        <  4:229978/141‑1 (MQ=255)
                                                                                                                               ccGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGcc               >  4:256289/1‑149 (MQ=255)
                                                                                                                                  gCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCg            <  8:296274/149‑1 (MQ=255)
                                                                                                                                  gCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCg            <  3:205535/149‑1 (MQ=255)
                                                                                                                                   cACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCgg           <  6:28078/149‑1 (MQ=255)
                                                                                                                                     ccattgcCCCGCAGTTCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCg            >  1:259263/6‑146 (MQ=255)
                                                                                                                                     ccattgcCCCGCAGTTCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCg            <  2:259263/141‑1 (MQ=255)
                                                                                                                                      ccGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGCCTCAGCTCCTGGCGGCCCAACGCCTCGGc          >  5:154204/1‑147 (MQ=255)
                                                                                                                                      cattgcCCCGCAGTTCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGc          >  2:59484/5‑147 (MQ=255)
                                                                                                                                             ccGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCCAGCAGGc  >  6:63373/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
TACGCGTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCCAGCAGGC  >  NC_002947/5129182‑5129470

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: