Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,899,207:1 +C 6.0% coding (1389/1413 nt) PP_2552 → DOPA decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,899,2071.C 6.0% 350.5 / 14.4 84coding (1389/1413 nt)PP_2552DOPA decarboxylase
Reads supporting (aligned to +/- strand):  ref base . (34/45);  new base C (3/2);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

TGCGCATACCAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCT  >  NC_002947/2899060‑2899351
                                                                                                                                                    |                                                                                                                                                
tGCGCATACCAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑Ct                                                                                                                                                <  1:557345/150‑1 (MQ=255)
     aTACCAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGGTGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAgg                                                                                                                                           >  2:1522108/1‑150 (MQ=255)
      tACCAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTNCCGACCNAGCGGGGGGATGTGCAGCGGCTGGGGGCACGT‑CTGCAGGa                                                                                                                                          <  1:1264295/150‑1 (MQ=255)
      tACCAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGa                                                                                                                                          <  1:494427/150‑1 (MQ=255)
        ccAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACg                                                                                                                                        <  2:763035/150‑1 (MQ=255)
           aGGGCTGGGCCGAGCGGCTGAATGCATCCTGCGCTTCTTATGTGACGCCGGCTACACTGGACGTGCGGTGCAGGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCCTGCA‑GACGTgc                                                                                                                                     <  1:1852589/150‑2 (MQ=255)
            gggCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCNGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGAt                                                                                                                                    >  1:1129918/1‑150 (MQ=255)
                 gggCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAgg                                                                                                                               <  1:860318/150‑1 (MQ=255)
                  ggCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAggg                                                                                                                              <  2:1948457/150‑1 (MQ=255)
                       aGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGaa                                                                                                                         <  2:1635050/150‑1 (MQ=255)
                           gCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCt                                                                                                                     <  2:2227683/150‑1 (MQ=255)
                            cTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTg                                                                                                                    >  2:2126368/1‑150 (MQ=255)
                              gAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTa                                                                                                                  <  1:1480839/150‑1 (MQ=255)
                              gAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTa                                                                                                                  >  2:2408470/1‑150 (MQ=255)
                              gAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTa                                                                                                                  <  2:2397346/150‑1 (MQ=255)
                               aaTGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTAc                                                                                                                 >  2:1376321/1‑150 (MQ=255)
                                   cATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTa                                                                                                             <  1:2126368/150‑1 (MQ=255)
                                     tCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTatcagatt                                                                                                           >  1:650011/1‑143 (MQ=255)
                                       cGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTg                                                                                                         <  2:542715/150‑1 (MQ=255)
                                       cGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTg                                                                                                         >  2:226306/1‑150 (MQ=255)
                                         gcgcTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTa                                                                                                       >  2:199827/1‑150 (MQ=255)
                                           gctgctTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTAt                                                                                                      <  2:1943564/149‑1 (MQ=255)
                                           gctgctTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTAt                                                                                                      >  1:1943564/1‑149 (MQ=255)
                                           gctgctTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGAGCTGAAGCAAGTACGGTACCTGTAtt                                                                                                     >  2:1343853/1‑149 (MQ=255)
                                            ctgctTATGTGCCGCCCGCTACACTGGACGGGCGGTGGATGCAGCGGGGTTCGATGGGTGCGCTGCCGACCGAGCGGGGGGATTTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGa                                                                                                    <  2:920354/150‑1 (MQ=255)
                                            ctgctTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGa                                                                                                    >  2:788818/1‑150 (MQ=255)
                                              gctTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGc                                                                                                  <  1:2110806/150‑1 (MQ=255)
                                              gctTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGAGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGc                                                                                                  >  1:1888188/1‑150 (MQ=255)
                                                ttATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGAc                                                                                                                                         <  1:1509447/109‑1 (MQ=255)
                                                ttATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGAc                                                                                                                                         >  2:1509447/1‑109 (MQ=255)
                                                ttATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGc                                                                                                                       >  2:805511/1‑127 (MQ=255)
                                                ttATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGc                                                                                                                       <  1:805511/127‑1 (MQ=255)
                                                ttATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCgg                                                                                                <  2:1853062/150‑1 (MQ=255)
                                                 tATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGAACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCggg                                                                                               <  1:1159291/150‑1 (MQ=255)
                                                   tgtgACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCt                                                                                             <  1:1130359/150‑1 (MQ=255)
                                                    gtgACGCCGGCTACCTTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCAGCCCGATGGTGTGGATGTCCAAAGGCTGTGGGCACGT‑CTGCAGGAAGTGATCAAGGGCTTAAGCCTGTACGGTACCTGTATGAGCGGGCtt                                                                                            <  1:1376321/150‑1 (MQ=255)
                                                       aCGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGcc                                                                                         >  2:1155537/1‑150 (MQ=255)
                                                              cTACACTGGACGGGCGGTGGATCGTGCGCGTTTTGCTTCGTGAGTTCCCCCCTGATTGGCGTGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAg                                                                                  <  2:1129918/150‑1 (MQ=255)
                                                               tACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGTGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAggg                                                                                                                              >  1:39409/1‑105 (MQ=255)
                                                               tACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGTGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAggg                                                                                                                              <  2:39409/105‑1 (MQ=255)
                                                               tACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTAt                                                                                                      >  2:1599644/1‑129 (MQ=255)
                                                               tACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTAt                                                                                                      <  1:1599644/129‑1 (MQ=255)
                                                               tACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAgg                                                                                 <  1:278355/150‑1 (MQ=255)
                                                                    tGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTg                                                                            >  2:1539188/1‑150 (MQ=255)
                                                                      gACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGACGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCa                                                                          >  1:2068219/1‑150 (MQ=255)
                                                                        cGGGCGGTGGATGGTGCGGGTTTCGATTGGTTCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCaaa                                                                        <  2:1103567/150‑1 (MQ=255)
                                                                         gggCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCaa                                                                         <  1:1529517/148‑1 (MQ=255)
                                                                         gggCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATAAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCaa                                                                         >  2:1529517/1‑148 (MQ=255)
                                                                         gggCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGACTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAag                                                                       >  1:1982068/1‑150 (MQ=255)
                                                                             ggTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagcagc                                                                   <  1:1725361/150‑1 (MQ=255)
                                                                              gTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGc                                                                                                                       >  2:1912604/1‑97 (MQ=255)
                                                                              gTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGc                                                                                                                       <  1:1912604/97‑1 (MQ=255)
                                                                                 gATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCt                                                               <  1:1539188/150‑1 (MQ=255)
                                                                                      tGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAAt                                                          <  2:440749/150‑1 (MQ=255)
                                                                                          ggTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTAt                                                      <  1:2408470/150‑1 (MQ=255)
                                                                                            tttCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCa                                                    <  2:852769/150‑1 (MQ=255)
                                                                                                gATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTgg                                                <  2:2297211/150‑1 (MQ=255)
                                                                                                     gTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagca                                                                     <  1:1877946/124‑1 (MQ=255)
                                                                                                     gTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagca                                                                     >  2:1877946/1‑124 (MQ=255)
                                                                                                     gTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCAc                                           >  2:1258288/1‑150 (MQ=255)
                                                                                                           tGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCACAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTg                                     >  1:101899/1‑150 (MQ=255)
                                                                                                              cgaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaaa                                  <  1:1155537/150‑1 (MQ=255)
                                                                                                                  cgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTgg                                                <  1:1102859/132‑1 (MQ=255)
                                                                                                                  cgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTgg                                                >  2:1102859/1‑132 (MQ=255)
                                                                                                                        gggggATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTg                                     >  2:1556258/1‑137 (MQ=255)
                                                                                                                        gggggATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTg                                     <  1:1556258/137‑1 (MQ=255)
                                                                                                                           ggATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACgg                     >  2:172725/1‑150 (MQ=255)
                                                                                                                                    gCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGc                                        >  2:918586/1‑122 (MQ=255)
                                                                                                                                    gCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGc                                        <  1:918586/122‑1 (MQ=255)
                                                                                                                                    gCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCa            <  1:1258288/150‑1 (MQ=255)
                                                                                                                                        cTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaaa                                  >  2:1169325/1‑124 (MQ=255)
                                                                                                                                        cTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaaa                                  <  1:1169325/124‑1 (MQ=255)
                                                                                                                                        cTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAgg        <  2:1416876/150‑1 (MQ=255)
                                                                                                                                         tgagggCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGTCTGGCTTCATAGGc       >  1:2401812/4‑150 (MQ=255)
                                                                                                                                            gggCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGc    <  2:1310660/150‑1 (MQ=255)
                                                                                                                                             atcACGTCCTGCA‑GACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGcc   >  2:2125168/3‑150 (MQ=255)
                                                                                                                                              tcACGTCCTGCA‑GACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagca                                                                     >  2:1633034/2‑83 (MQ=37)
                                                                                                                                              tcACGTCCTGCA‑GACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagca                                                                     <  1:1633034/82‑1 (MQ=37)
                                                                                                                                              tcACGTCCTGCA‑GACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCt  >  1:764150/2‑150 (MQ=255)
                                                                                                                                              gCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggctggc                <  1:2121782/136‑1 (MQ=255)
                                                                                                                                              gCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggctggc                >  2:2121782/1‑136 (MQ=255)
                                                                                                                                              gCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCt  >  2:894538/1‑150 (MQ=255)
                                                                                                                                                  gT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCt                               <  1:1695987/117‑1 (MQ=255)
                                                                                                                                                  gT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCt                               >  2:1695987/1‑117 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TGCGCATACCAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCT  >  NC_002947/2899060‑2899351

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