Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,442,053 C→A 6.0% P383Q (CCG→CAG)  topA → DNA topoisomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,442,0530CA 6.0% 209.3 / 2.7 67P383Q (CCG→CAG) topADNA topoisomerase I
Reads supporting (aligned to +/- strand):  ref base C (37/26);  new base A (2/2);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

CGGCTGACGCCCTGGATATGGCACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGCCGGCGATTTCGAGCTGCGCGCCAAGGGCC  >  NC_002947/2441906‑2442202
                                                                                                                                                   |                                                                                                                                                     
cGGCTGACGCCCTGGATATGGCACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGAGCCTTGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGa                                                                                                                                                     <  1:1788865/150‑1 (MQ=255)
  gCTGACGCCCTGGATATGGCACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGAcc                                                                                                                                                   >  1:1008269/1‑150 (MQ=255)
          ccTGGATATGGCACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGa                                                                                                                                                     >  1:1633643/1‑140 (MQ=255)
          ccTGGATATGGCACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGa                                                                                                                                                     <  2:1633643/140‑1 (MQ=255)
          ccTGGATATGGCACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAg                                                                                                                                           >  1:934867/1‑150 (MQ=255)
          ccTGGATATGGCACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAg                                                                                                                                           >  2:1300123/1‑150 (MQ=255)
           cTGGATATGGCACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCCGCAACGAGGGCGCCCAGGAGGCGCAAGAAGAGACTTCTCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGc                                                                                                                                          >  1:888126/1‑150 (MQ=255)
                 ttGGGACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATg                                                                                                                                    <  2:941980/149‑1 (MQ=255)
                   ggCACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGa                                                                                                                                  >  1:198661/1‑150 (MQ=255)
                   ggCACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGa                                                                                                                                  <  2:1008269/150‑1 (MQ=255)
                     cACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGc                                                                                                                                <  2:1532553/150‑1 (MQ=255)
                      aCGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACACGCTAGCGGCATGCATTGTTTACGGCAGCAAGGGAGTCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATAACCACCCGACCAAGCTCAGCGGCATGGAGCg                                                                                                                               >  1:799116/1‑150 (MQ=255)
                          aGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCAGACCAAGCTCAGCGGCATGGAGCGTGAc                                                                                                                           >  2:1788097/1‑150 (MQ=255)
                             tACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGTAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGcc                                                                                                                        >  2:1437831/1‑150 (MQ=255)
                             tACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGcc                                                                                                                        >  2:1117205/1‑150 (MQ=255)
                               cATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGa                                                                                                                      >  1:124721/1‑150 (MQ=255)
                                  cGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAgcg                                                                                                                   >  1:1805524/1‑150 (MQ=255)
                                       gggAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGt                                                                                                              >  1:658055/1‑150 (MQ=255)
                                         gAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTAc                                                                                                            <  2:1669008/150‑1 (MQ=255)
                                            ttCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAg                                                                                                         >  2:1439567/1‑150 (MQ=255)
                                                 cAAGCAGTACCTGCCGGAGGTGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGAt                                                                                                    >  1:839871/1‑150 (MQ=255)
                                                 cAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGAt                                                                                                    >  1:1136750/1‑150 (MQ=255)
                                                   aGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCt                                                                                                  <  2:1136750/150‑1 (MQ=255)
                                                    gCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTg                                                                                                 >  1:1130725/1‑150 (MQ=255)
                                                    gCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGAGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCTACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTg                                                                                                 >  2:689084/1‑150 (MQ=255)
                                                     cAGGACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCCCGCAAGAAGCGCACGATGCGCTCCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTgg                                                                                                >  2:505469/1‑150 (MQ=255)
                                                         aCCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGcc                                                                                            <  2:1805524/150‑1 (MQ=255)
                                                           cTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAg                                                                                          <  1:393028/150‑1 (MQ=255)
                                                                ggAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCt                                                                                     <  2:1734395/150‑1 (MQ=255)
                                                                 gAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTg                                                                                    <  1:1510219/150‑1 (MQ=255)
                                                                   ggCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGc                                                                                  >  1:954328/1‑150 (MQ=255)
                                                                   ggCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAATGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGAGGCATGGAGCGTGACGCCGAGCGCCGGTACGAGCTGATCTGGCGCCGGGTCCTGGc                                                                                  >  1:1293967/1‑150 (MQ=255)
                                                                    gcgcCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCg                                                                                 >  2:1247418/1‑150 (MQ=255)
                                                                     cgcCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGt                                                                                <  2:954328/150‑1 (MQ=255)
                                                                     cgcCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGt                                                                                <  1:1351348/150‑1 (MQ=255)
                                                                      gcCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTg                                                                               >  2:1292221/1‑150 (MQ=255)
                                                                        cACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGcc                                                                             >  2:1404259/1‑150 (MQ=255)
                                                                        cACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCAGACCAAGCTCAGCGGCATGGAGCGTGGCGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGcc                                                                             >  2:640684/1‑150 (MQ=255)
                                                                         aCTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCa                                                                            <  1:1423429/150‑1 (MQ=255)
                                                                         aCTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCAGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCa                                                                            <  1:1788097/150‑1 (MQ=255)
                                                                           tGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAGGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGa                                                                          <  2:151228/150‑1 (MQ=255)
                                                                            ggTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCaa                                                                                                                                                 <  1:1200543/78‑1 (MQ=255)
                                                                            ggTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCaa                                                                                                                                                 >  2:1200543/1‑78 (MQ=255)
                                                                               ttACGGCAGTAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATgcc                                                                      <  1:1437831/150‑1 (MQ=255)
                                                                                 aCGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATgccgc                                                                    <  1:333596/150‑1 (MQ=255)
                                                                                   ggCAGCAAGGAGGGCGTCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCa                                                                  <  2:1293967/150‑1 (MQ=255)
                                                                                    gcagcaAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAg                                                                 >  2:1523026/1‑150 (MQ=255)
                                                                                        caAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCa                                                             >  2:900551/1‑150 (MQ=255)
                                                                                         aaGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCt                                                                                                                <  2:1091186/98‑1 (MQ=255)
                                                                                         aaGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCt                                                                                                                >  1:1091186/1‑98 (MQ=255)
                                                                                              gggCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCt                                                       <  2:839871/150‑1 (MQ=255)
                                                                                               ggCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTg                                                      >  2:482254/1‑150 (MQ=255)
                                                                                               ggCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTAAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATCACGCCAGCGCAGTAcccc                                                      >  2:764552/1‑148 (MQ=255)
                                                                                                  gcCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTcc                                                   <  1:680461/150‑1 (MQ=255)
                                                                                                         aggCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCAGACCAAGCTCAGCGGCATGGAGCGTGGCGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCg                                            <  1:640684/150‑1 (MQ=255)
                                                                                                           gcgcACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAgcgc                                                                                                                  <  2:169801/78‑1 (MQ=255)
                                                                                                           gcgcACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAgcgc                                                                                                                  >  1:169801/1‑78 (MQ=255)
                                                                                                             gcACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGa                                                                          >  2:715637/1‑116 (MQ=255)
                                                                                                             gcACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGa                                                                          <  1:715637/116‑1 (MQ=255)
                                                                                                              cACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAgcgc                                                                                                                  <  2:1233942/75‑1 (MQ=255)
                                                                                                              cACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAgcgc                                                                                                                  >  1:1233942/1‑75 (MQ=255)
                                                                                                                 gAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTc                                    >  1:1913297/1‑150 (MQ=255)
                                                                                                                 gAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTc                                    >  2:27331/1‑150 (MQ=255)
                                                                                                                       aTTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGcc                              >  2:263417/1‑150 (MQ=255)
                                                                                                                        ttCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGCCg                             <  2:1130725/150‑1 (MQ=255)
                                                                                                                          cGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGCCGGc                           >  2:1120083/1‑150 (MQ=255)
                                                                                                                                      gTGGATACGCACGCGACGAAGCTCAGCGGCAGGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTGACCGTCGCTGCCGGCGATTTCGAGCTg               >  2:1012877/1‑150 (MQ=255)
                                                                                                                                              ccacccGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGCCGGCGATTTCGAGCTGCGCGCCaa       >  1:1048596/1‑150 (MQ=255)
                                                                                                                                               cacccGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGgcg                                                                                              <  1:901782/62‑1 (MQ=255)
                                                                                                                                               cacccGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGgcg                                                                                              >  2:901782/1‑62 (MQ=255)
                                                                                                                                               cacccGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGCCTGCGATTTCGAGCTGCGCGCCAAg      >  2:1408535/1‑150 (MQ=255)
                                                                                                                                               cacccGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGCCGGCGATTTCGAGCTGCGCGCCAAg      <  1:27331/150‑1 (MQ=255)
                                                                                                                                                acccGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTTATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGCCGGCGATTTCGAGCTGCGCGCCAAgg     >  1:1397578/1‑150 (MQ=255)
                                                                                                                                                acccGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGCCGGCGATTTCGAGCTGCGCGCCAAgg     <  1:505469/150‑1 (MQ=255)
                                                                                                                                                acccGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCCAGGTACGCGTCCACCAGCGTTACCGTCGCTGCCGGCGATTTCGAGCTGCGCGCCAAgg     >  2:1791515/1‑150 (MQ=255)
                                                                                                                                                 cccGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGCCGGCGATTTCGAGCTGCGCGCCAAggg    >  2:909438/1‑150 (MQ=255)
                                                                                                                                                  ccGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGCCGGCGATTTCGAGCTGCGCGCCAAGGGc   <  1:1117205/150‑1 (MQ=255)
                                                                                                                                                   cGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGCCGGCGATTTCGAGCTGCGCGCCAAGGGcc  <  1:1523026/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGGCTGACGCCCTGGATATGGCACGCAGCTACATCGAACGGGAGTTCGGCAAGCAGTACCTGCCGGAGGCGCCACTGGTTTACGGCAGCAAGGAGGGCGCCCAGGAGGCGCACGAAGCGATTCGGCCTTCTGACGTGAATACCCACCCGACCAAGCTCAGCGGCATGGAGCGTGACGCCGAGCGCCTGTACGAGCTGATCTGGCGCCAGTTCCTGGCGTGCCAGATGCCGCCAGCGCAGTACCTGTCCACCAGCGTTACCGTCGCTGCCGGCGATTTCGAGCTGCGCGCCAAGGGCC  >  NC_002947/2441906‑2442202

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