Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 396,856 C→G 6.8% L30L (CTG→CTC) ‡ PP_0329 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947396,8560CG 6.8% 284.1 / 8.6 88L30L (CTG→CTC) ‡PP_0329hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/42);  new base G (4/2);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

TATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCAC  >  NC_002947/396710‑396997
                                                                                                                                                   |                                                                                                                                              
tATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAg                                                                                                                                              <  7:173519/149‑1 (MQ=255)
  tCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAgg                                                                                                                                             <  5:367319/148‑1 (MQ=255)
    cAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGAtt                                                                                                                                          <  1:229788/149‑1 (MQ=255)
     aaaGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGAtt                                                                                                                                          >  5:425594/1‑148 (MQ=255)
     aaaGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGAtt                                                                                                                                          >  3:111535/1‑148 (MQ=255)
     aaaGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGAtt                                                                                                                                          >  8:392804/1‑148 (MQ=255)
          gggTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGc                                                                                                                                   <  7:391373/149‑1 (MQ=255)
            gtgtAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCAATCCTGCGGATTGGAACGCtc                                                                                                                                 <  1:176971/149‑2 (MQ=255)
             tgtAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTc                                                                                                                                <  6:269855/149‑1 (MQ=255)
                       aTTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGttt                                                                                                                       <  3:115837/148‑1 (MQ=255)
                         ttttttGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTcc                                                                                                                     >  5:105844/1‑148 (MQ=255)
                           ttttGCTTTCCGTCAGGCGTTCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGTAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCaa                                                                                                                   <  6:187997/148‑1 (MQ=255)
                            tttGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAAt                                                                                                                  <  3:30405/148‑1 (MQ=255)
                              tGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGa                                                                                                                >  5:251835/1‑148 (MQ=255)
                              tGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGa                                                                                                                <  6:425594/148‑1 (MQ=255)
                               gCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGa                                                                                                                >  3:400383/1‑147 (MQ=255)
                               gCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCAAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGa                                                                                                                >  2:237210/1‑147 (MQ=255)
                                cTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGTCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGa                                                                                                                >  1:119751/1‑147 (MQ=255)
                                cTTTCCGTCAGGCGTGCTAACGACAAGCACTGAGTCCATCTCATCTGTTGCTAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCTTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACa                                                                                                              <  2:352560/148‑1 (MQ=255)
                                 tttCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGt                                                                                                            >  2:438044/1‑149 (MQ=255)
                                  ttCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTg                                                                                                           >  5:325310/1‑149 (MQ=255)
                                   tCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGACCTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTg                                                                                                           <  1:30587/148‑1 (MQ=255)
                                      gTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                       <  5:157267/149‑1 (MQ=255)
                                           gCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTg                                                                                                  >  3:339231/1‑149 (MQ=255)
                                            cGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGc                                                                                                 <  7:392804/149‑1 (MQ=255)
                                              tGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCa                                                                                                                                               <  4:236207/102‑1 (MQ=255)
                                              tGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCa                                                                                                                                               >  3:236207/1‑102 (MQ=255)
                                              tGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctg                                                                                               <  2:119751/149‑1 (MQ=255)
                                                 tCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTTGTCAGTTTGctggtc                                                                                            >  4:262728/1‑149 (MQ=255)
                                                 tCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtc                                                                                            >  5:64892/1‑149 (MQ=255)
                                                 tCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtc                                                                                            >  2:379175/1‑149 (MQ=255)
                                                   aCGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctg                                                                                          <  6:251835/149‑1 (MQ=255)
                                                      aCAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctgg                                                                                         >  5:387002/1‑147 (MQ=255)
                                                            aCTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAata                                                                                 >  1:435956/1‑149 (MQ=255)
                                                              tGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATAtt                                                                               >  1:153784/1‑149 (MQ=255)
                                                              tGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTCGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctggt                                                                                        <  5:176618/140‑1 (MQ=255)
                                                              tGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTCGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctggt                                                                                        >  6:176618/1‑140 (MQ=255)
                                                                 gTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgg                                                                            <  6:64892/149‑1 (MQ=255)
                                                                    catctcatctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgggg                                                                          <  1:379175/148‑1 (MQ=255)
                                                                    catctcatctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgggg                                                                          <  4:400383/148‑1 (MQ=255)
                                                                          atctGTTGCGAAAAGCTGGTAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCa                                                                    <  2:322787/148‑1 (MQ=255)
                                                                          atctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACa                                                                                                              <  7:249401/106‑1 (MQ=255)
                                                                          atctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACa                                                                                                              >  8:249401/1‑106 (MQ=255)
                                                                          atctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAg                                                                   >  4:71318/1‑149 (MQ=255)
                                                                             tGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGAc                                                                >  8:144608/1‑149 (MQ=255)
                                                                              gTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACt                                                               >  1:407518/1‑149 (MQ=255)
                                                                               ttGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACt                                                               <  2:435956/148‑1 (MQ=255)
                                                                                  cGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGttt                                                           <  4:111535/149‑1 (MQ=255)
                                                                                  cGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGttt                                                           <  6:325310/149‑1 (MQ=255)
                                                                                    aaaaGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGttttt                                                         <  6:105844/149‑1 (MQ=255)
                                                                                     aaaGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTc                                                        <  4:169878/149‑1 (MQ=255)
                                                                                           tGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGa                                                   >  2:124149/1‑149 (MQ=255)
                                                                                            gggAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGAc                                                  >  6:223847/1‑149 (MQ=255)
                                                                                             ggAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGAc                                                  <  3:71318/148‑1 (MQ=255)
                                                                                              gAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAAACCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGt                                                >  1:102967/1‑149 (MQ=255)
                                                                                               aaCTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTg                                               >  1:424483/1‑149 (MQ=255)
                                                                                                  ttCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                            >  1:414529/1‑149 (MQ=255)
                                                                                                  ttCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                            >  2:69636/1‑149 (MQ=255)
                                                                                                    ccTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                            <  5:223847/147‑1 (MQ=255)
                                                                                                     cTGCTTCTTTGCCAGCGGACTAAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCg                                         >  4:3466/1‑149 (MQ=255)
                                                                                                       gCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGtt                                       >  8:70548/1‑149 (MQ=255)
                                                                                                             ttGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                       >  1:268663/1‑79 (MQ=255)
                                                                                                             ttGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                       <  2:268663/79‑1 (MQ=255)
                                                                                                              tGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCt                                   <  4:339231/146‑1 (MQ=255)
                                                                                                                 ctgCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGtt                             >  6:430617/3‑149 (MQ=255)
                                                                                                                  aGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTg                            <  6:387002/149‑1 (MQ=255)
                                                                                                                   gCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTgg                           <  1:438044/149‑1 (MQ=255)
                                                                                                                       aCTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                            <  6:195568/128‑1 (MQ=255)
                                                                                                                       aCTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                            >  5:195568/1‑128 (MQ=255)
                                                                                                                         tCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCgg                     <  4:423302/149‑1 (MQ=255)
                                                                                                                           aGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCgg                     <  3:262728/147‑1 (MQ=255)
                                                                                                                            gCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaaa                  >  4:10894/1‑149 (MQ=255)
                                                                                                                             cAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGt                                     <  6:368227/129‑1 (MQ=255)
                                                                                                                             cAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGt                                     >  5:368227/1‑129 (MQ=255)
                                                                                                                             cAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAAt                 <  2:407518/149‑1 (MQ=255)
                                                                                                                               agagGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAg               <  2:424483/149‑1 (MQ=255)
                                                                                                                                gagGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaaa                  >  6:387164/1‑145 (MQ=255)
                                                                                                                                  gggAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGtt                                       <  7:321488/122‑1 (MQ=255)
                                                                                                                                  gggAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGtt                                       >  8:321488/1‑122 (MQ=255)
                                                                                                                                  gggAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATg            <  2:342285/149‑1 (MQ=255)
                                                                                                                                   gaaGCGTTCAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATgg           >  7:64209/3‑149 (MQ=255)
                                                                                                                                   gaaGCGTTCAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATgg           >  3:324127/3‑149 (MQ=255)
                                                                                                                                    aaGCGTTCAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaa                   >  4:86229/2‑140 (MQ=255)
                                                                                                                                    aaGCGTTCAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaa                   <  3:86229/139‑1 (MQ=255)
                                                                                                                                    aaGCGTTCAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCGCCCGGTGTCTGCAGTTGGCCCCGGAAATAGATGGc          >  8:113789/2‑149 (MQ=255)
                                                                                                                                     aGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGcc         >  2:26220/1‑149 (MQ=255)
                                                                                                                                      gCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCt        <  3:10894/149‑1 (MQ=255)
                                                                                                                                             aaTCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCAc  >  2:76940/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATTGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCAC  >  NC_002947/396710‑396997

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: