Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,459,989 A→T 5.3% L139Q (CTG→CAG)  PP_3075 ← transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,459,9890AT 5.3% 294.0 / 5.6 95L139Q (CTG→CAG) PP_3075transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (39/51);  new base T (3/2);  total (42/53)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

TTCCAGGCTGGGCGCGCTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGCCGGGG  >  NC_002947/3459841‑3460135
                                                                                                                                                    |                                                                                                                                                  
ttCCAGGCTGGGCGCGCTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAg                                                                                                                                                   >  2:374870/1‑150 (MQ=255)
ttCCAGGCTGGGCGCGCTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAg                                                                                                                                                   >  1:1491118/1‑150 (MQ=255)
 tCCAGGCTGGGCGCGCTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAgc                                                                                                                                                  <  2:465323/150‑1 (MQ=255)
 tCCAGGCTGGGCGCGCTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAgc                                                                                                                                                  <  1:784471/150‑1 (MQ=255)
   cAGGCTGGGCGCGCTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAgcgc                                                                                                                                                >  2:296624/1‑150 (MQ=255)
    aGGCTGGGCGCGCTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCTGTTCGACGAACTGCGCa                                                                                                                                               <  1:809627/150‑1 (MQ=255)
      gCTGGGCGCGCTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCaaa                                                                                                                                             >  1:472740/1‑150 (MQ=255)
        tGGGCGCGCTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGc                                                                                                                                           <  2:53361/150‑1 (MQ=255)
           gcgcgcTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCaa                                                                                                                                        <  2:227121/150‑1 (MQ=255)
           gcgcgcTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCaa                                                                                                                                        >  2:2375456/1‑150 (MQ=255)
             gcgcTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATg                                                                                                                                      >  1:1258951/1‑150 (MQ=255)
             gcgcTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATg                                                                                                                                      <  1:851102/150‑1 (MQ=255)
             gcgcTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATg                                                                                                                                      <  1:1856185/150‑1 (MQ=255)
              cgcTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATgg                                                                                                                                     <  2:310024/150‑1 (MQ=255)
                 tgctgcCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGAAGTACGTCGAACAGCGCAAAGCCAATGgcgc                                                                                                                                  >  1:1183399/1‑150 (MQ=255)
                    tgcCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGa                                                                                                                               >  2:1979403/1‑150 (MQ=255)
                     gcCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGAt                                                                                                                              >  1:2145104/1‑150 (MQ=255)
                                 ggCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTg                                                                                                                  <  2:1081730/150‑1 (MQ=255)
                                 ggCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTg                                                                                                                  >  2:2278397/1‑150 (MQ=255)
                                 ggCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTg                                                                                                                  <  2:1523350/150‑1 (MQ=255)
                                      aGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCGTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCCTCGAACAGCGCAACGCCAATGGCGCCGATCAGCGCGCCCTGGTCAt                                                                                                             <  2:571010/150‑1 (MQ=255)
                                        cTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCg                                                                                                           >  2:1532712/1‑150 (MQ=255)
                                             tttGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGAt                                                                                                      <  1:950320/150‑1 (MQ=255)
                                             tttGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGAt                                                                                                      <  1:132280/150‑1 (MQ=255)
                                               tGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATgg                                                                                                    <  1:332751/150‑1 (MQ=255)
                                                   cTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAg                                                                                                >  1:779780/1‑150 (MQ=255)
                                                         ggcgCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCAt                                                                                          <  1:1047913/150‑1 (MQ=255)
                                                         ggcgCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCAt                                                                                          <  1:326629/150‑1 (MQ=255)
                                                           cgcgcAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGa                                                                                        <  2:782920/150‑1 (MQ=255)
                                                           cgcgcAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGa                                                                                        >  1:1788179/1‑150 (MQ=255)
                                                           cgcgcAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCAGGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAGTCCACTCGCGTCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGa                                                                                        <  2:1046790/150‑1 (MQ=255)
                                                            gcgcAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGac                                                                                       <  1:51730/150‑1 (MQ=255)
                                                             cgcAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGacc                                                                                      <  1:434092/150‑1 (MQ=255)
                                                               taGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGaccac                                                                                    <  1:2061845/149‑1 (MQ=255)
                                                                aGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGAt                                                                                                      <  2:1627173/131‑1 (MQ=255)
                                                                aGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGAt                                                                                                      >  1:1627173/1‑131 (MQ=255)
                                                                aGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGaccacc                                                                                   >  1:2320275/1‑150 (MQ=255)
                                                                  ttGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTAACGAACACGCAACACACAGCGCTGAGAGCGGAAGAGCGCAGTGCCGGGGGCGCGTGTGGGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCaa                                                                                 >  2:1463524/1‑150 (MQ=255)
                                                                   tGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGAGACGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAAc                                                                                <  2:690726/150‑1 (MQ=255)
                                                                    gCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACg                                                                               <  2:779780/150‑1 (MQ=255)
                                                                      gAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACGACCAACGGc                                                                             >  1:363392/1‑150 (MQ=255)
                                                                        gCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTc                                                                           <  2:472740/150‑1 (MQ=255)
                                                                            gTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTc                                                                       <  1:316462/150‑1 (MQ=255)
                                                                                aGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTgg                                                                                                        >  2:469450/1‑113 (MQ=255)
                                                                                aGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTgg                                                                                                        <  1:469450/113‑1 (MQ=255)
                                                                                 ggCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCgt                                                                  <  2:2145104/150‑1 (MQ=255)
                                                                                  gCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCgtg                                                                 >  2:779184/1‑150 (MQ=255)
                                                                                  gCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCgtg                                                                 >  2:179348/1‑150 (MQ=255)
                                                                                   ccAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCgtgt                                                                <  2:112351/150‑1 (MQ=255)
                                                                                   ccAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCgtgt                                                                <  2:803958/150‑1 (MQ=255)
                                                                                    cAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGt                                                                                                          >  1:1549180/1‑107 (MQ=255)
                                                                                    cAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGt                                                                                                          <  2:1549180/107‑1 (MQ=255)
                                                                                        tCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGc                                                           <  2:680532/150‑1 (MQ=255)
                                                                                        tCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGc                                                           <  2:2245524/150‑1 (MQ=255)
                                                                                             ctgcATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCCTGACCACCACCGCCTCGTGCGGCGTGTGCAGCATGTc                                                      >  2:380757/1‑150 (MQ=255)
                                                                                             ctgcATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTc                                                      <  1:179348/150‑1 (MQ=255)
                                                                                             ctgcATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTc                                                      <  2:1762455/150‑1 (MQ=255)
                                                                                               gcATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCAt                                                         <  1:245708/145‑1 (MQ=255)
                                                                                               gcATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCAt                                                         >  2:245708/1‑145 (MQ=255)
                                                                                               gcATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCca                                                    <  1:713417/150‑1 (MQ=255)
                                                                                                cATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCcag                                                   >  1:159406/1‑150 (MQ=255)
                                                                                                cATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCcag                                                   >  1:210711/1‑150 (MQ=255)
                                                                                                  tGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCATCGTGATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCcagca                                                 <  2:1786357/150‑1 (MQ=255)
                                                                                                     tGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGa                                              >  2:877476/1‑150 (MQ=255)
                                                                                                     tGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGa                                              >  1:1809626/1‑150 (MQ=255)
                                                                                                     tGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCACTGGTCCTCGTGGATGGACTGACGCAGGACCACCAGCGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGa                                              <  1:647644/150‑1 (MQ=255)
                                                                                                      ggCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCaa                                                                                 >  1:2151176/1‑114 (MQ=255)
                                                                                                      ggCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCaa                                                                                 <  2:2151176/114‑1 (MQ=255)
                                                                                                      ggCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGAt                                             >  1:1490334/1‑150 (MQ=255)
                                                                                                      ggCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGAt                                             <  2:2320275/150‑1 (MQ=255)
                                                                                                        cGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTTGAGCATGTCCAGCAGGATgg                                           <  1:2353656/150‑1 (MQ=255)
                                                                                                         gTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATggg                                          <  2:159406/150‑1 (MQ=255)
                                                                                                          tAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATgg                                                                                                    >  1:1943736/1‑91 (MQ=255)
                                                                                                          tAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATgg                                                                                                    <  2:1943736/91‑1 (MQ=255)
                                                                                                              gcTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGt                                                                                                          >  1:2034777/1‑81 (MQ=255)
                                                                                                              gcTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGt                                                                                                          <  2:2034777/81‑1 (MQ=255)
                                                                                                                  gAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGggccgg                                      >  1:1966419/1‑145 (MQ=255)
                                                                                                                  gAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGggccgg                                      <  2:1966419/145‑1 (MQ=255)
                                                                                                                        cggcgACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTcac                           <  1:2260006/150‑1 (MQ=255)
                                                                                                                        cggcgACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTcac                           <  1:912521/150‑1 (MQ=255)
                                                                                                                                 tttGCGCTGTTCGACGAACTGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGa                                                             <  1:113480/105‑1 (MQ=255)
                                                                                                                                 tttGCGCTGTTCGACGAACTGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGa                                                             >  2:113480/3‑107 (MQ=255)
                                                                                                                                 tttGCGCTGTTCGACGAACTGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATcgcg                  >  2:410988/3‑150 (MQ=255)
                                                                                                                                 tttGCGCTGTTCGACGAACTGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATcgcg                  >  1:2273357/3‑150 (MQ=255)
                                                                                                                                 gCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATcgcg                  >  2:781010/1‑150 (MQ=255)
                                                                                                                                 gCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATcgcg                  >  2:195351/1‑150 (MQ=255)
                                                                                                                                   tGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCa                >  1:2210977/1‑150 (MQ=255)
                                                                                                                                     cgcAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATc              <  2:1491118/150‑1 (MQ=255)
                                                                                                                                       cAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAg            <  2:1333123/150‑1 (MQ=255)
                                                                                                                                          ttCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCt         <  2:1788179/150‑1 (MQ=255)
                                                                                                                                             gtcgAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGcc      >  2:1147359/1‑150 (MQ=255)
                                                                                                                                             gtcgAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGcc      >  2:1747694/1‑150 (MQ=255)
                                                                                                                                             gtcgAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGcc      >  2:1814099/1‑150 (MQ=255)
                                                                                                                                             gtcgAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGcc      >  1:1247494/1‑150 (MQ=255)
                                                                                                                                                 aaCAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGCCgggg  <  1:1601500/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCCAGGCTGGGCGCGCTGCTGCCGACGAACTGGGCGAAGCTGTATTTGGCCTGGCGGGCGCGCAGTTGCGAGCGGGTCGAGGCCAGTTCCTGCTGCATGCTGGCGTAGCGCTTGAGCAGCGGCGACAGGCTGCGCAGTTCGTCGAACAGCGCAAAGCCAATGGCGCCGATCAGCGCGCCCTGGTCATCGTGGATGGGCAGGCGCATGACCACCAACGGCTCGTTCGGCGTGTCGAGCATGTCCAGCAGGATGGGCCGGCCATTGCTCACCACATCGCGCATCAGGCTGCCGGGG  >  NC_002947/3459841‑3460135

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