Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 607,608 C→A 7.5% P170Q (CCG→CAG)  ribA → GTP cyclohydrolase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947607,6080CA 7.5% 288.9 / 7.4 96P170Q (CCG→CAG) ribAGTP cyclohydrolase II
Reads supporting (aligned to +/- strand):  ref base C (46/42);  new base A (6/2);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCTGGTG  >  NC_002947/607460‑607757
                                                                                                                                                    |                                                                                                                                                     
cTGGGCTGCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCGGACGACCAACAATCCGTGCAATTGCAAATGTCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCg                                                                                                                                                      <  1:659868/150‑1 (MQ=255)
  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                                                                                                                                   <  2:313405/131‑1 (MQ=255)
                                      cctgccGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAAcc                                                                                                                                   >  1:313405/1‑131 (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                                                                                                 <  2:1415629/140‑1 (MQ=255)
                                                               gTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCgg                                                                                                 >  1:1415629/1‑140 (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                                                                              <  2:1610498/137‑1 (MQ=255)
                                                                                     ccAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGCGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTg                                                                              >  1:1610498/1‑137 (MQ=255)
                                                                                      caacaaTCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGc                                                                                                                                                     >  2:118379/1‑65 (MQ=255)
                                                                                      caacaaTCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGc                                                                                                                                                     <  1:118379/65‑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                                                                        >  2:1207002/1‑136 (MQ=255)
                                                                                            tCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAAc                                                                        <  1:1207002/136‑1 (MQ=255)
                                                                                                 gcAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGc                                                     >  1:187068/1‑150 (MQ=255)
                                                                                                      gTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGcta                                                                                                                      <  1:661807/80‑1 (MQ=255)
                                                                                                      gTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGcta                                                                                                                      >  2:661807/1‑80 (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                                                                                              >  1:2029255/1‑96 (MQ=255)
                                                                                                              ccTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCaa                                                                                              <  2:2029255/96‑1 (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                                                                                >  1:15608/1‑99 (MQ=255)
                                                                                                                         tGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACTTGc                                                                                <  2:15608/99‑1 (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)
                                                                                                                                                    |                                                                                                                                                     
CTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCTGGTG  >  NC_002947/607460‑607757

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: