Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,457,102 T→G 5.8% E77A (GAG→GCG)  hbdH ← 3‑hydroxybutyrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,457,1020TG 5.8% 196.6 / 4.6 69E77A (GAG→GCG) hbdH3‑hydroxybutyrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (34/31);  new base G (2/2);  total (36/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 = 6.33e-01

CATGCCCGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCC  >  NC_002947/3456954‑3457247
                                                                                                                                                    |                                                                                                                                                 
cATGCCCGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTc                                                                                                                                                  <  1:1257829/150‑1 (MQ=255)
                   ggcgTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAAc                                                                                                                               <  2:325677/149‑1 (MQ=255)
                     cgTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCa                                                                                                                                       >  1:646679/1‑140 (MQ=255)
                     cgTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCa                                                                                                                                       <  2:646679/140‑1 (MQ=255)
                         ccATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGcc                                                                                                                         <  1:142377/150‑1 (MQ=255)
                           aTGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCt                                                                                                                        >  2:130831/1‑149 (MQ=255)
                           aTGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCt                                                                                                                        <  1:130831/149‑1 (MQ=255)
                           aTGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTc                                                                                                                       >  2:58293/1‑150 (MQ=255)
                             ggAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGa                                                                                                                     <  1:421434/150‑1 (MQ=255)
                                aTACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACGCGAGTTCgg                                                                                                                                          >  1:493899/1‑126 (MQ=255)
                                aTACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACGCGAGTTCgg                                                                                                                                          <  2:493899/126‑1 (MQ=255)
                                  aCGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTgg                                                                                                                <  2:172069/150‑1 (MQ=255)
                                   cGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGa                                                                                                               <  1:58293/150‑1 (MQ=255)
                                   cGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGa                                                                                                               <  2:88576/150‑1 (MQ=255)
                                    ggCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGtt                                                                                                                                   >  1:80615/1‑129 (MQ=255)
                                    ggCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGtt                                                                                                                                   <  2:80615/129‑1 (MQ=255)
                                    ggCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGAc                                                                                                              <  1:991434/150‑1 (MQ=255)
                                         aCAGGTTCATGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACAt                                                                                                         >  1:282342/1‑150 (MQ=255)
                                          cAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATc                                                                                                        <  1:1220855/150‑1 (MQ=255)
                                           aGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATcc                                                                                                       >  1:229285/1‑150 (MQ=255)
                                             gTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGa                                                                                                     >  2:392163/1‑150 (MQ=255)
                                              ttCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGAc                                                                                                    >  1:820104/1‑150 (MQ=255)
                                                     cGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGCCCAGGATGTCGACGCCGCCGCACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGa                                                                                                               <  1:729623/132‑1 (MQ=255)
                                                     cGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGa                                                                                                               >  2:729623/1‑132 (MQ=255)
                                                      gatgatCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGc                                                                                            <  1:840445/150‑1 (MQ=255)
                                                      gatgatCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGc                                                                                            >  2:1142289/1‑150 (MQ=255)
                                                                tCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGa                                                                                                                     >  2:1004493/1‑115 (MQ=255)
                                                                tCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGa                                                                                                                     <  1:1004493/115‑1 (MQ=255)
                                                                    aGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACaa                                                                              >  2:962491/1‑150 (MQ=255)
                                                                             gCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCAccc                                                                     >  2:196965/1‑150 (MQ=255)
                                                                                     tGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCgg                                                             <  1:196965/150‑1 (MQ=255)
                                                                                     tGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACACGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCgg                                                             <  1:1032369/150‑1 (MQ=255)
                                                                                       cTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGc                                                           >  2:273904/1‑150 (MQ=255)
                                                                                       cTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGc                                                           >  1:951224/1‑150 (MQ=255)
                                                                                            ccGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAAttt                                                      >  2:110440/1‑150 (MQ=255)
                                                                                                tGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGc                                                  <  1:490860/150‑1 (MQ=255)
                                                                                                    acaTGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGTTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATg                                              >  1:1255900/1‑150 (MQ=255)
                                                                                                    acaTGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATg                                              >  2:884145/1‑150 (MQ=255)
                                                                                                        gCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCgg                                          <  1:392163/150‑1 (MQ=255)
                                                                                                        gCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCgg                                          <  2:820104/150‑1 (MQ=255)
                                                                                                            gATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGcc                                      >  2:90558/1‑150 (MQ=255)
                                                                                                            gATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGcc                                      >  1:193495/1‑150 (MQ=255)
                                                                                                               aCCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCgg                                   <  2:951224/150‑1 (MQ=255)
                                                                                                                 cGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCgg                                          >  1:1275628/1‑141 (MQ=255)
                                                                                                                 cGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCgg                                          <  2:1275628/141‑1 (MQ=255)
                                                                                                                 cGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGt                                 >  1:513821/1‑150 (MQ=255)
                                                                                                                  ggCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCgg                                   >  1:813200/1‑147 (MQ=255)
                                                                                                                  ggCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCgg                                   <  2:813200/147‑1 (MQ=255)
                                                                                                                    cGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCgg                                   <  1:553794/145‑1 (MQ=255)
                                                                                                                    cGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCgg                                   >  2:553794/1‑145 (MQ=255)
                                                                                                                       tgttgACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCggg                                         <  1:541194/136‑1 (MQ=255)
                                                                                                                       tgttgACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCggg                                         >  2:541194/1‑136 (MQ=255)
                                                                                                                         ttgACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGa                                                                                                     <  2:239528/74‑1 (MQ=255)
                                                                                                                         ttgACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGa                                                                                                     >  1:239528/1‑74 (MQ=255)
                                                                                                                           gACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGcc                       >  1:86460/1‑150 (MQ=255)
                                                                                                                             ccAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCg                                           <  1:174200/128‑1 (MQ=255)
                                                                                                                             ccAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCg                                           >  2:174200/1‑128 (MQ=255)
                                                                                                                              cAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGtt                    <  1:1142289/150‑1 (MQ=255)
                                                                                                                              cAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGtt                    <  1:110440/150‑1 (MQ=255)
                                                                                                                              cAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGtt                    >  2:518614/1‑150 (MQ=255)
                                                                                                                              cAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGtt                    <  1:962491/150‑1 (MQ=255)
                                                                                                                                 gATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCTTTGAg                 >  2:200279/1‑150 (MQ=255)
                                                                                                                                       gACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGAt           >  2:705368/1‑150 (MQ=255)
                                                                                                                                       gACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCCGGTCACCAAAGCCGTTGAGCACGAt           >  2:726975/1‑150 (MQ=255)
                                                                                                                                          ctggccgAACGCGAGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGtt        >  1:1110627/4‑150 (MQ=255)
                                                                                                                                            cgccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCgg                                   <  2:868232/121‑1 (MQ=255)
                                                                                                                                            cgccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCgg                                   >  1:868232/1‑121 (MQ=255)
                                                                                                                                             gccgAACGCGAGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGgc     <  2:1110627/150‑1 (MQ=255)
                                                                                                                                                gAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGgcgca  >  1:83795/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CATGCCCGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCC  >  NC_002947/3456954‑3457247

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: