Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,450,415 C→A 7.2% R223L (CGC→CTC)  PP_2146 ← SNF2/RAD54 family helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,450,4150CA 7.2% 217.7 / 4.2 69R223L (CGC→CTC) PP_2146SNF2/RAD54 family helicase
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base A (2/3);  total (33/36)
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.64e-01

GCAGCGCAGCATCGCTTCGCCGCTGTACCACCGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTACGAAGCGCGGTGGCTCGTAGAGGAAATCGTAAAGG  >  NC_002947/2450267‑2450561
                                                                                                                                                     |                                                                                                                                                  
gcagcgcagcATCGCTTCGCCGCTGTACCACCGCCCCCGGTTATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGgcgc                                                                                                                                                    >  6:351450/1‑149 (MQ=255)
     gcagcaTCGCTTCGCCGCTGTACCACCGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAACCATCCGCCCGGTTTCGAGGGCGCGCTg                                                                                                                                                <  1:201542/149‑1 (MQ=255)
       agcaTCGCATCGCCGCTGTACCACCGCCCCCGTTAATTTCCTTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGATGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGaa                                                                                                                                             <  4:372268/149‑1 (MQ=255)
        gcaTCGCTTCGCCGCTGAACCACCGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGAGCTgttc                                                                                                                                             >  4:255404/1‑145 (MQ=255)
         caTCGCTTCGCCGCTGTACCACCGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGAGCTgttcca                                                                                                                                           >  8:19593/1‑144 (MQ=255)
               ttCGCCGCTGTACCACCGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGAGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAC‑CATCCGCCCGGGTTCGAGGGCGCGCTGGAACAGCTCGc                                                                                                                                     >  7:145884/1‑149 (MQ=255)
                  gccgcTGTACCACCGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCAc                                                                                                                                  >  6:315229/1‑149 (MQ=255)
                  gccgcTGTACCACCGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGAGCTGTTCCAGCGCGccct                                                                                                                                  <  7:19593/149‑3 (MQ=255)
                  gccgcTGTACCACCGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGAGCTGTTCCAGCGCGccct                                                                                                                                  <  4:42779/149‑3 (MQ=255)
                  gccgcTGAACCACCGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGAGCTGTTCCAGCGCGccc                                                                                                                                   <  3:255404/148‑2 (MQ=255)
                   ccgcTGTACCACCGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGTACAGCTCGCCAcc                                                                                                                                 >  6:32775/1‑149 (MQ=255)
                    cgcTGTACCACCGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCt                                                                                                                                >  8:164375/1‑149 (MQ=255)
                           ccaccGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGc                                                                                                                         >  3:320510/1‑149 (MQ=255)
                                    ccGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGcc                                                                                                                                    >  2:268978/1‑129 (MQ=255)
                                    ccGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGcc                                                                                                                                    <  1:268978/129‑1 (MQ=255)
                                     cGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCg                                                                                                               <  1:294436/149‑1 (MQ=255)
                                          aTTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGc                                                                                                          >  8:8847/1‑149 (MQ=255)
                                             tCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCcgc                                                                                                       >  1:293026/1‑149 (MQ=255)
                                                gTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTTGAGCGGTAGCCGCGAt                                                                                                    >  1:328089/1‑149 (MQ=255)
                                                   gTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGCACAGCTCGCCACCTTCATGGCCCTACAGTCGGTAGCCGCGATCCt                                                                                                 >  6:156048/1‑149 (MQ=255)
                                                    tCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCt                                                                                                 <  7:145549/148‑1 (MQ=255)
                                                    tCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCt                                                                                                 <  5:166923/148‑1 (MQ=255)
                                                    tCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCt                                                                                                 <  2:293026/148‑1 (MQ=255)
                                                      cAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGgc                                                                                              >  4:310269/1‑149 (MQ=255)
                                                       aGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCGTTCAGTCGGTAGCCGCGATCCTGgcc                                                                                             >  6:229286/1‑149 (MQ=255)
                                                       aGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGgcc                                                                                             >  1:191476/1‑149 (MQ=255)
                                                           gTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGgccgccg                                                                                         <  6:263660/149‑1 (MQ=255)
                                                             aCCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGCCCCTTCAGTCGGTAGCCGCGATCCTGgccgccgcc                                                                                       >  2:360699/1‑149 (MQ=255)
                                                                      aaCTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCaa                                                                              <  4:343687/149‑1 (MQ=255)
                                                                      aaCTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCaa                                                                              <  5:351450/149‑1 (MQ=255)
                                                                       aCTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAg                                                                             <  5:315229/149‑1 (MQ=255)
                                                                       aCTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAg                                                                             <  1:222967/149‑1 (MQ=255)
                                                                        cTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAgg                                                                            <  3:153782/149‑1 (MQ=255)
                                                                           ggCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGcc                                                                          <  7:127626/148‑1 (MQ=255)
                                                                            gCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGTACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCga                                                                        <  5:32775/149‑1 (MQ=255)
                                                                            gCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCga                                                                        <  4:292080/149‑1 (MQ=255)
                                                                             cGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCgag                                                                       <  3:47053/149‑1 (MQ=255)
                                                                               tGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGagag                                                                     >  7:300926/1‑149 (MQ=255)
                                                                                  cgcAGCGGGCCGTCATGCAGCAGGGGGAGTGCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGTCGCTCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCTAGAGCAg                                                                  <  7:343734/149‑1 (MQ=255)
                                                                                    cACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGgc                                                                <  7:164375/149‑1 (MQ=255)
                                                                                       cGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCa                                                             <  1:110072/149‑1 (MQ=255)
                                                                                       cGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCa                                                             <  3:310269/149‑1 (MQ=255)
                                                                                       cGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCa                                                             <  4:210816/149‑1 (MQ=255)
                                                                                        gggCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGgcg                                                               >  8:266496/1‑146 (MQ=255)
                                                                                            cTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCa                                                             <  5:164113/144‑1 (MQ=255)
                                                                                            cTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCa                                                             >  6:164113/1‑144 (MQ=255)
                                                                                            cTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGAAAATTGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGAAGAATg                                                        >  7:330926/1‑149 (MQ=255)
                                                                                                aTGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCgg                                                     <  1:360699/148‑1 (MQ=255)
                                                                                                 tGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTTGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCAGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGTGGAt                                                   <  2:191476/149‑1 (MQ=255)
                                                                                                 tGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGAt                                                   >  8:9146/1‑149 (MQ=255)
                                                                                                     gCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTc                                                <  5:156048/148‑1 (MQ=255)
                                                                                                     gCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGAGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCt                                               >  1:271194/1‑149 (MQ=255)
                                                                                                               gTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGa                                                                                   >  6:225080/1‑103 (MQ=255)
                                                                                                               gTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGtt                                      <  8:300926/148‑1 (MQ=255)
                                                                                                               gTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTa                                     <  4:320510/149‑1 (MQ=255)
                                                                                                                 cGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTAc                                    >  7:143119/1‑148 (MQ=255)
                                                                                                                   gtgtCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTACg                                   >  4:310067/1‑147 (MQ=255)
                                                                                                                       cGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCaa                                                                              <  3:153357/100‑1 (MQ=255)
                                                                                                                       cGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCaa                                                                              >  4:153357/1‑100 (MQ=255)
                                                                                                                           gaa‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTc                                             <  6:47103/129‑1 (MQ=255)
                                                                                                                           gaa‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTc                                             >  5:47103/1‑129 (MQ=255)
                                                                                                                           gaa‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTACGAAGCGCGGTg                         >  4:235471/1‑149 (MQ=255)
                                                                                                                               cATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTACGAAGCGCGGTGGCt                      >  6:361626/1‑149 (MQ=255)
                                                                                                                                   cGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTc                                             <  8:164003/122‑1 (MQ=255)
                                                                                                                                   cGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTc                                             >  7:164003/1‑122 (MQ=255)
                                                                                                                                     cccGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCGTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTACGAAGCGCGGTGGCTCGTa                  >  7:267011/1‑147 (MQ=255)
                                                                                                                                     cccGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTACGAAGCGCGGTGGCTCGTaga                <  8:143119/149‑1 (MQ=255)
                                                                                                                                      ccGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTACGAAGCGCGGTGGCTCGTagag               <  2:328089/149‑1 (MQ=255)
                                                                                                                                        ggTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTACGAAGCGCGGTGGCTCGTagagc              >  7:262674/1‑147 (MQ=255)
                                                                                                                                        ggTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTACGAAGCGCGGTGGCTCGTAGAgg              >  4:341480/1‑148 (MQ=255)
                                                                                                                                                   cgcgcTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTACGAAGCGCGGTGGCTCGTAGAGGAAATCGTAAAgg  >  2:325771/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCAGCGCAGCATCGCTTCGCCGCTGTACCACCGCCCCCGGTAATTTCCGTTGTCCAGGCGTACCCAGCGAAACTCGGCGTGGCGCACCGGGCCTTCATGAAGCAGGGGGAGGTCGGTGTCGAAGAA‑CATCCGCCCGGTTTCGAGGGCGCGCTGGAACAGCTCGCCACCTTCATGGCCCTTCAGTCGGTAGCCGCGATCCTGGCCGCCGCCCTGATTCAAGGCCGAGAGCAGGCGCAGAATGCGGATGTCTTCGTCTGTTACGAAGCGCGGTGGCTCGTAGAGGAAATCGTAAAGG  >  NC_002947/2450267‑2450561

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: ‑