Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 837,811 T→C 6.2% K52E (AAA→GAA)  pth ← peptidyl‑tRNA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947837,8110TC 6.2% 234.0 / 4.4 80K52E (AAA→GAA) pthpeptidyl‑tRNA hydrolase
Reads supporting (aligned to +/- strand):  ref base T (33/42);  new base C (2/3);  total (35/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GGGGCAGATCGAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGA  >  NC_002947/837664‑837954
                                                                                                                                                       |                                                                                                                                               
ggggCTGATGGAGTTCGTCTTGTGCGACCAGTATCGCTTCCGGCTTTATGCGTAAGAAATTGGCCAAGGCCGCTACGGTCTGGGCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTTGGGAGCATCAAACGCTCGTCGTTGCCCTGATGGCTGAATTTa                                                                                                                                                <  3:320725/149‑1 (MQ=255)
ggggCAGATCGAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTa                                                                                                                                                <  3:274961/149‑1 (MQ=255)
ggggCAGATCGAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTa                                                                                                                                                <  6:416598/149‑1 (MQ=255)
 gggCAGATCGAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAg                                                                                                                                               >  8:366483/1‑149 (MQ=255)
   gCAGATCGAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGACGGCTAAATTCAGc                                                                                                                                              >  1:335037/1‑148 (MQ=255)
     aGATCGAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCTAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGt                                                                                                                                           >  3:5336/1‑149 (MQ=255)
     aGATCGAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTTGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGt                                                                                                                                           <  1:325402/149‑1 (MQ=255)
       aTCGAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGt                                                                                                                                           >  7:216386/1‑147 (MQ=255)
        tCGAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGACGGCTAAATTCAGCCATCAg                                                                                                                                        <  2:335037/149‑1 (MQ=255)
         cGAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCATCATAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCa                                                                                                                                         <  1:125376/148‑1 (MQ=255)
          gAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGc                                                                                                                                      <  3:414003/149‑1 (MQ=255)
          gAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGACGGCTAAATTCAGCCATCAgg                                                                                                                                       <  8:218718/148‑1 (MQ=255)
                gtcgtGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCTAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGaaa                                                                                                                                 <  4:5336/148‑1 (MQ=255)
                       gCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACttttt                                                                                                                         >  5:207525/1‑149 (MQ=255)
                        ccaccaGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAAGGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACttttt                                                                                                                         <  1:204923/148‑1 (MQ=255)
                        ccaccaGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTg                                                                                                                        <  1:373100/149‑1 (MQ=255)
                        ccaccaGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTg                                                                                                                        >  4:267460/1‑149 (MQ=255)
                        ccaccaGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTg                                                                                                                        >  4:168918/1‑149 (MQ=255)
                         caccaGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTgt                                                                                                                       <  1:205755/149‑1 (MQ=255)
                         caccaGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTg                                                                                                                        >  7:244599/1‑148 (MQ=255)
                          accaGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGTGTATCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACttttt                                                                                                                         >  1:136368/1‑149 (MQ=255)
                          accaGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTgt                                                                                                                       >  4:283916/1‑148 (MQ=255)
                          accaGTATCGCTTACGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTgt                                                                                                                       >  7:107269/1‑148 (MQ=255)
                           ccaGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTgtca                                                                                                                     <  8:253995/149‑1 (MQ=255)
                            caGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGTCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTgtcag                                                                                                                    <  5:323953/149‑1 (MQ=255)
                                 tCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGTCCGAAATACTTTTTgtcagcggt                                                                                                                <  8:107269/148‑1 (MQ=255)
                                 tCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTgtcagcggtc                                                                                                               <  3:255856/149‑1 (MQ=255)
                                 tCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTgtcagcggt                                                                                                                <  7:366483/148‑1 (MQ=255)
                                   gCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTgtcagcggtcag                                                                                                             <  4:211560/149‑1 (MQ=255)
                                    cTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTgtcagcggtcagc                                                                                                            <  2:426385/149‑1 (MQ=255)
                                     ttCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATAc                                                                                                                              >  2:68433/1‑130 (MQ=255)
                                     ttCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATAc                                                                                                                              <  1:68433/130‑1 (MQ=255)
                                        cGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGaaa                                                                                                        >  5:75930/1‑149 (MQ=255)
                                          gCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAaca                                                                                                      >  2:58560/1‑149 (MQ=255)
                                          gCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAaca                                                                                                      <  3:168918/149‑1 (MQ=255)
                                               aTGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCt                                                                                                  <  8:244599/148‑1 (MQ=255)
                                               aTGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCt                                                                                                  <  7:374156/148‑1 (MQ=255)
                                                tGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTg                                                                                                 <  3:429708/148‑1 (MQ=255)
                                                   ggAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGgcg                                                                                             >  5:401268/1‑149 (MQ=255)
                                                    gaagaaATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGgcgc                                                                                            >  8:347540/1‑149 (MQ=255)
                                                    gaagaaATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACCCTCGgcgc                                                                                            >  6:208959/1‑149 (MQ=255)
                                                      agaaATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTg                                                                                          >  4:219088/1‑149 (MQ=255)
                                                      agaaATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTg                                                                                          <  3:358308/149‑1 (MQ=255)
                                                               ccAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGt                                                                                 >  7:137116/1‑149 (MQ=255)
                                                                    gCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAg                                                                                     <  5:314547/140‑1 (MQ=255)
                                                                    gCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAg                                                                                     >  6:314547/1‑140 (MQ=255)
                                                                    gCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAAc                                                                            <  6:75930/149‑1 (MQ=255)
                                                                    gCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAAc                                                                            <  3:283916/149‑1 (MQ=255)
                                                                                tGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAg                                                                                     <  2:303995/128‑1 (MQ=255)
                                                                                tGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAg                                                                                     >  1:303995/1‑128 (MQ=255)
                                                                                   ccGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGc                                                             <  3:73188/149‑1 (MQ=255)
                                                                                     gCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGt                                                           <  8:216386/149‑1 (MQ=255)
                                                                                       tGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTa                                                         >  3:160140/1‑149 (MQ=255)
                                                                                         cGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTAt                                                        <  3:219088/148‑1 (MQ=255)
                                                                                         cGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATg                                                       >  7:216043/1‑149 (MQ=255)
                                                                                           gT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGcc                                                     >  3:160343/1‑149 (MQ=255)
                                                                                            t‑‑‑tCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCACTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCg                                                    <  6:207525/149‑1 (MQ=255)
                                                                                                 cAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCgg                                                   <  6:107082/148‑1 (MQ=255)
                                                                                                 cAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCgg                                                   <  4:160343/148‑1 (MQ=255)
                                                                                                   t‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTc                                                >  1:143951/1‑149 (MQ=255)
                                                                                                        ggtggtGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGt                                             <  7:347540/148‑1 (MQ=255)
                                                                                                         gttgtggGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTc                                           <  3:355878/146‑1 (MQ=255)
                                                                                                         gtggtgGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTc                                           <  2:136368/149‑1 (MQ=255)
                                                                                                           ggtggGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGt                                         >  7:353411/1‑149 (MQ=255)
                                                                                                               gggATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTc                                     >  2:149058/1‑149 (MQ=255)
                                                                                                                ggATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCg                                    >  3:321334/1‑149 (MQ=255)
                                                                                                                    caacaaACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGc                                >  3:201763/1‑149 (MQ=255)
                                                                                                                     aacaaACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGTCGGGTCTGTTCGTATTCGGGGcc                               >  7:157566/1‑149 (MQ=255)
                                                                                                                       caaACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCgg                             <  1:58560/149‑1 (MQ=255)
                                                                                                                        aaaCGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCggg                            >  4:348778/1‑149 (MQ=255)
                                                                                                                        aaaCGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCggg                            <  6:401268/149‑1 (MQ=255)
                                                                                                                            gaacgTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTa                        >  3:36571/1‑149 (MQ=255)
                                                                                                                                 tcgtTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAg                   <  8:216043/149‑1 (MQ=255)
                                                                                                                                       gccTGACGGCTAAATTCAGCCATCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAAc             >  5:149665/2‑149 (MQ=255)
                                                                                                                                               gCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACt     <  2:143951/149‑1 (MQ=255)
                                                                                                                                               gCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACt     <  2:337766/149‑1 (MQ=255)
                                                                                                                                                 tGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCCGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTg    >  4:307693/1‑148 (MQ=255)
                                                                                                                                                  gAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGa  >  2:295339/1‑149 (MQ=255)
                                                                                                                                                   aaTTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGa  <  5:288298/148‑1 (MQ=255)
                                                                                                                                                   aaTTCAGCCATCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGa  <  6:149665/148‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
GGGGCAGATCGAGTTCGTCGTGTGCCACCAGTATCGCTTCCGGCTTGATGCGGAAGAAATTGGCCAAGGCCGCTACGGACTGGCCGCTGCGGT‑‑‑TCAT‑GTAGGTGGTGGGGATCAACAAACGAACGTCGTTGCCCTGATGGCTGAATTTAGCCGTCAGGCCGAAATACTTTTTGTCAGCGGTCAGCGAAACACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGA  >  NC_002947/837664‑837954

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