Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,257,471 G→T 5.7% T161N (ACC→AAC)  PP_2854 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,257,4710GT 5.7% 269.3 / 3.3 84T161N (ACC→AAC) PP_2854hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/41);  new base T (3/2);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CCAAGCAGCGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAGGC  >  NC_002947/3257324‑3257618
                                                                                                                                                     |                                                                                                                                                   
ccAAGCAGCGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGt                                                                                                                                                    <  6:179234/149‑1 (MQ=255)
 cAAGCAGCGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTc                                                                                                                                                   >  4:334181/1‑149 (MQ=255)
         gTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggc                                                                                                                                           <  4:239149/149‑1 (MQ=255)
            gAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCgccggc                                                                                                                                        <  6:129211/149‑1 (MQ=255)
             aGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCggcggc                                                                                                                                        >  4:194154/1‑148 (MQ=255)
               ttCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAg                                                                                                                                     >  3:134150/1‑149 (MQ=255)
                gacgctgcagGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCa                                                                                                                                      <  8:2871/141‑1 (MQ=255)
                 cTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTa                                                                                                                                   >  3:149056/1‑149 (MQ=255)
                  tGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggcg                                                                                                                                          >  1:116993/1‑141 (MQ=255)
                  tGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggcg                                                                                                                                          <  2:116993/141‑1 (MQ=255)
                   gcagcaGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTAcc                                                                                                                                 >  2:304903/1‑149 (MQ=255)
                    cagcagGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCGCCGGCCAgaccct                                                                                                                                <  4:197112/149‑6 (MQ=255)
                    cagcagGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcgcc                                                                                                                                           <  2:149207/138‑3 (MQ=255)
                    cagcagGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcgcc                                                                                                                                           >  1:149207/1‑136 (MQ=255)
                       cagGCCACCTTGGGCCAGGCGTTGCACGGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTgg                                                                                                                              >  1:204947/1‑149 (MQ=255)
                          gCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCAGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGGTACTGGccg                                                                                                                           <  1:84394/149‑2 (MQ=255)
                          gCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCa                                                                                                                                      >  8:4071/1‑137 (MQ=255)
                          gCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCa                                                                                                                                      <  7:4071/137‑1 (MQ=255)
                          gCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCa                                                                                                                                      <  7:276453/137‑1 (MQ=255)
                          gCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCa                                                                                                                                      >  8:276453/1‑137 (MQ=255)
                            cACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGCAACACACACTGAGCGCGCCGGCGCAATTTTCCCCGACCAGAGTCTGGCCGGCGGCAAGTACCTGGCCTTGc                                                                                                                        >  4:254603/1‑149 (MQ=255)
                             aCCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCt                                                                                                                       >  2:35355/1‑149 (MQ=255)
                              ccTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGc                                                                                                                        <  4:238820/147‑1 (MQ=255)
                              ccTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGc                                                                                                                        >  3:238820/1‑147 (MQ=255)
                                    gCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCt                                                                                                                                                  <  6:287484/115‑1 (MQ=255)
                                      cAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTAcgg                                                                                                              >  7:154547/1‑149 (MQ=255)
                                        ggCGTTGCACGCTGCTGCCCTGGTTTGTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGCCTGGCCGGCGGCCAGTACCTGGCCTTGCt                                                                                                                       <  7:170550/138‑1 (MQ=255)
                                        ggCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTCCTGGCCGGCGGCCAGTACATGGACTTgat                                                                                                                       >  8:170550/1‑136 (MQ=255)
                                            ttGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTAc                                                                                                                                  <  7:346525/123‑1 (MQ=255)
                                            ttGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTAc                                                                                                                                  >  8:346525/1‑123 (MQ=255)
                                              gCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGcc                                                                                                      >  6:296968/1‑149 (MQ=255)
                                              gCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGc                                                                                                       <  3:194154/148‑1 (MQ=255)
                                                aCGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTg                                                                                                    <  4:134150/149‑1 (MQ=255)
                                                   gctgctGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGt                                                                                                  <  1:304903/149‑1 (MQ=255)
                                                   gctgctGCCCTGGTTGTTGCTGGCATGGATGTTGAGCTGCACACCCACACCGCTGGAAGTCACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGt                                                                                                  <  7:281409/149‑1 (MQ=255)
                                                     tgctgcCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTg                                                                                                >  4:19637/1‑149 (MQ=255)
                                                       ctgcCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTgg                                                                                               <  8:213868/148‑1 (MQ=255)
                                                        tgcCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCt                                                                                             >  3:66851/1‑149 (MQ=255)
                                                            cTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTgg                                                                                         <  1:35355/149‑1 (MQ=255)
                                                                   ttgCTGGCATTGATGTTGAGCTTCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGGCGGCGGCCTGTACCTGGGCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGt                                                                                   <  3:254603/148‑1 (MQ=255)
                                                                     gCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTg                                                                                 <  2:204947/148‑1 (MQ=255)
                                                                      cTGGCATTGATGTTGATCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGat                                                                               <  5:5053/149‑1 (MQ=255)
                                                                      cTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTgg                                                                                               <  3:193911/133‑1 (MQ=255)
                                                                      cTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTgg                                                                                               >  4:193911/1‑133 (MQ=255)
                                                                         gCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGtt                                                                                  <  7:177994/143‑1 (MQ=255)
                                                                         gCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGtt                                                                                  >  8:177994/1‑143 (MQ=255)
                                                                         gCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGatat                                                                             <  1:129721/148‑1 (MQ=255)
                                                                          cATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATc                                                                            >  4:97054/1‑148 (MQ=255)
                                                                                tGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttg                                                                     <  2:347538/149‑1 (MQ=255)
                                                                                 gTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgt                                                                    >  8:91503/1‑149 (MQ=255)
                                                                                  ttGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgtt                                                                   <  3:19637/149‑1 (MQ=255)
                                                                                       cTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAg                                                              >  4:306436/1‑149 (MQ=255)
                                                                                          cacacccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTGCGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAgcgg                                                           >  5:221932/1‑149 (MQ=255)
                                                                                           acacccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAgcggc                                                          >  2:30913/1‑149 (MQ=255)
                                                                                           acacccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTAGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgtt                                                                   <  1:147312/140‑1 (MQ=255)
                                                                                           acacccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTAGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgtt                                                                   >  2:147312/1‑140 (MQ=255)
                                                                                           acacccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATATTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAgcggc                                                          <  4:321918/149‑1 (MQ=255)
                                                                                              cccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCTGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCgtt                                                       <  8:154547/149‑1 (MQ=255)
                                                                                              cccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCgtt                                                       <  3:306436/149‑1 (MQ=255)
                                                                                               ccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCgttg                                                      >  1:359744/1‑149 (MQ=255)
                                                                                                       cTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGAGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCAcc                                              >  8:54044/1‑149 (MQ=255)
                                                                                                           aaGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCt                                          <  2:359744/149‑1 (MQ=255)
                                                                                                               cacacTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCAc                                      >  4:278774/1‑149 (MQ=255)
                                                                                                                acacTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACg                                     >  8:158908/1‑149 (MQ=255)
                                                                                                                  acTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGa                                   >  6:366062/1‑149 (MQ=255)
                                                                                                                      aGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGaccac                               >  4:333925/1‑149 (MQ=255)
                                                                                                                          cgcCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCTTTGTTGTCACCGGCTGCACGGACCACCTgg                           <  1:30913/149‑1 (MQ=255)
                                                                                                                          cgcCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTgg                           <  4:149056/149‑1 (MQ=255)
                                                                                                                                gcgcCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCACTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTa                       >  3:24424/1‑148 (MQ=255)
                                                                                                                                  gcCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTacac                    >  5:6685/1‑149 (MQ=255)
                                                                                                                                     aTTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCtt                 <  5:366062/149‑1 (MQ=255)
                                                                                                                                          gcaaggccAGGTACTGGCCGCCGGCCAGACCCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtg            >  3:327261/7‑149 (MQ=39)
                                                                                                                                          ccccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACTTACCACCTGGGTTACACCTTCggt             <  7:54044/148‑1 (MQ=255)
                                                                                                                                           cccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtg            <  7:91503/148‑1 (MQ=255)
                                                                                                                                           cccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGATGTTGTAAGCGGCGTTGTTGTCAc                                               <  2:227137/113‑1 (MQ=255)
                                                                                                                                           cccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGATGTTGTAAGCGGCGTTGTTGTCAc                                               >  1:227137/1‑113 (MQ=255)
                                                                                                                                           cccGACCAGGGTCTGGCCGGCCGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGGTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCAGGGTTACACCTTCggtgg           >  4:114714/1‑149 (MQ=255)
                                                                                                                                            ccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTacac                    >  3:98559/1‑139 (MQ=255)
                                                                                                                                            ccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTacac                    <  4:98559/139‑1 (MQ=255)
                                                                                                                                             aggccAGGTACTGGCCGCCGGCCAGACCCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggtg         >  1:253976/4‑149 (MQ=255)
                                                                                                                                               aCCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCtgtggtggt       <  4:24424/149‑1 (MQ=255)
                                                                                                                                                 cAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggtggt       <  3:251563/147‑1 (MQ=255)
                                                                                                                                                 cAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggtggt       >  4:251563/1‑147 (MQ=255)
                                                                                                                                                  aGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAg    <  4:66851/149‑1 (MQ=255)
                                                                                                                                                    ggTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAGGc  <  6:6685/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCAAGCAGCGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAGGC  >  NC_002947/3257324‑3257618

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