Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 823,548 A→C 5.3% H352Q (CAT→CAG)  PP_0708 ← betaine‑aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947823,5480AC 5.3% 232.1 / 2.6 74H352Q (CAT→CAG) PP_0708betaine‑aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base C (2/2);  total (41/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 = 9.98e-01

TTGGCTTCGGCAATGGCCTGTTCCTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCGCAGGTTTTCCAGTGCCGCC  >  NC_002947/823402‑823695
                                                                                                                                                   |                                                                                                                                                    
ttGGCTTCGGCAATGGCCTGTTCCTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ATATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATg                                                                                                                                                    >  7:273209/1‑149 (MQ=255)
     ttCGGCAATGGCCTGTTCCTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGgccgc                                                                                                                                               >  5:270149/1‑149 (MQ=255)
      tCGGCAATGGCCTGTTCCTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGgccgc                                                                                                                                               >  2:208245/1‑148 (MQ=255)
        ggCAATGGCCTGTTCCTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCgcca                                                                                                                                             <  7:286134/148‑1 (MQ=255)
         gCAATGGCCTGTTCCTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCgccagc                                                                                                                                           <  1:208245/149‑1 (MQ=255)
          cAATGGCCTGTTCCTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCgccagcc                                                                                                                                          <  2:337349/149‑1 (MQ=255)
                ccTGTTCCTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGCACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCAc                                                                                                                                    >  6:185978/1‑149 (MQ=255)
                 cTGTTCCTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACt                                                                                                                                   <  3:256798/149‑1 (MQ=255)
                  tGTTCCTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTg                                                                                                                                  >  2:333612/1‑149 (MQ=255)
                    ttCCTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGc                                                                                                                                 >  6:111424/1‑148 (MQ=255)
                       cTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCaaa                                                                                                                              <  5:244698/148‑1 (MQ=255)
                       cTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAAc                                                                                                                             <  6:190563/149‑1 (MQ=255)
                       cTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAAc                                                                                                                             >  5:190563/1‑149 (MQ=255)
                       cTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAAc                                                                                                                             >  3:103165/1‑149 (MQ=255)
                        tCGGTGGAAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAAcc                                                                                                                            >  6:106373/1‑149 (MQ=255)
                          ggTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑AAATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTggaccgt                                                                                                                            <  8:273209/147‑7 (MQ=255)
                           gTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCCCGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACcttc                                                                                                                         >  2:155698/1‑149 (MQ=255)
                           gTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGTCGCCAGCCAGGCACTGCAAACcttc                                                                                                                         <  4:57027/149‑1 (MQ=255)
                           gTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACcttc                                                                                                                         <  1:360524/149‑1 (MQ=255)
                             ggCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACcttctt                                                                                                                       >  3:31610/1‑149 (MQ=255)
                                aaaGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCCTGGCTggcggccatg                                                                                                                                       >  1:335953/1‑128 (MQ=255)
                                aaaGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCCTGGCTggcggccatg                                                                                                                                       <  2:335953/130‑3 (MQ=255)
                                aaaGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTccc                                                                                                                    >  3:146998/1‑149 (MQ=255)
                                aaaGCGCCGTGCGCAGAGCACATGGCCAAAGAAATCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGAATACAGAGTCTGGCAGACGAACCAGCCTTCACGGTCCAGAGCATGGCCCCCAGCCAGGAACTGCAAACCTTCTTccc                                                                                                                    >  6:99800/1‑149 (MQ=255)
                                      ccGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCCTGGCTGGCGGCCATGCACTggacc                                                                                                                              >  1:263997/1‑129 (MQ=255)
                                      ccGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCCTGGCTGGCGGCCATGCACTggacc                                                                                                                              <  2:263997/133‑5 (MQ=255)
                                        gTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTg                                                                                                            <  5:106373/149‑1 (MQ=255)
                                          gcgcAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGc                                                                                                          <  6:288948/149‑1 (MQ=255)
                                           cgcAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGc                                                                                                          >  2:31273/1‑148 (MQ=255)
                                             cAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTg                                                                                                            <  5:10155/144‑1 (MQ=255)
                                             cAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTg                                                                                                            >  6:10155/1‑144 (MQ=255)
                                               ttcccacGGGGCAAATTTTTTCTTGGTCCACAGGGGGGTTGCTTTCGGCAATACCTGTATAGCCCATGAGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCaaa                                                                                                       <  1:312195/142‑1 (MQ=255)
                                                       ggCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGGCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGcc                                                                                             <  6:261611/149‑1 (MQ=255)
                                                            aaGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCtt                                                                                        >  7:349381/1‑149 (MQ=255)
                                                                  tCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCa                                                                                 >  8:281905/1‑149 (MQ=255)
                                                                        gTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTgcgc                                                                           >  6:371984/1‑149 (MQ=255)
                                                                        gTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACGGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTgcgc                                                                           >  6:58188/1‑149 (MQ=255)
                                                                              aGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCctggct                                                                     <  4:146998/149‑1 (MQ=255)
                                                                                gcggcTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCctggctgg                                                                   <  8:349381/149‑1 (MQ=255)
                                                                                  ggcTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGt                                                                  >  6:208208/1‑148 (MQ=255)
                                                                                       tCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCgg                                                            >  5:360757/1‑149 (MQ=255)
                                                                                       tCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCgg                                                            <  6:360757/149‑1 (MQ=255)
                                                                                          ttCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGcc                                                         >  3:235699/1‑149 (MQ=255)
                                                                                                 acaTCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTc                                                   >  8:345511/1‑149 (MQ=255)
                                                                                                   aTCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGa                                                 >  2:380513/1‑149 (MQ=255)
                                                                                                    tCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGAc                                                <  4:31610/149‑1 (MQ=255)
                                                                                                     cGGTATACAGCGTCGGGCAGACGAACCAGCCTTCACGGTCCAGTGCATGGCCCCCACCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGAc                                                >  8:252676/1‑148 (MQ=255)
                                                                                                          taCAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTc                                          <  1:380513/149‑1 (MQ=255)
                                                                                                          taCAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTc                                          <  5:111424/149‑1 (MQ=255)
                                                                                                             aGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGt                                       >  7:39653/1‑149 (MQ=255)
                                                                                                              gCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTc                                      <  2:204142/149‑1 (MQ=255)
                                                                                                                 tCGGGCTGACGAACCAGGCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTTCATACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGGCAGCTTC‑AGGCATTGCGCGTGGCTGGTCAGCGGGCCCATGTCGAGCTCTTCGGTCaaa                                   <  5:99800/149‑1 (MQ=255)
                                                                                                                    ggCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGa                                <  6:270149/149‑1 (MQ=255)
                                                                                                                     gCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGAt                               >  2:35648/1‑149 (MQ=255)
                                                                                                                       tGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTa                             >  5:365190/1‑149 (MQ=255)
                                                                                                                        gACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTac                            >  7:30410/1‑149 (MQ=255)
                                                                                                                         aCGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGCCCTCTTCGGTCAAAGGATTacc                           >  8:4195/1‑149 (MQ=255)
                                                                                                                         aCGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTacc                           <  4:103165/149‑1 (MQ=255)
                                                                                                                         aCGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTacc                           <  1:31273/149‑1 (MQ=255)
                                                                                                                             aCCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTacc                           >  7:255113/1‑145 (MQ=255)
                                                                                                                             aCCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTacc                           <  8:255113/145‑1 (MQ=255)
                                                                                                                               cAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCa                                                                                                                                     <  6:306386/38‑1 (MQ=255)
                                                                                                                               cAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCa                                                                                                                                     >  5:306386/1‑38 (MQ=255)
                                                                                                                               cAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCg                     >  2:93522/1‑149 (MQ=255)
                                                                                                                                 gCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCg                     >  7:12354/1‑147 (MQ=255)
                                                                                                                                    ttCACGGTCCAGTGCGTGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCGCAGGt                >  4:106228/1‑149 (MQ=255)
                                                                                                                                    ttCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCGCAGGt                >  4:194816/1‑149 (MQ=255)
                                                                                                                                        cGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCGCAGGTTTTc            <  6:34100/149‑1 (MQ=255)
                                                                                                                                         ggTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCGCAGGTTTTcc           <  3:194816/149‑1 (MQ=255)
                                                                                                                                          gTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCGCAGGTTTTCCa          <  8:39653/149‑1 (MQ=255)
                                                                                                                                              aGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCGCAGGTTTTCCAGTgc      >  4:239444/1‑149 (MQ=255)
                                                                                                                                                tGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCGCAGGTTTTCCAGTgccg    >  5:160402/1‑149 (MQ=255)
                                                                                                                                                  cATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCGCAGGTTTTCCAGTgccgcc  >  6:309237/1‑149 (MQ=255)
                                                                                                                                                  cATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCGCAGGTTTTCCAGTgccgc   <  1:333612/148‑1 (MQ=255)
                                                                                                                                                   aTGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTCTAGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCaaa                                   <  1:283088/116‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTGGCTTCGGCAATGGCCTGTTCCTCGGTGGCAAAGCGCCGTGCGCAGAGCACCGGGCCAAAGATTTCCTCGGTCCACAGGCGGCTGTCTTTCGGC‑ACATCGGTATACAGCGTCGGGCTGACGAACCAGCCTTCACGGTCCAGTGCATGGCCGCCAGCCAGGCACTGCAAACCTTCTTCCCGCGCCGTGGCAAAGTAGCTGGCCACCTTC‑AGCCATTGCGCCTGGCTGGTCAGCGGGCCCATGTCGACCTCTTCGGTCAAAGGATTACCCACCCGCAGGTTTTCCAGTGCCGCC  >  NC_002947/823402‑823695

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