Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,631,188 C→G 11.1% Q204E (CAG→GAG)  lon‑II → DNA‑binding ATP‑dependent protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,631,1880CG11.1% 228.5 / 15.3 80Q204E (CAG→GAG) lon‑IIDNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base C (31/40);  new base G (3/6);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

ATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGAAAAAGC  >  NC_002947/2631041‑2631335
                                                                                                                                                   |                                                                                                                                                   
aTCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAg                                                                                                                                                   >  2:2170642/1‑150 (MQ=255)
    acgaGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCg                                                                                                                                               >  2:2343513/1‑150 (MQ=255)
      gaGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTGCAGGAGGTCGat                                                                                                                                             <  1:1218317/150‑2 (MQ=255)
      gaGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCCAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTGCAGGAGGTCGat                                                                                                                                             <  2:425486/150‑2 (MQ=255)
        gCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGaaaa                                                                                                                                           <  2:2139962/150‑1 (MQ=255)
         ccAGGGCGCCTTGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGAGGCCGAAATCGACCTGCAGGAGGTCGatttc                                                                                                                                          <  2:1022196/150‑5 (MQ=255)
         ccAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTGCAGGAGGTCGatttc                                                                                                                                          <  2:2302480/150‑5 (MQ=255)
         ccAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCATCGACTTGACTACCCGTGTCGAGCAAGTGCCAGCAGTGGGTGTGGCCGAAAGCGACCTCCTGCAGGTCGAAAAg                                                                                                                                          >  1:1738520/1‑150 (MQ=255)
         ccAGGGCGCCTGGTCGACACCATGGCCACGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTGCAGGAGGTCGatttc                                                                                                                                          <  2:2221799/150‑5 (MQ=255)
              gcgcCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCAt                                                                                                                                     >  2:2336048/1‑150 (MQ=255)
                  cTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAgc                                                                                                                                         >  1:1040439/1‑142 (MQ=255)
                  cTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAgc                                                                                                                                         <  2:1040439/142‑1 (MQ=255)
                       cGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCg                                                                                                                            <  1:2336048/150‑1 (MQ=255)
                        gACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCgg                                                                                                                           <  1:1645542/150‑1 (MQ=255)
                           acCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTGCAGGAGGTCGa                                                                                                                                              >  1:2216226/1‑128 (MQ=255)
                           acCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTGCAGGAGGTCGa                                                                                                                                              <  2:2216226/128‑1 (MQ=255)
                             cATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCaa                                                                                                                      >  2:501703/1‑150 (MQ=255)
                             cATGGCCGCGCACATGGCCCGGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACGTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCaa                                                                                                                      >  2:1274362/1‑150 (MQ=255)
                              aTGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATTCCGAAATCGCCCTCCGGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCaag                                                                                                                     <  2:1014657/150‑1 (MQ=255)
                              aTGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCaag                                                                                                                     >  2:1992718/1‑150 (MQ=255)
                                    gcgcACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAg                                                                                                               >  2:738945/1‑150 (MQ=255)
                                     cgcACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGa                                                                                                              <  1:1509712/150‑1 (MQ=255)
                                      gcACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGAt                                                                                                             >  2:749352/1‑150 (MQ=255)
                                                     gATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAgc                                                                                                                                         <  1:652229/107‑1 (MQ=255)
                                                     gATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAgc                                                                                                                                         >  2:652229/1‑107 (MQ=255)
                                                        cGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGa                                                                                           <  1:749352/150‑1 (MQ=255)
                                                        cGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGa                                                                                           <  1:2397073/150‑1 (MQ=255)
                                                           gCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGta                                                                                        >  1:1961823/1‑150 (MQ=255)
                                                           gCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGta                                                                                        >  1:1961835/1‑150 (MQ=255)
                                                               aaGCAGGAAATTCTTGAAATCGTCGACTGGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGtactac                                                                                    <  2:118752/150‑1 (MQ=255)
                                                                aGCAGGTAATTCTTTAAATCGTCGACTTGACTGCCCGTGTCGAGCACGGGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAAGTCGAAAAGCGCATTCGCGGCCGGTTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTAcc                                                                                   <  1:396175/150‑1 (MQ=255)
                                                                       aaTTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGa                                                                            <  2:119999/150‑1 (MQ=255)
                                                                        aTTCTTGAAATCGACGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAg                                                                           <  1:1696663/150‑1 (MQ=255)
                                                                         ttCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGNTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGc                                                                          <  1:1486453/150‑1 (MQ=255)
                                                                         ttCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGc                                                                          <  1:501703/150‑1 (MQ=255)
                                                                          tCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCa                                                                         >  2:480648/1‑150 (MQ=255)
                                                                                 aTCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAATAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAg                                                                  <  1:301445/150‑1 (MQ=255)
                                                                                 aTCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAg                                                                  <  1:2170642/150‑1 (MQ=255)
                                                                                    gtcgACTTGACTACCCGTGTCGAGCAAGTGCTGGCGCTGCTGGATGCCGAAATCGACATCCAGCTGGTAGGAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGAAGCCAGCGCGAGGCCTACCTGCAAGAGCAGATGAAGGcc                                                               >  2:817141/1‑150 (MQ=255)
                                                                                        aCTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTc                                                           >  1:2328475/1‑150 (MQ=255)
                                                                                        aCTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTc                                                           >  2:1228733/1‑150 (MQ=255)
                                                                                          ttGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCag                                                         >  2:1501806/1‑150 (MQ=255)
                                                                                             aCTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCagaaa                                                      <  1:2043503/150‑1 (MQ=255)
                                                                                               tACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGTGCATTCagaaaga                                                    >  1:1099953/1‑150 (MQ=255)
                                                                                               tACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCagaaaga                                                    >  2:1924107/1‑150 (MQ=255)
                                                                                                aCCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCagaaagaa                                                   <  2:1694056/150‑1 (MQ=255)
                                                                                                     tgtCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTgg                                              <  1:1501806/150‑1 (MQ=255)
                                                                                                      gtCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGc                                             <  1:1358726/150‑1 (MQ=255)
                                                                                                            cACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTACTCCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGt                                       <  1:738945/150‑1 (MQ=255)
                                                                                                               gTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGAc                                    <  2:1858719/150‑1 (MQ=255)
                                                                                                               gTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGAc                                    <  2:1859244/150‑1 (MQ=255)
                                                                                                               gTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGCACTACCTGAATGAGCAGACGAAGGCCATTCAGAAAGAACTTGGCGATGGTGAc                                    >  1:2292314/1‑150 (MQ=255)
                                                                                                                   tGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAg                                <  1:480648/150‑1 (MQ=255)
                                                                                                                     gcgcTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGACGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGc                              <  2:2292314/150‑1 (MQ=255)
                                                                                                                         tgctgGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAg                                                                  <  2:1456165/110‑1 (MQ=255)
                                                                                                                         tgctgGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAg                                                                  >  1:1456165/1‑110 (MQ=255)
                                                                                                                         tgctgGATGCAGAAATCGACCTCCTGGAGGTCGAAAATCGCATTCTTGGCCGGGTCAAGGAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCCGAAAGAACGTGGCGATGGTGACGAAGGCcaca                          >  2:760009/1‑150 (MQ=255)
                                                                                                                           ctgGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAAt                        <  2:1829128/150‑1 (MQ=255)
                                                                                                                                 gCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTc                  <  2:1961835/150‑1 (MQ=255)
                                                                                                                                 gCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTc                  <  2:1961823/150‑1 (MQ=255)
                                                                                                                                   gcttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGa                >  2:213238/6‑150 (MQ=255)
                                                                                                                                   cGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTAcc                                                                                   >  2:167835/1‑83 (MQ=255)
                                                                                                                                   cGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTAcc                                                                                   <  1:167835/83‑1 (MQ=255)
                                                                                                                                   cGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGa                >  1:959318/1‑150 (MQ=255)
                                                                                                                                   cGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGa                >  1:1619688/1‑150 (MQ=255)
                                                                                                                                   cGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACATGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGa                <  2:164722/150‑1 (MQ=255)
                                                                                                                                    gAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGaa               <  1:1924107/150‑1 (MQ=255)
                                                                                                                                    gAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGaa               >  2:1563006/1‑150 (MQ=255)
                                                                                                                                     aaaTCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCa                                                                                                                <  2:1686583/52‑1 (MQ=255)
                                                                                                                                     aaaTCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCa                                                                                                                >  1:1686583/1‑52 (MQ=255)
                                                                                                                                             ctcctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAg                                >  1:2344162/1‑124 (MQ=255)
                                                                                                                                             ctcctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAg                                <  2:2344162/124‑1 (MQ=255)
                                                                                                                                             ctcctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGaaa      >  2:379154/1‑150 (MQ=255)
                                                                                                                                               cctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGaaaaa    <  1:565960/150‑1 (MQ=255)
                                                                                                                                                ctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGTCGATGGAGCGAAGCCAGCGTGAGTACTACCTGAATGGGCAGATGAAGGCCATTCAGAAAGGACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGAAAAAg   <  1:2343513/150‑1 (MQ=255)
                                                                                                                                                ctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGAGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATg                       <  2:294320/130‑1 (MQ=255)
                                                                                                                                                ctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGAGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATg                       >  1:294320/1‑130 (MQ=255)
                                                                                                                                                 tGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGAAAAAgc  <  1:839157/150‑1 (MQ=255)
                                                                                                                                                  gCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCa                                                              <  1:780106/89‑1 (MQ=255)
                                                                                                                                                  gCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCa                                                              >  2:780106/1‑89 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGAAAAAGC  >  NC_002947/2631041‑2631335

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