Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,899,213 Δ1 bp 6.0% coding (1395/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,2130G. 6.0% 266.8 / 10.9 83coding (1395/1413 nt)PP_2552DOPA decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (34/44);  new base . (3/2);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

ATACCAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGG  >  NC_002947/2899065‑2899361
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                                <  1:1599644/129‑1 (MQ=255)
                                                          tACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTAt                                                                                                                >  2:1599644/1‑129 (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                                                                                                                                 <  1:1912604/97‑1 (MQ=255)
                                                                         gTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGc                                                                                                                                 >  2:1912604/1‑97 (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                                            <  1:1169325/124‑1 (MQ=255)
                                                                                                                                   cTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaaa                                            >  2:1169325/1‑124 (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)
                                                                                                                                                    aGGACGTGATCAAGGGCTGAAGCCTGTACGGTACATGTATGAGCGGGCTTGCTCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCgg  >  1:673974/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATACCAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGG  >  NC_002947/2899065‑2899361

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