Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,063,627 C→A 7.2% D238Y (GAC→TAC)  PP_3581 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,063,6270CA 7.2% 221.8 / 5.9 69D238Y (GAC→TAC) PP_3581membrane protein
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base A (3/2);  total (36/33)
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.53e-01

TGGCGATGGCCTGTTCGATGCGATCGGCCAGCAGGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTGCCATGCCGGTGTAA  >  NC_002947/4063480‑4063772
                                                                                                                                                   |                                                                                                                                                 
tGGCGATGGCCTGTTCGATGCGATCGGCCAGCAGGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTc                                                                                                                                                   <  8:2452/148‑1 (MQ=255)
       ggCCTGTTCGATGCGATCGGCCAGCAGGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCAcc                                                                                                                                            >  5:35366/1‑148 (MQ=255)
             ttCGATGCGATCGGCCAGCAGGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGt                                                                                                                                       <  2:280908/147‑1 (MQ=255)
                gATGCGATCGGCCAGCAGGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGc                                                                                                                                   <  6:35366/148‑1 (MQ=255)
                 aTGCGATCGGCCAGCAGGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCa                                                                                                                                  >  7:270088/1‑148 (MQ=255)
                     gATCGGCCAGCAGGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAg                                                                                                                                    >  6:90142/1‑142 (MQ=255)
                     gATCGGCCAGCAGGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAg                                                                                                                                    <  5:90142/142‑1 (MQ=255)
                      aTCGGCCAGCAGGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGcc                                                                                                                             >  1:233579/1‑148 (MQ=255)
                       tCGGCCAGCAGGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCAc                                                                                                                           >  8:205746/1‑149 (MQ=255)
                          gCCAGCAGGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGc                                                                                                                        >  7:158186/1‑149 (MQ=255)
                            cagcagGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCa                                                                                                                                  >  4:45166/1‑137 (MQ=255)
                             agcagGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAg                                                                                                                     <  5:102604/149‑1 (MQ=255)
                                   gTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGc                                                                                                               >  3:192866/1‑149 (MQ=255)
                                       ggTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGc                                                                                                                        <  5:77111/136‑1 (MQ=255)
                                       ggTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGc                                                                                                                        >  6:77111/1‑136 (MQ=255)
                                        gTTGGGCAGATCGGTGAGTGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCa                                                                                                             <  4:50904/146‑1 (MQ=255)
                                        gTTGGGCAGATCGGTGAGTGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCa                                                                                                             >  3:50904/1‑146 (MQ=255)
                                           gggCAGATCGGTGAGGGTGTCGTGCAAGGGCAGTTGGGTCAGTTCCTGGTTGGACAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCgtg                                                                                                        <  1:288197/148‑1 (MQ=255)
                                           gggCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCgtg                                                                                                        <  2:233579/148‑1 (MQ=255)
                                            ggCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCgtggt                                                                                                      <  6:211923/149‑1 (MQ=255)
                                            ggCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCgtgg                                                                                                       <  3:283083/148‑1 (MQ=255)
                                             gCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCgtggt                                                                                                      >  3:304719/1‑148 (MQ=255)
                                             gCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCgtggt                                                                                                      >  7:188020/1‑148 (MQ=255)
                                              cAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGa                                                                                                    >  6:99914/1‑149 (MQ=255)
                                                 aTCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCa                                                                                                 <  5:99914/149‑1 (MQ=255)
                                                    ggTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCa                                                                                              >  7:158187/1‑149 (MQ=255)
                                                             gtCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGc                                                                                                                        <  6:319613/114‑1 (MQ=255)
                                                             gtCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGc                                                                                                                        >  5:319613/1‑114 (MQ=255)
                                                                gTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAAc                                                                                    >  6:39954/1‑147 (MQ=255)
                                                                 tGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCa                                                                                                             <  5:146253/121‑1 (MQ=255)
                                                                 tGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCa                                                                                                             >  6:146253/1‑121 (MQ=255)
                                                                  gCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGTACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAgg                                                                                >  1:200664/1‑149 (MQ=255)
                                                                   cAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAgg                                                                                <  1:115356/148‑1 (MQ=255)
                                                                   cAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAgg                                                                                <  4:304719/148‑1 (MQ=255)
                                                                    aaGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCt                                                                              <  5:6533/149‑1 (MQ=255)
                                                                      ggCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGt                                                                            >  7:255872/1‑149 (MQ=255)
                                                                         cAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTccc                                                                         >  6:87189/1‑149 (MQ=255)
                                                                         cAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTccc                                                                         <  6:114183/149‑1 (MQ=255)
                                                                           gTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTccccc                                                                       >  5:149145/1‑149 (MQ=255)
                                                                               ggTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGc                                                                    >  6:273626/1‑148 (MQ=255)
                                                                                        cTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACccgccg                                                           <  6:149145/148‑1 (MQ=255)
                                                                                           gTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGTACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCa                                                       <  2:200664/149‑1 (MQ=255)
                                                                                           gTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCgccgcc                                                         >  7:6003/1‑147 (MQ=255)
                                                                                              ggCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAgcg                                                    >  6:107684/1‑149 (MQ=255)
                                                                                               gCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAgcgc                                                   >  5:199321/1‑149 (MQ=255)
                                                                                                      tCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCa                                             >  6:232704/1‑148 (MQ=255)
                                                                                                             ccGGGCCAACTCGGCTGTACGAGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGc                                     <  4:174651/148‑1 (MQ=255)
                                                                                                                 gCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTc                                 >  3:105574/1‑149 (MQ=255)
                                                                                                                 gCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTc                                 >  3:110267/1‑149 (MQ=255)
                                                                                                                 gCCAACTCGGCTGTACGTGCCTCGAGACGCGCCTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCTCCCAGCGCGCCGCAAAAGCTGCCCTc                                 <  1:10178/149‑1 (MQ=255)
                                                                                                                   cAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAgg                                                                                >  4:109017/1‑100 (MQ=255)
                                                                                                                   cAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAgg                                                                                <  3:109017/100‑1 (MQ=255)
                                                                                                                   cAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCgg                               <  4:6125/149‑1 (MQ=255)
                                                                                                                   cAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCgg                               <  7:205746/149‑1 (MQ=255)
                                                                                                                   cAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCgg                               <  8:158186/149‑1 (MQ=255)
                                                                                                                       ggtgCTGGACGCGCGTCTCGAGGCACGTACAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGCCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGa                            >  5:303092/4‑148 (MQ=255)
                                                                                                                         tgCTGGACGCGCGTCTCGAGGCACGTACAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAAt                         >  3:202336/2‑149 (MQ=255)
                                                                                                                         tgCTGGACGCGCGTCTCGAGGCACGTACAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAAt                         >  2:271095/2‑149 (MQ=255)
                                                                                                                              tACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTcgc                     <  5:107684/148‑1 (MQ=255)
                                                                                                                                  tGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTg                <  5:273626/149‑1 (MQ=255)
                                                                                                                                   gCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTGc               >  2:158556/1‑149 (MQ=255)
                                                                                                                                      tCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTGCCAt            >  3:227064/1‑149 (MQ=255)
                                                                                                                                      tCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTGCCAt            >  8:275165/1‑149 (MQ=255)
                                                                                                                                        gagaCGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTGCCATGc          <  5:39954/149‑1 (MQ=255)
                                                                                                                                          gaggcacgTACAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTGCCATGCCg        <  1:271095/143‑1 (MQ=255)
                                                                                                                                          gaCGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTGCCATGCCg        <  7:41466/149‑1 (MQ=255)
                                                                                                                                            cgcgcgTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTGCCATGCCGgt      <  5:87189/149‑1 (MQ=255)
                                                                                                                                               gcgTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTGCCATGcc         <  7:6073/143‑1 (MQ=255)
                                                                                                                                               gcgTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTGCCATGcc         >  8:6073/1‑143 (MQ=255)
                                                                                                                                                cgTACAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGCCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTGCCATGCCGGTGTaa  <  6:303092/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGGCGATGGCCTGTTCGATGCGATCGGCCAGCAGGGTACGGTTGGGCAGATCGGTGAGGGTGTCGTGCAAGGCCAGTTGGGTCAGTTCCTGGTTGGCCAAGGTCAGCGAACGGGCCAACTCGGCTGTACGTGCCTCGAGACGCGCGTCCAGCACCGAAGTCAGCAAAGCCACTGCCAGTACTGCCAGCGTGGTGATCAGCACTAGGTAAACCAGGCTGTCCCCCTGCAACCCGCCGCCCAGCGCGCCGCAAAAGCTGCCCTCGGGGAAATTCGCCGCTGCCATGCCGGTGTAA  >  NC_002947/4063480‑4063772

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