Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,898,158 C→T 5.0% L114L (CTG→TTG)  PP_2552 → DOPA decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,898,1580CT 5.0% 362.8 / 3.5 121L114L (CTG→TTG) PP_2552DOPA decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (52/62);  new base T (3/4);  total (55/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

ATGCCCGGCCTGTCCCATTGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAACGCGCCACTGACTACGCCCTGGTACGTGGT  >  NC_002947/2898011‑2898307
                                                                                                                                                   |                                                                                                                                                       
aTGCCCGGCCTGTCCCATTGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTg                                                                                                                                                       >  2:1492519/1‑150 (MQ=255)
   accGGCCTGTCCCATTGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCGCGGCCCTGAGCGAACTggca                                                                                                                                                    <  2:1372999/149‑3 (MQ=255)
      ggCCTGTCCCATTGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGgaagaa                                                                                                                                                 >  2:1710345/1‑150 (MQ=255)
          tGTCCCATTGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAacca                                                                                                                                             <  2:658250/150‑1 (MQ=255)
           gTCCCATTGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCCGAGCGAACTGGAAGAAaccac                                                                                                                                            <  1:955413/150‑1 (MQ=255)
            tCCCATTGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAaccacc                                                                                                                                           <  2:942890/150‑1 (MQ=255)
             cccATTGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCAccc                                                                                                                                          <  1:1282700/150‑1 (MQ=255)
              ccATTGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCt                                                                                                                                         >  1:161617/1‑150 (MQ=255)
               cATTGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGgaa                                                                                                                                                    <  2:1844478/138‑1 (MQ=255)
               cATTGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGgaa                                                                                                                                                    >  1:1844478/1‑138 (MQ=255)
                 ttGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGa                                                                                                                                      <  2:133805/150‑1 (MQ=255)
                   ggCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACt                                                                                                                                    <  1:1899562/150‑1 (MQ=255)
                     cagcaCCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCGTGGGTGAACTGGAAGAAACCACCCTCGACTgg                                                                                                                                  <  1:989473/150‑1 (MQ=255)
                           ccGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGAGTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTgcgc                                                                                                                            <  2:562452/150‑1 (MQ=255)
                                    tATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGGCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACt                                                                                                                                                        >  2:1311710/1‑113 (MQ=255)
                                    tATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACt                                                                                                                                                        <  1:1311710/113‑1 (MQ=255)
                                       ggCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGc                                                                                                                <  1:1910534/150‑1 (MQ=255)
                                        gCTCTTCCCCTACAAACGGCTCCCTGTCTTAGGTGCTGGGGGACTTACTCCGCCCTGGCCTGGTCGTGCTGTGCCTGCCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGcc                                                                                                               <  1:1716192/150‑1 (MQ=255)
                                         cTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACATTTAGATTTTGGCTGCGCCAGTTGCTTGGCCt                                                                                                              <  1:1297031/150‑1 (MQ=255)
                                          tATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACGACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTg                                                                                                             >  1:773098/1‑150 (MQ=255)
                                          tATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTg                                                                                                             >  1:261693/1‑150 (MQ=255)
                                              tCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTg                                                                                                         <  1:1004233/150‑1 (MQ=255)
                                                 cTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGTTCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGcc                                                                                                      >  1:497441/1‑150 (MQ=255)
                                                 cTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAATTGGAAGAAACCAccc                                                                                                                                          >  1:1658820/1‑114 (MQ=255)
                                                 cTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAATTGGAAGAAACCAccc                                                                                                                                          >  1:1658824/1‑114 (MQ=255)
                                                 cTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAATTGGAAGAAACCAccc                                                                                                                                          <  2:1658824/114‑1 (MQ=255)
                                                  ttCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGgaa                                                                                                                                                    >  2:1151731/1‑103 (MQ=255)
                                                  ttCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGgaa                                                                                                                                                    <  1:1151731/103‑1 (MQ=255)
                                                  ttCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAg                                                                                                                         >  1:1667061/1‑130 (MQ=255)
                                                  ttCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAg                                                                                                                         <  2:1667061/130‑1 (MQ=255)
                                                  ttCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCa                                                                                                     <  2:161617/150‑1 (MQ=255)
                                                  ttCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCa                                                                                                     >  2:302426/1‑150 (MQ=255)
                                                  ttCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCa                                                                                                     >  1:1248277/1‑150 (MQ=255)
                                                     cAATTGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTCCCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCagtg                                                                                                  <  1:123138/150‑1 (MQ=255)
                                                     cAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCagtg                                                                                                  <  2:261693/150‑1 (MQ=255)
                                                       aTGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGgaaga                                                                                                                                                  >  2:1008350/1‑100 (MQ=255)
                                                       aTGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGgaaga                                                                                                                                                  <  1:1008350/100‑1 (MQ=255)
                                                       aTGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCagtgga                                                                                                <  1:1850901/150‑1 (MQ=255)
                                                       aTGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCagtgga                                                                                                <  1:1661158/150‑1 (MQ=255)
                                                         ggCCCCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGGTGCTTGGCCAGACAGGCCaagcaact                                                                                              <  2:568379/150‑8 (MQ=255)
                                                            aCCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTggg                                                                                           <  1:908474/150‑1 (MQ=255)
                                                            aCCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTggg                                                                                           >  1:1564557/1‑150 (MQ=255)
                                                            aCCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAAGCGAGCCCTGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTggg                                                                                           <  1:1277446/150‑1 (MQ=255)
                                                              ccTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGt                                                                                         >  2:1792945/1‑150 (MQ=255)
                                                              ccTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGt                                                                                         >  2:1700847/1‑150 (MQ=255)
                                                                tGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAGGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGa                                                                                       <  2:1466058/150‑1 (MQ=255)
                                                                tGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGa                                                                                       >  2:150020/1‑150 (MQ=255)
                                                                 gTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGAt                                                                                      <  1:1145953/150‑1 (MQ=255)
                                                                     tCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAg                                                                                  <  1:1089665/150‑1 (MQ=255)
                                                                     tCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAg                                                                                  <  2:1868616/150‑1 (MQ=255)
                                                                     tCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAg                                                                                  <  1:1529919/150‑1 (MQ=255)
                                                                     tCGGGGCTGGGGGTCTGCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAg                                                                                  >  2:1554979/1‑150 (MQ=255)
                                                                      cGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAgg                                                                                 <  1:303342/150‑1 (MQ=255)
                                                                         tGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGaca                                                                              <  2:165965/150‑1 (MQ=255)
                                                                         tGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCCGACGATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCAGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGaca                                                                              >  1:1789750/1‑150 (MQ=255)
                                                                         ggCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGaca                                                                              <  2:357248/149‑1 (MQ=255)
                                                                               gggACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCt                                                                        >  1:578488/1‑150 (MQ=255)
                                                                                   cTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCGCCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGAc                                                                    <  1:1073780/150‑1 (MQ=255)
                                                                                    ttCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGAcc                                                                   >  1:1636058/1‑150 (MQ=255)
                                                                                    ttCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGAcc                                                                   >  1:1158908/1‑150 (MQ=255)
                                                                                    ttCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCAGTCCCATAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGAcc                                                                   >  1:612015/1‑150 (MQ=255)
                                                                                      ccTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAg                                                                 >  2:1009472/1‑150 (MQ=255)
                                                                                       cTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAATTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTggg                                                                                           <  2:1136734/123‑1 (MQ=255)
                                                                                       cTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAATTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTggg                                                                                           >  1:1136734/1‑123 (MQ=255)
                                                                                       cTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGc                                                                >  2:1232972/1‑150 (MQ=255)
                                                                                           gTACCCGCCGGGGTGTGCTGTGCGTGGCCTGGCAAGCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCAccc                                                            <  2:346784/150‑1 (MQ=255)
                                                                                                ggTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCtggtg                                                       <  2:1636058/150‑1 (MQ=255)
                                                                                                    tGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATc                                                                                     <  2:651636/116‑1 (MQ=255)
                                                                                                    tGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATc                                                                                     >  1:651636/1‑116 (MQ=255)
                                                                                                      ggCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCACGACCAGCCCCCTGGTGGCGCTg                                                 >  2:182095/1‑150 (MQ=255)
                                                                                                       gCGTACTGGGCCTCTGCTTTCAAGCCTGGCCTGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGa                                                <  2:873170/150‑1 (MQ=255)
                                                                                                          tGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCa                                             <  2:1539424/150‑1 (MQ=255)
                                                                                                            cTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGt                                           <  1:1263025/150‑1 (MQ=255)
                                                                                                             tGGGCCTGTCCTGTCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAac                                                                                                                                               >  2:1352092/1‑49 (MQ=255)
                                                                                                             tGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCa                                                                                                                          >  1:1513156/1‑70 (MQ=255)
                                                                                                             tGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTg                                          >  2:344229/1‑150 (MQ=255)
                                                                                                              gggCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGc                                         >  1:1151582/1‑150 (MQ=255)
                                                                                                                gCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGccc                                       >  2:619111/1‑150 (MQ=255)
                                                                                                                gCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGccc                                       <  2:541310/150‑1 (MQ=255)
                                                                                                                  cTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑t                                   <  1:1710345/150‑1 (MQ=255)
                                                                                                                  cTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑t                                   <  1:1492519/150‑1 (MQ=255)
                                                                                                                  cTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑t                                   >  1:1192438/1‑150 (MQ=255)
                                                                                                                   tGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGgcgc                                                   >  1:984389/1‑135 (MQ=255)
                                                                                                                   tGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTCTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGgcgc                                                   <  2:984389/135‑1 (MQ=255)
                                                                                                                     tCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGaa                                >  1:812409/1‑150 (MQ=255)
                                                                                                                      ccTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAAc                               <  1:1009472/150‑1 (MQ=255)
                                                                                                                        tGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAAcgc                             >  1:821802/1‑150 (MQ=255)
                                                                                                                         ggCAATCCAGCCCGGCGCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATACGTGCCCG‑‑AGAAcgcg                            >  2:973320/1‑150 (MQ=255)
                                                                                                                          gCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAACgcgc                           <  1:973320/150‑1 (MQ=255)
                                                                                                                                cAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAg                                                                                  <  2:1272648/91‑1 (MQ=255)
                                                                                                                                cAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAACGCGCCACTGa                     <  1:619111/150‑1 (MQ=255)
                                                                                                                                cAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAACGCGCCACTGa                     >  1:301312/1‑150 (MQ=255)
                                                                                                                                cAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAACGCGCCACTGa                     <  1:1415059/150‑1 (MQ=255)
                                                                                                                                         ccTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCCCCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAACGCGCCACTGACTACGCCCt            <  2:845503/150‑1 (MQ=255)
                                                                                                                                               cGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAACGCGCCa                         <  1:777686/131‑1 (MQ=255)
                                                                                                                                               cGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAACGCGCCa                         >  2:777686/1‑131 (MQ=255)
                                                                                                                                                gAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAACGCGCCACTGACTACGCCCTGGTACgt     <  2:1158908/150‑1 (MQ=255)
                                                                                                                                                 aaCTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGCTCCTCGCCACTGCCTCCACCCGGACCCACGAACGGAGCCAGGGTCCGCATGAACGCTCCACTGACTACGCCCGGGTACgtg    <  2:497441/150‑1 (MQ=255)
                                                                                                                                                 aaCTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAACGCGCCACTGACTACGCCCTGGTACgtg    <  1:1344354/150‑1 (MQ=255)
                                                                                                                                                   cTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAACGCGCCACTGACTACGCCCTGGTACgtggt  <  1:338166/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
ATGCCCGGCCTGTCCCATTGGCAGCACCCGGACTTCTATGGCTATTTCCCTTCCAATGGCACCCTGTCCTCGGTGCTGGGGGACTTCCTCAGTACCGGTCTGGGCGTGCTGGGCCTGTCCTGGCAATCCAGCCCGGCCCTGAGCGAACTGGAAGAAACCACCCTCGACTGGCTGCGCCAGTTGCTTGGCCTGTCTGGCCAGTGGAGTGGGGTGATCCAGGACACTGCCTCGACCAGCACCCTGGTGGCGCTGATCAGTGCCCG‑‑TGAACGCGCCACTGACTACGCCCTGGTACGTGGT  >  NC_002947/2898011‑2898307

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: