Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,773,922 A→T 5.9% T209S (ACA→TCA)  PP_5743 → TonB‑dependent receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,773,9220AT 5.9% 202.2 / 2.1 68T209S (ACA→TCA) PP_5743TonB‑dependent receptor protein
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base T (2/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

TGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCTTTT  >  NC_002947/3773776‑3774067
                                                                                                                                                   |                                                                                                                                                 
tGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACa                                                                                                                                                 >  4:68400/1‑149 (MQ=255)
    gCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGAGATggg                                                                                                                                             >  2:275224/1‑149 (MQ=255)
      ccAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCg                                                                                                                                           >  4:223978/1‑149 (MQ=255)
        aaGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTc                                                                                                                                         <  5:117629/149‑1 (MQ=255)
         aGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGTCATGGGcacgt                                                                                                                                        <  6:85405/149‑5 (MQ=255)
         aGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCt                                                                                                                                        >  7:77223/1‑149 (MQ=255)
           gTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTg                                                                                                                                       >  1:295997/1‑148 (MQ=255)
               ccGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATAc                                                                                                                                  <  3:110730/149‑1 (MQ=255)
                        gCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGt                                                                                                                         <  2:394984/149‑1 (MQ=255)
                            cAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCt                                                                                                                     <  6:326022/149‑1 (MQ=255)
                             aaCCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTg                                                                                                                    <  4:125827/149‑1 (MQ=255)
                             aaCCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTg                                                                                                                    <  5:194728/149‑1 (MQ=255)
                                  ccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGt                                                                                                                         >  1:57109/1‑139 (MQ=255)
                                  ccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGt                                                                                                                         <  2:57109/139‑1 (MQ=255)
                                  ccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGgtgt                                                                                                                >  5:175845/1‑148 (MQ=255)
                                   caTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTaa                                                                                                              >  6:345470/1‑149 (MQ=255)
                                   caTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTaa                                                                                                              >  7:24909/1‑149 (MQ=255)
                                   caTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTaa                                                                                                              >  7:30760/1‑149 (MQ=255)
                                   caTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTaa                                                                                                              <  8:77223/149‑1 (MQ=255)
                                      ggggAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACCCACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACgg                                                                                                           >  4:94604/1‑149 (MQ=255)
                                        ggAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGt                                                                                                          >  1:80056/1‑149 (MQ=255)
                                              gtaaCTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCg                                                                                                    >  3:121478/3‑148 (MQ=255)
                                               cAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGAGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGAAATGGGCTTCTTCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCa                                                                                                  <  8:72712/149‑1 (MQ=255)
                                                   tCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCt                                                                                              >  2:292972/1‑149 (MQ=255)
                                                   tCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCt                                                                                              >  2:113543/1‑149 (MQ=255)
                                                          cggcAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCgcggc                                                                                       <  5:363531/149‑1 (MQ=255)
                                                              aGCTTCCGCACGGCCTCCGAAC‑AAGCCAGCCTGCTTGGCCGCGAAGGCGAGTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTaa                                                                                   <  8:304218/149‑1 (MQ=255)
                                                               gCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAAt                                                                                  <  1:292972/149‑1 (MQ=255)
                                                                cTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATg                                                                                 <  5:345470/149‑1 (MQ=255)
                                                                   ccGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGc                                                                               <  7:259916/148‑1 (MQ=255)
                                                                   ccGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGc                                                                               <  3:223978/148‑1 (MQ=255)
                                                                    cGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCt                                                                              <  8:109848/148‑1 (MQ=255)
                                                                    cGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCt                                                                              <  2:295997/148‑1 (MQ=255)
                                                                     gCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGt                                                                                                          >  1:74626/1‑119 (MQ=255)
                                                                     gCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGt                                                                                                          <  2:74626/119‑1 (MQ=255)
                                                                            cACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTg                                                                      >  3:194898/1‑148 (MQ=255)
                                                                               cGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTg                                                                      >  8:179825/1‑145 (MQ=255)
                                                                                  aC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCa                                                                         >  6:373058/1‑139 (MQ=255)
                                                                                  aC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCa                                                                         <  5:373058/139‑1 (MQ=255)
                                                                                       gccagccTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGa                                                           <  3:68400/149‑1 (MQ=255)
                                                                                                aTTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTAccc                                                    <  7:179825/147‑1 (MQ=255)
                                                                                                     ccGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTcgcg                                             >  5:259694/1‑149 (MQ=255)
                                                                                                        cgAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGc                                          <  2:80056/149‑1 (MQ=255)
                                                                                                            ggCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGt                                      >  8:173788/1‑149 (MQ=255)
                                                                                                             gCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTc                                     >  6:224484/1‑149 (MQ=255)
                                                                                                                 tttCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGc                                 <  8:24909/149‑1 (MQ=255)
                                                                                                                 tttCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGc                                 <  8:30760/149‑1 (MQ=255)
                                                                                                                      aCTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTg                            >  2:339975/1‑149 (MQ=255)
                                                                                                                          ccACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTtgc                        >  7:338407/1‑149 (MQ=255)
                                                                                                                          ccACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTtgc                        >  2:358022/1‑149 (MQ=255)
                                                                                                                           cacacaCCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGTTCGCGGGCAGGTCCAGCAACTGTtgct                       >  8:151662/1‑149 (MQ=255)
                                                                                                                           cacacaCCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTtgct                       >  1:356448/1‑149 (MQ=255)
                                                                                                                              acacCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCa                    >  2:385393/1‑149 (MQ=255)
                                                                                                                                acCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACt                  >  6:234990/1‑149 (MQ=255)
                                                                                                                                acCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACt                  >  4:292733/1‑149 (MQ=255)
                                                                                                                                acCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACt                  >  4:129162/1‑149 (MQ=255)
                                                                                                                                  ccAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTc                <  8:338407/149‑1 (MQ=255)
                                                                                                                                   cAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTtgctagattgg                >  5:327964/1‑141 (MQ=255)
                                                                                                                                   cAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTtgc                        <  6:327964/140‑1 (MQ=255)
                                                                                                                                    aGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGc              >  5:384651/1‑149 (MQ=255)
                                                                                                                                         gacgCCCATGTCATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCg         >  1:175658/4‑149 (MQ=255)
                                                                                                                                          acgCCCATGTCATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCtctc          >  7:234352/3‑147 (MQ=255)
                                                                                                                                          acgCCCATGTCATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCtctc          <  8:234352/145‑1 (MQ=255)
                                                                                                                                           tGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGc       <  4:121478/149‑1 (MQ=255)
                                                                                                                                             cccATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCt     >  3:28884/1‑149 (MQ=255)
                                                                                                                                              ccatgccatgGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCCGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCtt    >  3:227293/1‑149 (MQ=255)
                                                                                                                                              ccATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCtt    >  5:13017/1‑149 (MQ=255)
                                                                                                                                                aTGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCtttt  <  6:259694/149‑1 (MQ=255)
                                                                                                                                                aTGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCtttt  <  1:385393/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAAC‑AAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCTTTT  >  NC_002947/3773776‑3774067

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