Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,133,981 T→A 8.7% Q350L (CAG→CTG)  PP_3638 ← acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,133,9810TA 8.7% 230.8 / 10.3 80Q350L (CAG→CTG) PP_3638acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (35/38);  new base A (4/3);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

ACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGC  >  NC_002947/4133834‑4134126
                                                                                                                                                       |                                                                                                                                                 
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                                                                      >  3:85470/1‑122 (MQ=255)
                                                                                                          aTGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                      <  4:85470/122‑1 (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                                          <  3:182260/145‑1 (MQ=255)
                                                                                                               aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGt                                          <  3:187100/145‑1 (MQ=255)
                                                                                                               aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGt                                          >  4:187100/1‑145 (MQ=255)
                                                                                                               aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGt                                          >  4:182260/1‑145 (MQ=255)
                                                                                                               aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGc                                      >  2:271844/1‑149 (MQ=255)
                                                                                                               aCGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGc                                      >  3:37611/1‑149 (MQ=255)
                                                                                                                        cAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccg                               <  8:322819/147‑1 (MQ=255)
                                                                                                                                ccGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCt                     >  7:184556/1‑149 (MQ=255)
                                                                                                                                ccGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCt                     >  4:134494/1‑149 (MQ=255)
                                                                                                                                 cgcgCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCTGGCAAACCCTTCCCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTg                    <  2:311912/149‑1 (MQ=255)
                                                                                                                                 cgcgCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTg                    >  2:64485/1‑149 (MQ=255)
                                                                                                                                 cgcgCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTg                    <  2:257709/149‑1 (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                                            <  5:273463/117‑1 (MQ=255)
                                                                                                                                         tGCA‑ACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCagag                                            >  6:273463/1‑117 (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           <  8:181520/142‑1 (MQ=255)
                                                                                                                                                  ctgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg           >  7:181520/1‑142 (MQ=255)
                                                                                                                                                  ctgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTAcc    >  2:177143/1‑149 (MQ=255)
                                                                                                                                                     cctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACcgc  <  1:315733/148‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
ACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGC‑GCCAGCACCT‑GCTGCCAACCACC‑CTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCA‑ACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGC  >  NC_002947/4133834‑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: ‑