Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 607,611 T→G 6.2% L171R (CTG→CGG)  ribA → GTP cyclohydrolase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947607,6110TG 6.2% 278.7 / 8.1 97L171R (CTG→CGG) ribAGTP cyclohydrolase II
Reads supporting (aligned to +/- strand):  ref base T (50/41);  new base G (4/2);  total (54/43)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

GGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCTGGTGG  >  NC_002947/607462‑607758
                                                                                                                                                     |                                                                                                                                                   
gggCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCt                                                                                                                                                     >  1:2280058/1‑150 (MQ=255)
 ggCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTg                                                                                                                                                    >  1:1899748/1‑150 (MQ=255)
    ttCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGcac                                                                                                                                                 >  2:110830/1‑150 (MQ=255)
       gCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACAcc                                                                                                                                              <  1:1300009/150‑1 (MQ=255)
        cTGCCGACCAGAGTGACTACGAAATATGACTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCg                                                                                                                                             <  2:1831951/150‑1 (MQ=255)
          gCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTAATGCCCAACACTCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCggg                                                                                                                                           <  2:2267193/150‑1 (MQ=255)
            cGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCa                                                                                                                                         <  1:1326676/150‑1 (MQ=255)
                cAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAAc                                                                                                                                     >  1:831115/1‑150 (MQ=255)
                   cGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCg                                                                                                                                  <  1:1055639/150‑1 (MQ=255)
                     cttggtcGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGca                                                                                                                                >  2:451343/7‑150 (MQ=255)
                        cTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACaa                                                                                                                             <  1:110830/150‑1 (MQ=255)
                         tACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAAc                                                                                                                            >  2:2146633/1‑150 (MQ=255)
                         tACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAAc                                                                                                                            <  2:2280058/150‑1 (MQ=255)
                         tACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAAAATCGTGGTTGCCGAGCGTGTGCAGCTGCACACCGGGCATAACCCGCACAAc                                                                                                                            >  1:89187/1‑150 (MQ=255)
                         tACGCCATATACCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAAc                                                                                                                            <  1:1063288/150‑1 (MQ=255)
                         tACGCCATATACCTGCCGATCCTCAAGCACCAGGACGTCAAGTCGCTGCTCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGATGCCGAGCATGTGCCGCTGCACACCGGGCATAACCCGCACAAc                                                                                                                            >  2:670010/1‑150 (MQ=255)
                            gCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGc                                                                                                                         <  2:1503266/150‑1 (MQ=255)
                              cATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGcta                                                                                                                       <  1:2113364/150‑1 (MQ=255)
                                   gcctgccGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCGGATGACCAACAATAACAGCAAGGTCAAAGAATTGACCGACGGGAAGATCGTGGATGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTAcc                                                                                                                  >  1:2098635/1‑150 (MQ=255)
                                    cctgccGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAAcc                                                                                                                                    >  1:313405/1‑131 (MQ=255)
                                    cctgccGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAAcc                                                                                                                                    <  2:313405/131‑1 (MQ=255)
                                          gATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCggg                                                                                                                                           <  2:2098658/118‑1 (MQ=255)
                                          gATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGAACGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCggg                                                                                                                                           >  1:2098658/1‑118 (MQ=255)
                                              cTCGAGCACCTGGGCGTCAAGTCTATGCGCCTGATGATCAACAATCCGCTCAAGGGCAAAGCCCTGACCAACATGAACATCCTGGTTGCCGAGCGGGTGACGCTGCACCCCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAg                                                                                                       <  1:451343/150‑1 (MQ=255)
                                              cTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAg                                                                                                       >  1:1744236/1‑150 (MQ=255)
                                               tCGAGCACCTGGGCGTCAAGTCGCTGCTCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAgg                                                                                                      >  2:790035/1‑150 (MQ=255)
                                                  aGCACCTGGGCGTCAAGTCGATGCGCCTGGTGAGAAATAAGCCGTGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTTCCTGGCGACCAAGGCCg                                                                                                   <  2:1230833/150‑1 (MQ=255)
                                                   gCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCgg                                                                                                  >  2:1115296/1‑150 (MQ=255)
                                                   gCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCgg                                                                                                  >  1:1298067/1‑150 (MQ=255)
                                                    cACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGc                                                                                                 >  1:984685/1‑150 (MQ=255)
                                                       cTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCAGCGCAAGGTCAAAGCCCTGACCGACATGAAAATCGTGGTTGCCGATCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGGAAg                                                                                              >  1:971165/1‑150 (MQ=255)
                                                         gggCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTCCCGCTACATAGCGAGCATAAGACGGAGAACCGCTACTACCTGGCGACAAAGGCCGGCAAGtt                                                                                            >  1:24142/1‑150 (MQ=255)
                                                           gCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTgg                                                                                          <  1:972067/150‑1 (MQ=255)
                                                           gCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTgg                                                                                          <  1:230169/150‑1 (MQ=255)
                                                            cGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTggg                                                                                         >  2:88409/1‑150 (MQ=255)
                                                             gTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCgg                                                                                                  >  1:1415629/1‑140 (MQ=255)
                                                             gTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCgg                                                                                                  <  2:1415629/140‑1 (MQ=255)
                                                               cAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGAAGCTGCAAACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCca                                                                                      >  1:1489576/1‑150 (MQ=255)
                                                                aaGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCcac                                                                                     >  1:1511546/1‑150 (MQ=255)
                                                                aaGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCcac                                                                                     >  1:1608831/1‑150 (MQ=255)
                                                                  gTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCCTAGACCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACAt                                                                                   >  1:1352622/1‑150 (MQ=255)
                                                                       tGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTgg                                                                              <  2:1511546/150‑1 (MQ=255)
                                                                       tGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTgg                                                                              <  2:1668992/150‑1 (MQ=255)
                                                                        gcgcCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTggg                                                                             <  1:258790/150‑1 (MQ=255)
                                                                          gcCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCa                                                                           >  1:2400835/1‑150 (MQ=255)
                                                                            cTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAAc                                                                         <  2:831115/150‑1 (MQ=255)
                                                                             tgatgaCCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACg                                                                        >  2:2116896/1‑150 (MQ=255)
                                                                               aagaCCAACAATCCCCGCCAGCTCAAAGCCCTGCCCGACATCAACATCGTGGTTGCCGAGCGTGAGCCGCTGCACAACGCGCATAACCCGCACAACCGCTACTACCTGGCGAACAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGaa                                                                      <  1:790035/148‑1 (MQ=255)
                                                                                   ccAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGCGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTg                                                                               >  1:1610498/1‑137 (MQ=255)
                                                                                   ccAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGCGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTg                                                                               <  2:1610498/137‑1 (MQ=255)
                                                                                    caacaaTCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCa                                                                 <  1:2146633/150‑1 (MQ=255)
                                                                                    caacaaTCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCa                                                                 >  2:1154466/1‑150 (MQ=255)
                                                                                          tCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTg                                                                               >  1:1321021/1‑130 (MQ=255)
                                                                                          tCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTg                                                                               <  2:1321021/130‑1 (MQ=255)
                                                                                          tCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAAc                                                                         <  1:1207002/136‑1 (MQ=255)
                                                                                          tCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAAc                                                                         >  2:1207002/1‑136 (MQ=255)
                                                                                               gcAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGc                                                      >  1:187068/1‑150 (MQ=255)
                                                                                                    gTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGcta                                                                                                                       >  2:661807/1‑80 (MQ=255)
                                                                                                    gTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGcta                                                                                                                       <  1:661807/80‑1 (MQ=255)
                                                                                                       aaaGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCCGGCg                                              >  2:622139/1‑150 (MQ=255)
                                                                                                        aaGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGt                                             <  2:1489576/150‑1 (MQ=255)
                                                                                                        aaGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGCGCATAACTCGCACAAGCACGACGACCGGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGACGt                                             >  1:900782/1‑150 (MQ=255)
                                                                                                         aGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGCGGTGCCCCAGGCGTg                                            <  2:1608831/150‑1 (MQ=255)
                                                                                                           cccTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGAc                                          <  1:2022047/150‑1 (MQ=255)
                                                                                                           cccTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGAc                                          <  1:435001/150‑1 (MQ=255)
                                                                                                            ccTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCaa                                                                                               <  2:2029255/96‑1 (MQ=255)
                                                                                                            ccTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCaa                                                                                               >  1:2029255/1‑96 (MQ=255)
                                                                                                            ccTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAg                                                          <  2:1311742/133‑1 (MQ=255)
                                                                                                            ccTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAg                                                          >  1:1311742/1‑133 (MQ=255)
                                                                                                            ccTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGAcc                                         >  1:1791595/1‑150 (MQ=255)
                                                                                                                aCCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCcgcg                                     >  2:1840836/1‑150 (MQ=255)
                                                                                                                  cGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGcc                                   >  1:663045/1‑150 (MQ=255)
                                                                                                                       tGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACTTGc                                                                                 <  2:15608/99‑1 (MQ=255)
                                                                                                                       tGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACTTGc                                                                                 >  1:15608/1‑99 (MQ=255)
                                                                                                                           cATCGTGGTTGCCGAGCGTGTACCGCTGCACACCGGGCATAACCCGTACAACCGCTACTACCTGGCGACCAAGGCCGGAAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAg                          <  1:670010/150‑1 (MQ=255)
                                                                                                                             tCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAgcg                        >  1:1711631/1‑150 (MQ=255)
                                                                                                                              gttcgTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAac                                                                     <  1:2117733/100‑1 (MQ=255)
                                                                                                                              gttcgTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAac                                                                     >  2:2117733/5‑104 (MQ=255)
                                                                                                                               gTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAgcggc                      >  1:2028537/1‑150 (MQ=255)
                                                                                                                                tGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAg                                                          >  1:94792/1‑113 (MQ=255)
                                                                                                                                tGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAg                                                          <  2:94792/113‑1 (MQ=255)
                                                                                                                                       cGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGgcg              >  2:1615112/1‑150 (MQ=255)
                                                                                                                                       cGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGCGCCCCAGGCGTGACCCGCGCCGAGGTCAAgcgg                       >  2:1651370/1‑141 (MQ=255)
                                                                                                                                       cGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGCGCCCCAGGCGTGACCCGCGCCGAGGTCAAgcgg                       <  1:1651370/141‑1 (MQ=255)
                                                                                                                                         ccggtgtGCAGCGGCACACGGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCt            >  1:2072749/4‑150 (MQ=255)
                                                                                                                                         ccggtgtGCAGCGGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAg                                <  1:1471601/127‑1 (MQ=255)
                                                                                                                                         ccggtgtGCAGCGGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAg                                >  2:1471601/4‑130 (MQ=255)
                                                                                                                                         ccggtgtGCAGCGGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCt            >  2:551532/4‑150 (MQ=255)
                                                                                                                                            gtgtGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGCCCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGc         <  2:1352622/150‑1 (MQ=255)
                                                                                                                                            gtgtGCAGCGGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAgcggc                      >  1:961761/1‑137 (MQ=255)
                                                                                                                                            gtgtGCAGCGGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAgcggc                      <  2:961761/137‑1 (MQ=255)
                                                                                                                                              gtgCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCt       <  2:915477/150‑1 (MQ=255)
                                                                                                                                               tgCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCtg      >  2:1917115/1‑150 (MQ=255)
                                                                                                                                                  cGCTGCACACCGGGTATAACCCGCACAACCGCTACTACCTGTCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCtggtg   >  1:1997142/1‑150 (MQ=255)
                                                                                                                                                  cGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGTGGCGCCTGGCGCTGGCCtggtg   >  2:603471/1‑150 (MQ=255)
                                                                                                                                                  cGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCtggtg   >  2:329599/1‑150 (MQ=255)
                                                                                                                                                   gCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCtggtgg  <  1:88409/150‑1 (MQ=255)
                                                                                                                                                   acTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGATGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCtggtgg  >  2:455401/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCTGGTGG  >  NC_002947/607462‑607758

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