Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,350,084 C→T 5.6% A853V (GCC→GTC) ‡ PP_4704 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,350,0840CT 5.6% 325.7 / 5.1 107A853V (GCC→GTC) ‡PP_4704hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (47/54);  new base T (3/3);  total (50/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01

Only 100 of 109 total aligned reads displayed.

CTCAGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAGCAG  >  NC_002947/5349936‑5350232
                                                                                                                                                    |                                                                                                                                                     
cTCAGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGc                                                                                                                                                       <  5:291786/149‑1 (MQ=255)
   aGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCAc                                                                                                                                                    >  7:177357/1‑149 (MQ=255)
     cggcgAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCt                                                                                                                                                  <  6:218021/149‑1 (MQ=255)
     ccgcgAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGCAACTGCGCCAGGAAATCGAAGGCACCAGAGCCCAGGCCAGCGACGGCCCGATCATAGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCt                                                                                                                                                  >  8:190030/3‑149 (MQ=255)
      ggcgAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTg                                                                                                                                                 <  2:347494/149‑1 (MQ=255)
       gcgAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTgg                                                                                                                                                >  4:419402/1‑149 (MQ=255)
         gAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTa                                                                                                                                              <  7:24391/149‑1 (MQ=255)
         gAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTa                                                                                                                                              >  4:20374/1‑149 (MQ=255)
         gAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTa                                                                                                                                              <  7:233885/149‑1 (MQ=255)
            gCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACaa                                                                                                                                           >  3:43143/1‑149 (MQ=255)
              ccAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGa                                                                                                                                         <  6:278994/149‑1 (MQ=255)
              ccAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGa                                                                                                                                         <  8:211858/149‑1 (MQ=255)
                 gtcgtcCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACa                                                                                                                                            <  2:79675/143‑1 (MQ=255)
                 gtcgtcCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACa                                                                                                                                            >  1:79675/1‑143 (MQ=255)
                    gtcCTTGGACAAGCTGACATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGTGTACCAGGCCACCTGGTACAAGATCCTg                                                                                                                                    <  1:246745/148‑1 (MQ=255)
                     tcCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCa                                                                                                                                  <  6:78021/149‑1 (MQ=255)
                          ggACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCg                                                                                                                               <  3:152433/147‑1 (MQ=255)
                            aCAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACgg                                                                                                                           >  2:81486/1‑149 (MQ=255)
                             cAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGt                                                                                                                                               >  3:16127/1‑128 (MQ=255)
                             cAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGt                                                                                                                                               <  4:16127/128‑1 (MQ=255)
                             cAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACggg                                                                                                                          <  5:404293/149‑1 (MQ=255)
                              aaGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACgggg                                                                                                                         >  3:38140/1‑149 (MQ=255)
                              aaGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACgggg                                                                                                                         >  5:114897/1‑149 (MQ=255)
                              aaGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACgggg                                                                                                                         >  7:378601/1‑149 (MQ=255)
                               aGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGt                                                                                                                        <  6:114897/149‑1 (MQ=255)
                                  tGGCATCCGCGCTGAAAGAGCGCACGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGACCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGa                                                                                                                      >  8:11107/1‑148 (MQ=255)
                                     cATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGTGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAg                                                                                                                  <  1:333082/149‑1 (MQ=255)
                                      aTCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAgg                                                                                                                 >  4:193342/1‑149 (MQ=255)
                                      aTCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAgg                                                                                                                 >  7:71057/1‑149 (MQ=255)
                                       tCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGTGACGGCCCGCTGATTGCCCAACTGCGCCTTGAGCGTGAGTACCAGGTGGCCTGGTActc                                                                                                                                           >  1:112099/1‑120 (MQ=255)
                                       tCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGTGACGGCCCGCTGATTGCCCAACTGCGCCTTGAGCGTGAGTACCAGGTGGCCTGGTActc                                                                                                                                           <  2:112099/122‑3 (MQ=255)
                                       tCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGTGACGGCCCGCTGATTGCCCAACTGCGCCTTGAGCGTGAGTACCAGGTGGCCTGGTActc                                                                                                                                           >  1:106306/1‑120 (MQ=255)
                                       tCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGTGACGGCCCGCTGATTGCCCAACTGCGCCTTGAGCGTGAGTACCAGGTGGCCTGGTActc                                                                                                                                           <  2:106306/122‑3 (MQ=255)
                                       tCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGt                                                                                                                >  7:166019/1‑149 (MQ=255)
                                          gcgcTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGcc                                                                                                             >  5:43919/1‑149 (MQ=255)
                                           cgcTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCg                                                                                                            <  3:358094/149‑1 (MQ=255)
                                           cgcTGAAAGAGCGCCCGGAACGGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCg                                                                                                            <  4:297925/149‑1 (MQ=255)
                                            gcTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCgg                                                                                                           <  3:235024/149‑1 (MQ=255)
                                             cTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGc                                                                                                          <  4:43143/149‑1 (MQ=255)
                                             cTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGc                                                                                                          <  4:69260/149‑1 (MQ=255)
                                             cTGAAAGAGCGCACGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGc                                                                                                          <  5:303595/149‑1 (MQ=255)
                                              tGAAAGAGCGCCCGGAACTGCGCCTCGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGc                                                                                                          >  4:202116/1‑148 (MQ=255)
                                               gAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGcc                                                                                                         >  8:179092/1‑148 (MQ=255)
                                                aaaGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCaa                                                                                                       >  2:320463/1‑149 (MQ=255)
                                                         ccGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGAt                                                                                              <  6:43919/149‑1 (MQ=255)
                                                           ggAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGc                                                                                            <  2:283787/149‑1 (MQ=255)
                                                              aCTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTTAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtgg                                                                                         <  5:42061/149‑1 (MQ=255)
                                                                tGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtggtg                                                                                       >  3:346820/1‑149 (MQ=255)
                                                                tGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGCCCACCTGGTACAAGATCCTGCAGAGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtggtg                                                                                       >  5:227157/1‑149 (MQ=255)
                                                                 gcgcCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGAGAAGGTGCCGGCCAATGCCTCGATGCtggtgg                                                                                      >  8:308485/1‑149 (MQ=255)
                                                                  cgcgTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCTGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggt                                                                                    <  7:190030/149‑1 (MQ=255)
                                                                  cgcCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtggtg                                                                                       <  7:281111/147‑1 (MQ=255)
                                                                    ccTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGTCCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtggtggtgg                                                                                  <  1:391547/149‑1 (MQ=255)
                                                                       ggAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGgac                                                                                >  3:86516/1‑148 (MQ=255)
                                                                        gAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGgac                                                                                <  7:179092/147‑1 (MQ=255)
                                                                         aaaTCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGgacga                                                                              >  2:271863/1‑148 (MQ=255)
                                                                          aaTCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACa                                                                            <  1:406180/149‑1 (MQ=255)
                                                                            tCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGt                                                                                                                        >  7:71374/1‑104 (MQ=255)
                                                                            tCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGt                                                                                                                        <  8:71374/104‑1 (MQ=255)
                                                                            tCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGc                                                                          <  8:121399/149‑1 (MQ=255)
                                                                               aaGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGa                                                                        >  6:61153/1‑148 (MQ=255)
                                                                                 ggCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACaa                                                                     <  3:20374/149‑1 (MQ=255)
                                                                                  gCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGCTCCTGCAGCGACGGGGTGACAAGGTGCCGCCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAg                                                                    >  2:101066/1‑149 (MQ=255)
                                                                                     ccAGCGCCCAGTCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGccc                                                                 <  4:346820/149‑1 (MQ=255)
                                                                                        gcgcCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCcgcg                                                              <  3:419402/149‑1 (MQ=255)
                                                                                        gcgcCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCcgc                                                               >  3:387822/1‑148 (MQ=255)
                                                                                         cgcCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGa                                                             <  3:193342/149‑1 (MQ=255)
                                                                                           cccAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATg                                                           <  4:348996/149‑1 (MQ=255)
                                                                                           cccAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATg                                                           <  1:271863/149‑1 (MQ=255)
                                                                                            ccAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGc                                                          >  1:383477/1‑149 (MQ=255)
                                                                                              aGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTc                                                        <  8:378601/149‑1 (MQ=255)
                                                                                                 ccAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGaa                                                     <  6:130799/149‑1 (MQ=255)
                                                                                                     cGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCa                                                 >  2:419936/1‑149 (MQ=255)
                                                                                                       aCGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCAtt                                               <  3:259271/149‑1 (MQ=255)
                                                                                                         ggCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTa                                             >  5:414743/1‑149 (MQ=255)
                                                                                                          gCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGTCCGCGATGCTCGAAGGCATTTAc                                            <  7:308485/149‑1 (MQ=255)
                                                                                                           cccGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTCCCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGCTGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCGCGAAGGCATTTAcc                                           <  2:327499/149‑1 (MQ=255)
                                                                                                                  aTTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCg                                    >  4:53005/1‑149 (MQ=255)
                                                                                                                   ttGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGCTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTAccc                                     <  1:320463/148‑1 (MQ=255)
                                                                                                                      cccAACAGCGCCTGGAGCGTGAGTACCACGCCACCTGGTACAAGATCCTGCAGCGACTGTGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTg                                <  8:366101/149‑1 (MQ=255)
                                                                                                                       ccAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTg                                >  2:252268/1‑148 (MQ=255)
                                                                                                                         aaCAGCGCCTTGAGCGTGAGTACCTGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCTATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGaaa                             <  3:53005/149‑1 (MQ=255)
                                                                                                                            aGCGCCTTGAGCGTGAGTACCAGGCCCCCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCtggtg                                                                                       >  8:173418/1‑89 (MQ=255)
                                                                                                                              cgcCTTGAGCGGGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAAcagca                        <  4:38019/149‑1 (MQ=255)
                                                                                                                                 cTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAgcc                     <  3:202116/149‑1 (MQ=255)
                                                                                                                                    gAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGGAACAGCAgccgc                   <  4:38140/148‑1 (MQ=255)
                                                                                                                                     aGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCt                 <  1:394727/149‑1 (MQ=255)
                                                                                                                                         tGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCAACAAGCCAGCGATGCTCGAAGGAATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGcc              >  1:140628/1‑148 (MQ=255)
                                                                                                                                          gAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCg             >  5:374987/1‑148 (MQ=255)
                                                                                                                                            gTACCAGGTGGCCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAgccgc                   >  6:303548/1‑140 (MQ=255)
                                                                                                                                            gTACCAGGTGGCCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAgccgc                   <  5:303548/140‑1 (MQ=255)
                                                                                                                                            gTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAAt          <  6:414743/149‑1 (MQ=255)
                                                                                                                                             tACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGATTg         >  4:388291/1‑149 (MQ=255)
                                                                                                                                             tACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATg         >  7:102776/1‑149 (MQ=255)
                                                                                                                                               ccAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAg                                                    >  6:187776/1‑104 (MQ=255)
                                                                                                                                               ccAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAg                                                    <  5:187776/104‑1 (MQ=255)
                                                                                                                                               ccAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGGCGCCTGCCGAATggg       <  2:383477/149‑1 (MQ=255)
                                                                                                                                               ccAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATggg       <  1:419936/149‑1 (MQ=255)
                                                                                                                                                  ggCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAg     >  5:143677/1‑148 (MQ=255)
                                                                                                                                                    ccACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCAGTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAAAGTGAGCAg  >  8:69838/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTCAGCGGCGAGGCCCAGTCGTCCTTGGACAAGCTGGCATCCGCGCTGAAAGAGCGCCCGGAACTGCGCCTGGAAATCGAAGGCACCAGCGCCCAGGCCAGCGACGGCCCGCTGATTGCCCAACAGCGCCTTGAGCGTGAGTACCAGGCCACCTGGTACAAGATCCTGCAGCGACGGGGTGACAAGGTGCCGGCCAATGCCTCGATGCTGGTGGTGGACGACAGCGACAAGCCCGCGATGCTCGAAGGCATTTACCGTACCCGCCTGAAACAGCAGCCGCCTGCCGAATGGGAGCAG  >  NC_002947/5349936‑5350232

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: