Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,562,245 C→T 6.5% G96D (GGT→GAT)  PP_3145 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,562,2450CT 6.5% 178.2 / 2.4 61G96D (GGT→GAT) PP_3145hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base T (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

TTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGAT  >  NC_002947/3562098‑3562390
                                                                                                                                                   |                                                                                                                                                  
ttGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCAcc                                                                                                                                                   >  3:330052/1‑149 (MQ=255)
            gTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCtt                                                                                                                                       <  3:35863/149‑1 (MQ=255)
             tGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTg                                                                                                                                      >  8:215877/1‑149 (MQ=255)
               aTTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTGGGCGGCGGCTTCGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGc                                                                                                                                    <  6:228506/149‑1 (MQ=255)
                       cAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGc                                                                                                                             <  7:215877/148‑1 (MQ=255)
                            tCTGCTGCGAGATACCCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCgg                                                                                                                       <  1:334528/149‑1 (MQ=255)
                            tCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCgg                                                                                                                       <  3:47563/149‑1 (MQ=255)
                                ctgcGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCt                                                                                                                     <  7:67956/147‑1 (MQ=255)
                                  gcGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAg                                                                                                                 >  1:369826/1‑149 (MQ=255)
                                  gcGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAg                                                                                                                 >  6:301888/1‑149 (MQ=255)
                                    gagaTAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTc                                                                                                               <  3:172086/149‑1 (MQ=255)
                                        tAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCttt                                                                                                            <  7:147352/148‑1 (MQ=255)
                                        tAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGCCTCTTGGCTTCTTGCAGCCGGCTGTAGTCttt                                                                                                            >  5:126184/1‑148 (MQ=255)
                                         aaCCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGcttctt                                                                                                                               >  4:240756/1‑128 (MQ=255)
                                         aaCCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGcttctt                                                                                                                               <  3:240756/128‑1 (MQ=255)
                                         aaCCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCttt                                                                                                            >  1:67/1‑147 (MQ=255)
                                          aCCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCATCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAg                                                                                                         >  2:349595/1‑149 (MQ=255)
                                           ccGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAgg                                                                                                        <  5:192741/149‑1 (MQ=255)
                                           ccGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAg                                                                                                         >  6:222820/1‑148 (MQ=255)
                                             gCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGaa                                                                                                      <  6:126184/149‑1 (MQ=255)
                                               tGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGc                                                                                                    >  5:39534/1‑149 (MQ=255)
                                                 ttGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCg                                                                                                  <  5:222820/149‑1 (MQ=255)
                                                    gCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGa                                                                                                <  4:330052/148‑1 (MQ=255)
                                                        cGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGctt                                                                                                                                  <  1:224062/110‑1 (MQ=255)
                                                        cGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGctt                                                                                                                                  >  2:224062/1‑110 (MQ=255)
                                                        cGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAAt                                                                                           <  6:39534/149‑1 (MQ=255)
                                                                ggTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCtt                                                                                                                                       <  8:224605/97‑1 (MQ=255)
                                                                ggTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCtt                                                                                                                                       >  7:224605/1‑97 (MQ=255)
                                                                     aaCTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCATCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTggg                                                                                    <  4:12956/143‑1 (MQ=255)
                                                                     aaCTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCATCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTggg                                                                                    >  3:12956/1‑143 (MQ=255)
                                                                               gctCGACGATCACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCCTCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTg                                                                    <  6:343660/149‑1 (MQ=255)
                                                                               gctCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTg                                                                    <  4:7096/149‑1 (MQ=255)
                                                                               gctCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTg                                                                    <  7:183030/149‑1 (MQ=255)
                                                                                 tCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCTCGCGGTTGCTg                                                                    >  1:338314/1‑147 (MQ=255)
                                                                                 tCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCTCGCGGTTGCTg                                                                    >  1:333868/1‑147 (MQ=255)
                                                                                      gaACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑gg                                                            >  3:306932/1‑149 (MQ=255)
                                                                                            ctagCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGACACACGACGGGCTCTTCGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCTCTTGCTGCAGAA‑GGCGACCa                                                      >  4:25289/4‑149 (MQ=255)
                                                                                            cTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCTGGTCACCCTACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCa                                                      <  1:305787/149‑1 (MQ=255)
                                                                                            cTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCa                                                      <  4:306932/149‑1 (MQ=255)
                                                                                            cTTGCTGAACTGGTTGATATTGACCTTGGGGTCGATCAAGCCCATCACCGGGTCACACGACGGGCTCATGGCTTCATGCAGCCCGCAGTAGTCTTTGAGGAAGCCGGAGTAATACTTCCGCTCCACGCGGTTCCTGGAGCA‑GGCGGCCa                                                      >  7:287725/1‑149 (MQ=255)
                                                                                                 tGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGc                                                  >  6:184531/1‑148 (MQ=255)
                                                                                                 tGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGc                                                  >  8:256543/1‑148 (MQ=255)
                                                                                                        gTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCATTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAgcg                                          >  2:293649/1‑149 (MQ=255)
                                                                                                          tGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑ggcg                                                          <  5:325707/131‑1 (MQ=255)
                                                                                                          tGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑ggcg                                                          >  6:325707/1‑131 (MQ=255)
                                                                                                           gATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCGACGGCAGGCAGCgcgg                                       >  4:198117/1‑149 (MQ=255)
                                                                                                            atatTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCgcggc                                      <  7:256543/149‑1 (MQ=255)
                                                                                                               ttGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGca                                   >  5:349235/1‑149 (MQ=255)
                                                                                                               ttGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGca                                   >  2:192742/1‑149 (MQ=255)
                                                                                                                     ttGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAg                             >  5:102534/1‑149 (MQ=255)
                                                                                                                          gTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTaa                        >  8:207139/1‑149 (MQ=255)
                                                                                                                           tCGATCCAGCGCATCACCGGGTCATCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAAt                       <  1:349595/149‑1 (MQ=255)
                                                                                                                             gATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAAGAGGGTAAtca                     >  4:68159/1‑147 (MQ=255)
                                                                                                                             gATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGc                     >  8:282724/1‑149 (MQ=255)
                                                                                                                               tCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGa                   <  6:102534/149‑1 (MQ=255)
                                                                                                                                 cAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAgg                 <  2:369826/149‑1 (MQ=255)
                                                                                                                                  aGCGCATCACCGGGTCACCCGACGGGCTCTTGGCGTCTTGCAGCCGGCTGTAGTCTTTGAGGCAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGAA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAgg                 >  6:81516/1‑148 (MQ=255)
                                                                                                                                       aTCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGa           >  4:176018/1‑149 (MQ=255)
                                                                                                                                        tCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGAt          >  3:205961/1‑149 (MQ=255)
                                                                                                                                              ggTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCatg    <  7:282724/149‑1 (MQ=255)
                                                                                                                                               gTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCATTGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCatg    <  1:293649/149‑1 (MQ=255)
                                                                                                                                                tCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCatgat  <  6:349235/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGAAGTAGCGGGTGATTTCCTGCAGCGTCTGCTGCGAGATAACCGCTGTTGGCTGCGGCTTGGGGTAGAACTGGCTTGGCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCACGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGAT  >  NC_002947/3562098‑3562390

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