Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 583,719 C→G 8.7% P529A (CCC→GCC) ‡ selB → selenocysteyl‑tRNA‑specific translation elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947583,7190CG 8.7% 271.5 / 14.1 92P529A (CCC→GCC) ‡selBselenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base C (43/40);  new base G (6/2);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

TGCTGGATGAACTGCTCGACGACAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCG  >  NC_002947/583573‑583865
                                                                                                                                                  |                                                                                                                                                  
tGCTGGATGAACTGCTCGACGACAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGaaccaccc                                                                                                                                                  >  1:370172/1‑149 (MQ=255)
   tGGATGAACTGCTCGACGACAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTgg                                                                                                                                               <  7:181913/149‑1 (MQ=255)
   tGGATGAACTGCTCGACGACAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTgg                                                                                                                                               >  6:146503/1‑149 (MQ=255)
     gATGAACTGCTCGACGACAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGt                                                                                                                                             >  3:412908/1‑149 (MQ=255)
      aTGAACTGCTCGACGACAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGt                                                                                                                                             >  5:389759/1‑148 (MQ=255)
       tGAACTGCTCGACGACAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCAGGGTGGGTCa                                                                                                                                           <  8:43557/149‑1 (MQ=255)
           cTGCTCGACGACAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACGCAGGGTGGGTCaaat                                                                                                                                        <  4:199772/148‑4 (MQ=255)
              cTCGACGACAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCt                                                                                                                                    >  7:239765/1‑149 (MQ=255)
                  acgacAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGc                                                                                                                                 >  1:155628/1‑148 (MQ=255)
                  acgacAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGc                                                                                                                                 >  3:351177/1‑148 (MQ=255)
                     acAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGt                                                                                                                             <  4:412908/149‑1 (MQ=255)
                      cAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTg                                                                                                                            <  1:12919/149‑1 (MQ=255)
                       aGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGt                                                                                                                             >  4:317514/1‑147 (MQ=255)
                        gcgcCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTTCAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAACCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTgag                                                                                                                          <  2:56039/149‑1 (MQ=255)
                        gcgcCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTgag                                                                                                                          <  1:205301/149‑1 (MQ=255)
                         agcCATGGCCCGCTGCTGGCGATGGCTGCATATGCCCGATCACAAGGTACAGTTGAGCGAAGACGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATAGCAGGGCAATTTGACCCACCCGGGGTCAGGCCGCTGTCGAGTgagg                                                                                                                         >  6:281126/2‑149 (MQ=255)
                             aTCGACAGCAGCGGGACATGGATGCATCTTACAGATCACAAGGTACAGTTGAGCGAAGCCTACAGTGCATTGTCGGCACGTATGAAGCCGAAGTTGATCGCAGTGCAATTTGACACACACTGGGTCATGACGATGGCGAGTGAGgaga                                                                                                                      >  4:417964/1‑148 (MQ=255)
                                 ccAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACtgt                                                                                                                 <  7:346608/149‑1 (MQ=255)
                                 ccAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCACCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACtgt                                                                                                                 <  3:184044/149‑1 (MQ=255)
                                  cagcagCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACtgtg                                                                                                                <  2:387687/149‑1 (MQ=255)
                                        cGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAg                                                                                                          >  3:359912/1‑149 (MQ=255)
                                        cGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAg                                                                                                          <  1:79086/149‑1 (MQ=255)
                                        cGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAg                                                                                                          <  3:154620/149‑1 (MQ=255)
                                         gggCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGc                                                                                                         <  8:241655/149‑1 (MQ=255)
                                           gCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGcc                                                                                                        >  2:435319/1‑148 (MQ=255)
                                           gCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCg                                                                                                       <  5:324544/149‑1 (MQ=255)
                                            ccATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGa                                                                                                      >  1:177441/1‑149 (MQ=255)
                                             cATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGAt                                                                                                     <  4:276419/149‑1 (MQ=255)
                                                ggCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGa                                                                                                      <  8:279016/145‑1 (MQ=255)
                                                ggCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGa                                                                                                      >  7:279016/1‑145 (MQ=255)
                                                 gCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTgc                                                                                                 <  1:435319/149‑1 (MQ=255)
                                                     cATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGAtgt                                                                                                   >  7:373235/1‑143 (MQ=255)
                                                     cATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGAtgt                                                                                                   <  8:373235/143‑1 (MQ=255)
                                                     cATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTgcgc                                                                                               <  8:52752/147‑1 (MQ=255)
                                                     cATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGcc                                                                                              <  3:317514/148‑1 (MQ=255)
                                                          gCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTAc                                                                                        >  8:245682/1‑149 (MQ=255)
                                                           cccGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACt                                                                                       >  2:43829/1‑149 (MQ=255)
                                                           cccGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACt                                                                                       <  2:278189/149‑1 (MQ=255)
                                                            ccGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGTGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTg                                                                                      >  7:220864/1‑149 (MQ=255)
                                                            ccGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTg                                                                                      >  4:143732/1‑149 (MQ=255)
                                                            ccGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTg                                                                                      <  6:11821/149‑1 (MQ=255)
                                                                tCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTgcgc                                                                                   <  4:359912/148‑1 (MQ=255)
                                                                      ggTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTgg                                                                            <  7:245682/149‑1 (MQ=255)
                                                                       gTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGc                                                                           <  1:43829/149‑1 (MQ=255)
                                                                       gTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGc                                                                           >  2:151779/1‑149 (MQ=255)
                                                                              tGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCt                                                                     <  2:370172/148‑1 (MQ=255)
                                                                                aGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTgg                                                                   <  2:155628/148‑1 (MQ=255)
                                                                                  cGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGcc                                                                >  8:57026/1‑149 (MQ=255)
                                                                                   gAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCt                                                               >  4:56786/1‑149 (MQ=255)
                                                                                   gAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCt                                                               >  3:77317/1‑149 (MQ=255)
                                                                                    aaGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCTGAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTg                                                              >  1:309200/1‑149 (MQ=255)
                                                                                     aGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTgg                                                             <  3:51711/149‑1 (MQ=255)
                                                                                      gCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGt                                                            <  8:239765/149‑1 (MQ=255)
                                                                                        cGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGc                                                          <  2:23568/149‑1 (MQ=255)
                                                                                         gACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCa                                                         <  4:351177/149‑1 (MQ=255)
                                                                                          aCAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCAt                                                        <  1:169497/149‑1 (MQ=255)
                                                                                           cAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCAt                                                        >  6:165516/1‑148 (MQ=255)
                                                                                              tGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAgg                                                    >  2:108422/1‑149 (MQ=255)
                                                                                                cATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAggtg                                                  <  3:10590/149‑1 (MQ=255)
                                                                                                cATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAggtg                                                  <  1:264752/149‑1 (MQ=255)
                                                                                                cATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAggt                                                   <  5:146503/148‑1 (MQ=255)
                                                                                                 aTTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAggtgg                                                 >  8:182794/1‑149 (MQ=255)
                                                                                                      gggCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTgcgc                                            >  5:83279/1‑149 (MQ=255)
                                                                                                       ggCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTgcgc                                            >  5:266941/1‑148 (MQ=255)
                                                                                                          aCGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACt                                        >  2:198684/1‑149 (MQ=255)
                                                                                                              cTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGt                                     >  3:122834/1‑148 (MQ=255)
                                                                                                               tGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTc                                   <  4:122834/149‑1 (MQ=255)
                                                                                                                gCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTc                                   >  3:218828/1‑148 (MQ=255)
                                                                                                                  aGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTAc                                <  6:389759/149‑1 (MQ=255)
                                                                                                                   gCCGAAGTTGATCGCAGGGCAATTTGACCCATCCTGGGTCAGGCCGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTAcc                               >  5:6821/1‑149 (MQ=255)
                                                                                                                    ccGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTAccc                              >  3:155248/1‑149 (MQ=255)
                                                                                                                        aGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAg                          <  7:182794/149‑1 (MQ=255)
                                                                                                                           tGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAggt                                                   >  3:268656/1‑121 (MQ=255)
                                                                                                                           tGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAggt                                                   <  4:268656/121‑1 (MQ=255)
                                                                                                                            gATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCACGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGa                      >  5:332903/1‑149 (MQ=255)
                                                                                                                              tCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACg                    >  4:353693/1‑149 (MQ=255)
                                                                                                                              tCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACg                    >  8:69845/1‑149 (MQ=255)
                                                                                                                               cGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGc                   >  8:222955/1‑149 (MQ=255)
                                                                                                                                gCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACGTGTTCTACCCCGAGGCGACGCt                  >  7:384392/1‑149 (MQ=255)
                                                                                                                                  aGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTg                 >  3:224383/1‑148 (MQ=255)
                                                                                                                                   gggCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCa               <  8:269409/149‑1 (MQ=255)
                                                                                                                                   gggAAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGGGAGGAGAACTGTGCCGAAGCCTATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCa               <  5:281126/149‑1 (MQ=255)
                                                                                                                                    ggCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCAGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAg              >  2:91610/1‑149 (MQ=255)
                                                                                                                                     gCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAgc             >  2:348754/1‑149 (MQ=255)
                                                                                                                                      cAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGAc                                         >  6:158910/1‑120 (MQ=255)
                                                                                                                                      cAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGAc                                         <  5:158910/120‑1 (MQ=255)
                                                                                                                                        tcctGACCCAGGGTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCa          >  4:5271/4‑149 (MQ=255)
                                                                                                                                        tcctGACCCAGGGTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCa          >  5:152361/4‑149 (MQ=255)
                                                                                                                                        tcctGACCCAGGGTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCa          >  8:231329/4‑149 (MQ=255)
                                                                                                                                        tcctGACCCAGGGTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCa          >  8:225217/4‑149 (MQ=255)
                                                                                                                                          ctGACCCAGGGTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGtt                                    >  1:98862/2‑121 (MQ=255)
                                                                                                                                          ctGACCCAGGGTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCagg        >  6:2548/2‑147 (MQ=255)
                                                                                                                                            gACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATgg       >  2:433737/1‑148 (MQ=255)
                                                                                                                                               ccacccTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGc   >  4:286584/1‑149 (MQ=255)
                                                                                                                                               ccacccTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGc   >  3:294867/1‑149 (MQ=255)
                                                                                                                                                cacccTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGTTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCg  >  1:113078/1‑149 (MQ=255)
                                                                                                                                                 acccTGGGTCAGGACGCTGGCGAGTGTGGAGAACTGTGGCGAAGCCTATGTGCGCCTGTTACTGCGCATGTTGGCCGGGGTGGGGCTGGTGGATCAGGTTGTGCGCGACTTGtt                                    <  7:78473/114‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGCTGGATGAACTGCTCGACGACAGCGCCATCGCCAGCAGCGGGCCATGGCTGCATCTGCCCGATCACAAGGTACAGTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCG  >  NC_002947/583573‑583865

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