Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,133,979 C→G 9.8% A351P (GCA→CCA) ‡ PP_3638 ← acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,133,9790CG 9.8% 233.4 / 13.0 80A351P (GCA→CCA) ‡PP_3638acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (35/37);  new base G (4/4);  total (39/41)
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.57e-01

CACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGC  >  NC_002947/4133831‑4134126
                                                                                                                                                        |                                                                                                                                                   
cacaGCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACgcctgc                                                                                                                                                     <  7:270866/149‑1 (MQ=255)
   accacTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGcctgcctg                                                                                                                                                  >  6:108889/1‑149 (MQ=255)
        tGGGGCCAAGGATGCTCACGACCTCACCAGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGtt                                                                                                                                             >  1:357863/1‑149 (MQ=255)
          gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGtt                                                                                                                                             <  5:137702/147‑1 (MQ=255)
          gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTg                                                                                                                                            <  8:71837/148‑1 (MQ=255)
          gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTg                                                                                                                                            >  2:129713/1‑148 (MQ=255)
          gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTg                                                                                                                                            <  8:282809/148‑1 (MQ=255)
          gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTg                                                                                                                                            <  1:129713/148‑1 (MQ=255)
          gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGc                                                                                                                                           <  7:217378/149‑1 (MQ=255)
            gCCAAGGCAGCACAAGACCTCCAGTGGTGCCAGCTTCATGTACAACTGC‑GCCATCACCTATCTGCAAACCCCCTATTACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCAGGCAGGCGTTGc                                                                                                                                           <  1:290042/149‑1 (MQ=255)
                   aTGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCACACGCCTGCCTGGCGTTGCAGTTCGgc                                                                                                                                   <  1:186050/149‑1 (MQ=255)
                       tCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGgcgcgc                                                                                                                               >  2:116557/1‑148 (MQ=255)
                                   cTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGc                                                                                                                  >  6:354165/1‑149 (MQ=255)
                                      gTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTgg                                                                                                               <  1:116557/149‑1 (MQ=255)
                                         ccAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGc                                                                                                             <  5:108889/148‑1 (MQ=255)
                                         ccAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGc                                                                                                             <  5:354165/148‑1 (MQ=255)
                                          ctggcgCAGTTCGACCTGCAGCTGGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCg                                                                                                            >  7:90219/6‑149 (MQ=255)
                                                   gTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTggg                                                                                                              >  2:310770/1‑137 (MQ=255)
                                                   gTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTggg                                                                                                              <  1:310770/137‑1 (MQ=255)
                                                     cGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                         >  2:171015/1‑140 (MQ=255)
                                                     cGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                         <  1:171015/140‑1 (MQ=255)
                                                      gAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCa                                                                                                <  2:357863/148‑1 (MQ=255)
                                                         cTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCAGGCAGGCGTTGCat                                                                                                                                         <  6:86979/104‑2 (MQ=255)
                                                          tGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCGATTGCACCAGGCTGGGCGTGACGATCGGCACGAAg                                                                                           >  5:38858/1‑149 (MQ=255)
                                                            c‑gcCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGCCGATCGTCACgccca                                                                                          <  2:70835/148‑5 (MQ=255)
                                                            c‑gcCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGc                                                                                         <  5:61785/149‑1 (MQ=255)
                                                               ccAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGcc                                                                                        >  4:180898/1‑148 (MQ=255)
                                                                cAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGa                                                                                      >  1:275276/1‑149 (MQ=255)
                                                                   cACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTc                                                                                   <  5:222205/149‑1 (MQ=255)
                                                                          ctgcCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCGATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGc                                                                               <  6:38858/147‑1 (MQ=255)
                                                                               aaCCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggc                                                                        >  2:41285/1‑149 (MQ=255)
                                                                               aaCCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggc                                                                        >  7:56937/1‑149 (MQ=255)
                                                                                 ccacc‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                      <  2:204768/149‑1 (MQ=255)
                                                                                  cacc‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCTACCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCggg                                                                     <  6:129059/149‑1 (MQ=255)
                                                                                  cacc‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCggg                                                                     >  5:147822/1‑149 (MQ=255)
                                                                                   acc‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                         <  8:242766/112‑1 (MQ=255)
                                                                                    cc‑cTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGc                                                                    >  4:44062/1‑148 (MQ=255)
                                                                                    cc‑cTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACACCGCTCATGCA‑ACGCCGGCCTGGCGGTGCAGGTCGGCGCGCAATTGAACCAGGCTGGCCTTGACGATCGGCACGACGGCCGACTCCCGCAATCGGCGGGc                                                                    >  6:356380/1‑148 (MQ=255)
                                                                                     c‑cTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCa                                                                   >  4:83177/1‑148 (MQ=255)
                                                                                        tCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                                >  1:311912/1‑149 (MQ=255)
                                                                                          aCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATc                                                                                                    <  2:298753/111‑1 (MQ=255)
                                                                                          aCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATc                                                                                                    >  1:298753/1‑111 (MQ=255)
                                                                                               ggTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGc                                                         <  1:108954/149‑1 (MQ=255)
                                                                                                gTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGcc                                                        <  7:216242/149‑1 (MQ=255)
                                                                                                 tAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCg                                                       >  1:257709/1‑149 (MQ=255)
                                                                                                    cccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                                <  5:5062/137‑1 (MQ=255)
                                                                                                    cccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                                >  6:5062/1‑137 (MQ=255)
                                                                                                     ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                                >  4:114962/1‑136 (MQ=255)
                                                                                                     ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                                <  3:114962/136‑1 (MQ=255)
                                                                                                     ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAtt                                                   >  7:322819/1‑149 (MQ=255)
                                                                                                        aGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCACCATCGGCGGGCAAACCCTTCGCCGTATTCGc                                                <  2:275276/149‑1 (MQ=255)
                                                                                                             aTGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                      <  4:85470/122‑1 (MQ=255)
                                                                                                             aTGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                      >  3:85470/1‑122 (MQ=255)
                                                                                                                tCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGt                                          >  6:286156/1‑147 (MQ=255)
                                                                                                                tCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGt                                          >  1:106476/1‑147 (MQ=255)
                                                                                                                tCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTaa                                        <  6:147822/149‑1 (MQ=255)
                                                                                                                  aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAAtt                                                                                                                           <  5:243821/64‑1 (MQ=255)
                                                                                                                  aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAAtt                                                                                                                           >  6:243821/1‑64 (MQ=255)
                                                                                                                  aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGt                                          >  4:182260/1‑145 (MQ=255)
                                                                                                                  aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGt                                          <  3:182260/145‑1 (MQ=255)
                                                                                                                  aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGt                                          <  3:187100/145‑1 (MQ=255)
                                                                                                                  aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGt                                          >  4:187100/1‑145 (MQ=255)
                                                                                                                  aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGc                                      >  3:37611/1‑149 (MQ=255)
                                                                                                                  aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGc                                      >  2:271844/1‑149 (MQ=255)
                                                                                                                           cAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccg                               <  8:322819/147‑1 (MQ=255)
                                                                                                                                   ccGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCt                     >  4:134494/1‑149 (MQ=255)
                                                                                                                                   ccGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCt                     >  7:184556/1‑149 (MQ=255)
                                                                                                                                    cgcgCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCTGGCAAACCCTTCCCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTg                    <  2:311912/149‑1 (MQ=255)
                                                                                                                                    cgcgCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTg                    <  2:257709/149‑1 (MQ=255)
                                                                                                                                    cgcgCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTg                    >  2:64485/1‑149 (MQ=255)
                                                                                                                                     gcgcTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTg                    >  4:152501/1‑148 (MQ=255)
                                                                                                                                     gcgcTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGc                   <  3:44062/149‑1 (MQ=255)
                                                                                                                                       gcTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCaa                 <  1:41285/149‑1 (MQ=255)
                                                                                                                                       cgaactGCA‑ACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCaa                 >  8:39904/6‑149 (MQ=255)
                                                                                                                                         tCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACt               >  2:315733/1‑149 (MQ=255)
                                                                                                                                          actGCA‑ACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc              >  2:323993/3‑149 (MQ=255)
                                                                                                                                           aTGGA‑ACTCGTGGCTGGCGTTGCAGTTCGGCGCGCAATTGCACCTGGGTGGGCGTGACTATCGGCACGAAGGCCGACTGGCGCGATCGGCGGGCTAACCCTTCGCCGTATTCGCAGAGGTATTCCTTGCCGCCACTGGCCTGCAACTc              <  2:106476/148‑1 (MQ=255)
                                                                                                                                            tgaa‑aCGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGTGCGTGCCGATCGGCACGAAGGCCGCATCCCGCCATCGGCGGGCAAACCCTTCGCCGTATGCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCa            <  5:356380/146‑1 (MQ=255)
                                                                                                                                            tGCA‑ACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCagag                                            >  6:273463/1‑117 (MQ=255)
                                                                                                                                            tGCA‑ACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCagag                                            <  5:273463/117‑1 (MQ=255)
                                                                                                                                            tGCA‑ACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc              >  5:248056/1‑147 (MQ=255)
                                                                                                                                              cA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg           <  4:37611/148‑1 (MQ=255)
                                                                                                                                                    cctgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTAc     <  1:64485/149‑1 (MQ=255)
                                                                                                                                                     ctgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg           >  7:181520/1‑142 (MQ=255)
                                                                                                                                                     ctgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg           <  8:181520/142‑1 (MQ=255)
                                                                                                                                                     ctgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTAcc    >  2:177143/1‑149 (MQ=255)
                                                                                                                                                        cctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACcgc  <  1:315733/148‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
CACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGC  >  NC_002947/4133831‑4134126

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