Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,633,152 C→A 7.5% N208K (AAC→AAA) ‡ PP_3201 → BNR domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,633,1520CA 7.5% 184.3 / 4.1 61N208K (AAC→AAA) ‡PP_3201BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (25/31);  new base A (3/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

TTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTCTGCTGCTCGGCGGCCTGCGC  >  NC_002947/3633005‑3633299
                                                                                                                                                   |                                                                                                                                                   
tttCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAAc                                                                                                                                                     >  6:24868/1‑148 (MQ=255)
    aCCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGgt                                                                                                                                                >  1:113320/1‑149 (MQ=255)
       gATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGGTTGCCTTGGCGCGCCAGGGCATACAGgtgca                                                                                                                                             <  4:104361/149‑3 (MQ=38)
        ctGCCCGCCACGGCCTGGTCGTCGGTGCCTACGGGGGGGCCATGCGCACCGACGACGGCGTGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTAc                                                                                                                                            <  6:77211/148‑1 (MQ=255)
            ccGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCg                                                                                                                                        >  3:297685/1‑149 (MQ=255)
            ccGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGGTTGCCTTGGCGCGCCAGGGCATACAGgtg                                                                                                                                               >  1:8503/1‑142 (MQ=39)
            ccGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGGTTGCCTTGGCGCGCCAGGGCATACAGgtg                                                                                                                                               <  2:8503/142‑1 (MQ=39)
             cGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGc                                                                                                                                       >  6:40846/1‑149 (MQ=255)
                 aCGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGc                                                                                                                                   <  5:24868/149‑1 (MQ=255)
                      cTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAg                                                                                                                                >  8:153774/1‑147 (MQ=255)
                      cTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAg                                                                                                                                >  3:353941/1‑147 (MQ=255)
                         gtcgtcGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAgg                                                                                                                            <  5:322281/148‑1 (MQ=255)
                          tcgtcgGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGc                                                                                                                          <  6:191098/149‑1 (MQ=255)
                          tcgtcgGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGc                                                                                                                          <  1:368240/149‑1 (MQ=255)
                            gtcgGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCct                                                                                                                        >  6:62993/1‑149 (MQ=255)
                                    cTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTATATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTtc                                                                                                                <  4:210464/149‑1 (MQ=255)
                                     tACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTtcg                                                                                                               >  5:184319/1‑149 (MQ=255)
                                            tGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGt                                                                                                        <  7:123129/149‑1 (MQ=255)
                                             ggCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGt                                                                                                        >  2:328932/1‑148 (MQ=255)
                                                cATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATg                                                                                                    <  5:40846/149‑1 (MQ=255)
                                                cATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATg                                                                                                    <  5:311934/149‑1 (MQ=255)
                                                 aTGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATgg                                                                                                   >  3:138630/1‑149 (MQ=255)
                                                    cgcACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCg                                                                                                 <  4:353941/148‑1 (MQ=255)
                                                     gcACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGt                                                                                               >  7:201289/1‑149 (MQ=255)
                                                                      cgcAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCt                                                                              <  2:95017/149‑1 (MQ=255)
                                                                           cTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACg                                                                         >  4:371487/1‑149 (MQ=255)
                                                                            tGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACgg                                                                        >  2:250264/1‑149 (MQ=255)
                                                                             ggCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACggg                                                                       >  2:169718/1‑149 (MQ=255)
                                                                              gCAGTCGGAGCTCGGCCGACGGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGc                                                                      <  5:265918/149‑1 (MQ=255)
                                                                              gCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACggg                                                                       <  8:201289/148‑1 (MQ=255)
                                                                                 gTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCt                                                                   <  4:297685/149‑1 (MQ=255)
                                                                                    ggAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTAc                                                                 <  5:212438/148‑1 (MQ=255)
                                                                                     gagaTCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGa                                                               <  6:184319/149‑1 (MQ=255)
                                                                                      agaTGGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCGTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGAt                                                              <  6:82689/149‑1 (MQ=255)
                                                                                      agaTCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGAt                                                              <  3:371487/149‑1 (MQ=255)
                                                                                      agaTCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGAt                                                              <  3:203129/149‑1 (MQ=255)
                                                                                        aTCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATgg                                                            >  7:270368/1‑149 (MQ=255)
                                                                                           ggCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCa                                                          <  3:22108/148‑1 (MQ=255)
                                                                                             ccGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGc                                                        >  4:94825/1‑148 (MQ=255)
                                                                                                       aaCGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTcgccgc                                             >  3:256030/1‑149 (MQ=255)
                                                                                                                 gACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGtt                                    <  3:94825/148‑1 (MQ=255)
                                                                                                                 gACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGtt                                    <  1:250264/148‑1 (MQ=255)
                                                                                                                  aCTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGc                                  >  6:57189/1‑149 (MQ=255)
                                                                                                                   cTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGcc                                 <  2:113320/149‑1 (MQ=255)
                                                                                                                    tGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGccc                                >  2:218380/1‑149 (MQ=255)
                                                                                                                        ccTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATg                            <  7:153774/149‑1 (MQ=255)
                                                                                                                         ctggttGCCTTGGCGCGCCAGGGCATACAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATgg                           >  1:229005/6‑149 (MQ=255)
                                                                                                                         ctggttGCCTTGGCGCGCCAGGGCATACAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCAGATgg                           >  6:211853/6‑149 (MQ=255)
                                                                                                                            tATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCg                        <  4:138630/149‑1 (MQ=255)
                                                                                                                              tGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTc                      <  1:218380/149‑1 (MQ=255)
                                                                                                                              tGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTc                      >  2:266212/1‑149 (MQ=255)
                                                                                                                                 ccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCAt                                                                                         >  2:54251/1‑79 (MQ=255)
                                                                                                                                 ccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCAt                                                                                         <  1:54251/79‑1 (MQ=255)
                                                                                                                                 ccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCACTGCACGGGCCCTACGATGGCAGCTACTTcgccg                                              <  7:226955/122‑1 (MQ=255)
                                                                                                                                 ccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCACTGCACGGGCCCTACGATGGCAGCTACTTcgccg                                              >  8:226955/1‑122 (MQ=255)
                                                                                                                                     gcgcgcCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTctgctg                 <  5:62993/147‑1 (MQ=255)
                                                                                                                                     gcgcgcCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTCTgctgct               >  7:72058/1‑149 (MQ=255)
                                                                                                                                         gcCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTCTGCTGCTcggc           >  1:83632/1‑149 (MQ=255)
                                                                                                                                              ggCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTCTGCTGCTCGGCGGCCt      <  8:270368/149‑1 (MQ=255)
                                                                                                                                               gCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTCTGCTGCTCGGCGGCCTg     >  5:226368/1‑149 (MQ=255)
                                                                                                                                               gCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTCTGCTGCTCGGCGGCCTg     >  2:209878/1‑149 (MQ=255)
                                                                                                                                                  aCCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTCTGCTGCTCGGCGGCCTgcgc  <  1:266212/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTCTGCTGCTCGGCGGCCTGCGC  >  NC_002947/3633005‑3633299

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