Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,324,930 A→C 5.7% M506L (ATG→CTG)  PP_4684 → kinase‑like domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,324,9300AC 5.7% 326.7 / 6.8 104M506L (ATG→CTG) PP_4684kinase‑like domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (46/52);  new base C (3/3);  total (49/55)
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.89e-01

Only 100 of 108 total aligned reads displayed.

CCAGTTGGCTGGAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTCCCAACGCCCGCTACACTC  >  NC_002947/5324781‑5325075
                                                                                                                                                     |                                                                                                                                                 
ccAGTTGGCTGGAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGgcagca                                                                                                                                                   <  1:1657994/150‑1 (MQ=255)
    ttGGCTGGAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCAt                                                                                                                                                  >  2:1350270/1‑147 (MQ=255)
    ttGGCTGGAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCAt                                                                                                                                                  <  1:1350270/147‑1 (MQ=255)
    ttGGCTGGAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGa                                                                                                                                               >  2:1128828/1‑150 (MQ=255)
     tGGCTGGAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGAt                                                                                                                                              <  2:495435/150‑1 (MQ=255)
     tGGCTGGAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGAt                                                                                                                                              <  1:1128828/150‑1 (MQ=255)
       gCTGGAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCATGCTGCCGGCAGCCTgggt                                                                                                                                              <  2:742163/148‑3 (MQ=255)
       gCTGGAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCATGCTGCCGGCAGCCTgggt                                                                                                                                              >  1:742163/1‑146 (MQ=255)
        cTGGAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAg                                                                                                                                           <  1:1087129/150‑1 (MQ=255)
         tGGAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAgg                                                                                                                                          >  1:764154/1‑150 (MQ=255)
          ggAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGt                                                                                                                                         <  1:2057265/150‑1 (MQ=255)
            agcagcGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGa                                                                                                                                       <  1:158398/150‑1 (MQ=255)
             gcagcGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGAt                                                                                                                                      <  1:1444681/150‑1 (MQ=255)
             gcagcGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGAt                                                                                                                                      >  1:1816667/1‑150 (MQ=255)
              cagcGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATc                                                                                                                                     <  2:1816667/150‑1 (MQ=255)
                gcgcCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGa                                                                                                                                   <  2:1817140/150‑1 (MQ=255)
                   ccAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCGGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCa                                                                                                                                <  1:2327230/150‑1 (MQ=255)
                   ccAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCa                                                                                                                                >  1:2191335/1‑150 (MQ=255)
                    cAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAg                                                                                                                               >  2:2104797/1‑150 (MQ=255)
                     aGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGTTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGa                                                                                                                              >  1:2036705/1‑150 (MQ=255)
                     aGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGa                                                                                                                              <  2:1822720/150‑1 (MQ=255)
                     aGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGa                                                                                                                              <  1:348205/150‑1 (MQ=255)
                     aGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGa                                                                                                                              <  2:1020364/150‑1 (MQ=255)
                       gCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGAtt                                                                                                                            >  1:792163/1‑150 (MQ=255)
                        ccGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTc                                                                                                                           >  1:2316167/1‑150 (MQ=255)
                             aaCAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCa                                                                                                                       <  2:258532/149‑1 (MQ=255)
                                 gCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCATCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCATGCTGCCGGCAGCCTgggtatc                                                                                                                                           <  1:1488754/125‑6 (MQ=255)
                                   gACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTg                                                                                                                >  2:2054679/1‑150 (MQ=255)
                                     cccGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTgcc                                                                                                              >  1:1053527/1‑150 (MQ=255)
                                      ccGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTc                                                                                                                           <  1:1246163/136‑1 (MQ=255)
                                      ccGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTc                                                                                                                           >  2:1246163/1‑136 (MQ=255)
                                      ccGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTgcct                                                                                                             <  1:1426965/150‑1 (MQ=255)
                                        gTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTgcctgg                                                                                                           >  2:1474947/1‑150 (MQ=255)
                                                 cACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGt                                                                                                  >  1:696905/1‑150 (MQ=255)
                                                 cACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGt                                                                                                  >  1:1768269/1‑150 (MQ=255)
                                                 cACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAAATTCGCTAGCGGTTGCCTGGCCTGGAAGt                                                                                                  >  2:1124289/1‑150 (MQ=255)
                                                 cACGCTGGAGGTGATAAAAACTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGt                                                                                                  >  1:2373023/1‑150 (MQ=255)
                                                  aCGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGtt                                                                                                 >  1:33674/1‑150 (MQ=255)
                                                       ggAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCa                                                                                            >  1:513423/1‑150 (MQ=255)
                                                        gAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCAc                                                                                           <  2:792163/150‑1 (MQ=255)
                                                         aGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAg                                                                                                                      <  2:41446/122‑1 (MQ=255)
                                                         aGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAg                                                                                                                      >  1:41446/1‑122 (MQ=255)
                                                             gATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGAt                                                                                                                                              <  2:1795549/94‑1 (MQ=255)
                                                              aacAAAGCGCAGCAGGCCAACCGTGACATGCCGGATTCCCAGAAGCTGGGGCACATCACCCGCTTCTATACCCATGCTGTCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCt                                                                                     <  2:419220/148‑1 (MQ=255)
                                                              aTCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCt                                                                                     <  2:696905/150‑1 (MQ=255)
                                                                cAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCttt                                                                                   <  2:405491/150‑3 (MQ=255)
                                                                  aGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGtgc                                                                                 <  2:513423/150‑1 (MQ=255)
                                                                   ggCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGtgcc                                                                                >  2:386825/1‑150 (MQ=255)
                                                                   ggCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGtgcc                                                                                >  1:1516400/1‑150 (MQ=255)
                                                                    gCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCATGCTGCCGGCAGCCTgg                                                                                                                                                >  2:1590249/1‑85 (MQ=255)
                                                                    gCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCATGCTGCCGGCAGCCTgg                                                                                                                                                <  1:1590249/85‑1 (MQ=255)
                                                                    gCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGGTGCCTGGCCTGTAAGTTCGTCACCTTCCTGtgcct                                                                               <  1:1124289/150‑1 (MQ=255)
                                                                      tCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTACCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGtgcctgc                                                                             >  2:1866510/1‑150 (MQ=255)
                                                                        agcagGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGccc                                                                           >  2:2151085/1‑150 (MQ=255)
                                                                        agcagGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGccc                                                                           <  2:2215050/150‑1 (MQ=255)
                                                                        agcagGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGccc                                                                           <  1:1474947/150‑1 (MQ=255)
                                                                          cagGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGc                                                                         >  2:1451961/1‑150 (MQ=255)
                                                                              ccAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAAt                                                                     <  2:1303187/150‑1 (MQ=255)
                                                                                aGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATgg                                                                   <  2:2373023/150‑1 (MQ=255)
                                                                                  ccGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTACCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGggc                                                                 <  1:1866510/150‑1 (MQ=255)
                                                                                   cGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGggca                                                                <  2:1053527/150‑1 (MQ=255)
                                                                                   cGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGggca                                                                >  2:1611572/1‑150 (MQ=255)
                                                                                              ggATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTaa                                                                                                    >  2:50177/1‑103 (MQ=255)
                                                                                              ggATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTaa                                                                                                    <  1:50177/103‑1 (MQ=255)
                                                                                               gATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGc                                                    <  2:2299082/150‑1 (MQ=255)
                                                                                                    gcAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGtgcctgc                                                                             >  2:1619325/1‑120 (MQ=255)
                                                                                                    gcAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGtgcctgc                                                                             <  1:1619325/120‑1 (MQ=255)
                                                                                                    gcAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCAcc                                               <  2:764154/150‑1 (MQ=255)
                                                                                                    gcAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCAcc                                               <  1:2054679/150‑1 (MQ=255)
                                                                                                    gcAGGAGCTGGTGCACAGCACCCGAGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCAcc                                               >  2:1599193/1‑150 (MQ=255)
                                                                                                           cTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGAt                                        <  2:1516400/150‑1 (MQ=255)
                                                                                                             ggTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATgg                                      <  1:2354552/150‑1 (MQ=255)
                                                                                                                  aaaGCACACGCATCGATACCAAGGCTACCGGCAGTATGGATCATGTGATCTAGTGGATTCGCCAGCGGTTGCCTGGACTGTAACTTCGTCACCTTCCCCTGCCAGCACGCCAATGGGCAGGTACAGGCCGCTCACCCTACGATGGCaaaa                                 <  2:266478/148‑1 (MQ=255)
                                                                                                                   caGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCTCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCaaaaa                                >  1:1121339/1‑150 (MQ=255)
                                                                                                                   caGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCaaaaa                                <  1:1185298/150‑1 (MQ=255)
                                                                                                                           gcgTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCata                        >  2:1389988/1‑150 (MQ=255)
                                                                                                                                gATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTcc                   <  1:2151085/150‑1 (MQ=255)
                                                                                                                                gATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTcc                   <  1:1686402/150‑1 (MQ=255)
                                                                                                                                 aTACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCt                                                                                     >  1:1225903/1‑83 (MQ=255)
                                                                                                                                 aTACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCt                                                                                     <  2:1225903/83‑1 (MQ=255)
                                                                                                                                 aTACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTccc                  >  1:2398832/1‑150 (MQ=255)
                                                                                                                                  tACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTcc                                                                                      <  1:1148255/81‑1 (MQ=255)
                                                                                                                                  tACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTcc                                                                                      >  2:1148255/1‑81 (MQ=255)
                                                                                                                                   aCCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTc                                                  <  2:146563/116‑1 (MQ=255)
                                                                                                                                   aCCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTc                                                  >  1:146563/1‑116 (MQ=255)
                                                                                                                                   aCCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTCCCaa                <  2:1768269/150‑1 (MQ=255)
                                                                                                                                    cccAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTCCCAAc               <  2:2191335/150‑1 (MQ=255)
                                                                                                                                     ccAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATggg                                                                  <  1:847934/98‑1 (MQ=255)
                                                                                                                                     ccAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATggg                                                                  >  2:847934/1‑98 (MQ=255)
                                                                                                                                     ccAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTCCCAACg              >  2:312166/1‑150 (MQ=255)
                                                                                                                                     ccAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTCCCAACg              >  1:752375/1‑150 (MQ=255)
                                                                                                                                     ccAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTCCCAACg              >  1:668702/1‑150 (MQ=255)
                                                                                                                                      cAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTCCCaa                >  2:157895/1‑147 (MQ=255)
                                                                                                                                      cAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTCCCaa                <  1:157895/147‑1 (MQ=255)
                                                                                                                                          cTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTCCCAACGCCCGc         <  1:1599193/150‑1 (MQ=255)
                                                                                                                                              cGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTcc                                           >  2:2192342/1‑112 (MQ=255)
                                                                                                                                              cGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTcc                                           <  1:2192342/112‑1 (MQ=255)
                                                                                                                                                 cagcaTGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTCCCAACGCCCGCTACActc  <  2:2036705/150‑1 (MQ=255)
                                                                                                                                                  agcaTGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCaaaaa                                >  2:1564654/1‑119 (MQ=255)
                                                                                                                                                  agcaTGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCaaaaa                                <  1:1564654/119‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCAGTTGGCTGGAGCAGCGCCAGGCCGAAAACAGCGACCCGTCGGATGCCACGCTGGAGGTGATCAAGGCTCAGCAGGCCAGCCGTGAGGCGCTGGATGCGCAGGAGCTGGTGCACAGCACCCGCGTCGATACCCAGGCTGCCGGCAGCATGGATCAGGTGATCGAGCAGATTCGCCAGCGGTTGCCTGGCCTGTAAGTTCGTCACCTTCCTGTGCCTGCCCGCCAATGGGCAGGCACAGGCCGCTCACCCTCCGATGGCAAAAAGCCGCATATTTCCCAACGCCCGCTACACTC  >  NC_002947/5324781‑5325075

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: