Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,533,413 A→G 7.0% E16E (GAA→GAG atoA → 3‑oxoacid CoA‑transferase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,533,4130AG 7.0% 246.8 / 6.6 86E16E (GAA→GAGatoA3‑oxoacid CoA‑transferase subunit A
Reads supporting (aligned to +/- strand):  ref base A (42/38);  new base G (3/3);  total (45/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.56e-01

TGACGAAAACAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGTTTTGGTCTGG  >  NC_002947/3533265‑3533558
                                                                                                                                                    |                                                                                                                                                 
tGACGAAAACAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAg                                                                                                                                                  <  2:1555371/150‑1 (MQ=255)
  aCGAAAACAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGc                                                                                                                                                <  2:1579446/150‑1 (MQ=255)
    gAAAACAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCTTCGAGGGCCt                                                                                                                                              >  2:299796/1‑150 (MQ=255)
    gAAAACAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCt                                                                                                                                              >  1:335809/1‑150 (MQ=255)
    gAAAACAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCt                                                                                                                                              <  2:2291064/150‑1 (MQ=255)
      aaaCAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCTTCGAGGGCCTgc                                                                                                                                            <  1:299796/150‑2 (MQ=255)
       aaCAAAATTAACAAATCAGTCGCTTGCTGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGAc                                                                                                                                           <  1:363693/150‑1 (MQ=255)
       aaCAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGAc                                                                                                                                           <  1:449672/150‑1 (MQ=255)
           aaaTTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACGACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGac                                                                                                                                       <  1:2221634/150‑1 (MQ=255)
                 aCAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACAt                                                                                                                                  <  2:2049577/149‑1 (MQ=255)
                 aCAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACAt                                                                                                                                  >  1:2049577/1‑149 (MQ=255)
                  cAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGa                                                                                                                                >  2:857699/1‑150 (MQ=255)
                   aaaTCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGa                                                                                                                                            >  1:639050/1‑137 (MQ=255)
                   aaaTCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGa                                                                                                                                            <  2:639050/137‑1 (MQ=255)
                   aaaTCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGAc                                                                                                                               <  1:544858/150‑1 (MQ=255)
                    aaTCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGCCAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACg                                                                                                                              >  1:423451/1‑150 (MQ=255)
                          tCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACCGCCTCGTCGGTTAAACAACTACAAAAGGCCGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGCCCTGACCGACAACATGACGGTACTg                                                                                                                        >  2:655157/1‑150 (MQ=255)
                            gCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGc                                                                                                                      >  1:216964/1‑150 (MQ=255)
                            gCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGc                                                                                                                      <  2:963101/150‑1 (MQ=255)
                              ttGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCTTCGAGGGTCTgc                                                                                                                                            >  1:1193/1‑125 (MQ=255)
                              ttGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCTTCGAGGGCCTgc                                                                                                                                            <  2:1193/126‑2 (MQ=255)
                              ttGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCg                                                                                                                    <  1:905539/150‑1 (MQ=255)
                              ttGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCg                                                                                                                    >  1:1569890/1‑150 (MQ=255)
                                 cGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCggtg                                                                                                                 <  1:857699/150‑1 (MQ=255)
                                   gTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGcc                                                                                                                                               <  1:1686921/118‑1 (MQ=255)
                                   gTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGcc                                                                                                                                               >  2:1686921/1‑118 (MQ=255)
                                   gTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCggtggt                                                                                                               <  2:216964/150‑1 (MQ=255)
                                       aGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAg                                                                                                                                                  >  2:156063/1‑111 (MQ=255)
                                       aGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAg                                                                                                                                                  <  1:156063/111‑1 (MQ=255)
                                         aTGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGc                                                                                                         >  1:1198029/1‑150 (MQ=255)
                                               gCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGc                                                                                                   <  1:1615158/150‑1 (MQ=255)
                                                     gCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGc                                                                                                                                                <  2:622878/99‑1 (MQ=255)
                                                     gCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGc                                                                                                                                                >  1:622878/1‑99 (MQ=255)
                                                       tCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATccc                                                                                           >  2:1781726/1‑150 (MQ=255)
                                                             aCAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGaaaa                                                                                     >  1:38246/1‑150 (MQ=255)
                                                                   tcgtcgGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCGCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTAt                                                                               <  1:1837639/150‑1 (MQ=255)
                                                                       cgGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGc                                                                           <  2:1198029/150‑1 (MQ=255)
                                                                        ggTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGAc                                                                                                                               <  1:829681/97‑1 (MQ=255)
                                                                        ggTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGAc                                                                                                                               >  2:829681/1‑97 (MQ=255)
                                                                        ggTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCg                                                                          <  2:335809/150‑1 (MQ=255)
                                                                          ttAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGaa                                                                        >  2:379118/1‑150 (MQ=255)
                                                                              acaacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCTTCGAGGGCCTgct                                                                                                                                           >  1:792064/1‑77 (MQ=255)
                                                                              acaacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCTTCGAGGGCCTgct                                                                                                                                           <  2:792064/79‑3 (MQ=255)
                                                                              acaacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCa                                                                    >  1:56081/1‑150 (MQ=255)
                                                                                aacTACAAAAGGCGGATCAGTATGGCCGGACGGGACGGGCGTGTTGCAAACAATGAACAGGCGCTCTAAAGCGTGAGCGACAACATGACGGTACAGGTCGGTGGTTTCGGCGTGTCCGTCATCCCGGAAAACCTTATCAGCGAAATCAAg                                                                  >  2:658188/1‑150 (MQ=255)
                                                                                aacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCg                                                                                                  >  2:1578513/1‑118 (MQ=255)
                                                                                aacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCg                                                                                                  <  1:1578513/118‑1 (MQ=255)
                                                                                aacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAg                                                                  <  2:1606880/150‑1 (MQ=255)
                                                                                aacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAg                                                                  <  1:1781726/150‑1 (MQ=255)
                                                                                aacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAg                                                                  >  1:1911040/1‑150 (MQ=255)
                                                                                   tACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAgcgg                                                               >  1:299060/1‑150 (MQ=255)
                                                                                   tACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAgcgg                                                               >  1:602271/1‑150 (MQ=255)
                                                                                      aaaaGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAg                                                                                                                                                  <  1:1660381/64‑1 (MQ=255)
                                                                                      aaaaGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAg                                                                                                                                                  >  2:1660381/1‑64 (MQ=255)
                                                                                       aaaGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGa                                                                                                                                            >  2:985351/1‑69 (MQ=255)
                                                                                       aaaGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGa                                                                                                                                            <  1:985351/69‑1 (MQ=255)
                                                                                        aaGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGCAAACCTTATCAGCGAAAGCAAgcggcgcgg                                                          >  2:21333/1‑150 (MQ=255)
                                                                                                  aGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGc                                                <  1:379118/150‑1 (MQ=255)
                                                                                                          cGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGt                                        >  1:1822312/1‑150 (MQ=255)
                                                                                                          cGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGt                                        >  1:912438/1‑150 (MQ=255)
                                                                                                          cGGACTGGACAAGCGCATTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATc                                                                                             >  1:2064326/1‑97 (MQ=255)
                                                                                                          cGGACTGGACAAGCGCATTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATc                                                                                             <  2:2064326/97‑1 (MQ=255)
                                                                                                            gACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTgg                                      <  2:1569890/150‑1 (MQ=255)
                                                                                                             aCTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGt                                     >  1:1807796/1‑150 (MQ=255)
                                                                                                                ggACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTc                                  <  2:38246/150‑1 (MQ=255)
                                                                                                                ggACAAGCGCGTTGAAACGTATCATCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTc                                  <  2:1159912/150‑1 (MQ=255)
                                                                                                                 gACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCg                                                                                                  <  1:1464392/85‑1 (MQ=255)
                                                                                                                 gACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCg                                                                                                  >  2:1464392/1‑85 (MQ=255)
                                                                                                                    aaGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCggtg                                                                                                                 <  1:516549/67‑1 (MQ=255)
                                                                                                                    aaGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCggtg                                                                                                                 >  2:516549/1‑67 (MQ=255)
                                                                                                                    aaGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATc                                                                     >  2:1963342/1‑111 (MQ=255)
                                                                                                                    aaGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATc                                                                     <  1:1963342/111‑1 (MQ=255)
                                                                                                                    aaGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATAAAGCGGCGGGGCGTCAAGGGCCTGACCGTGGTTTCcaac                              >  2:802086/1‑150 (MQ=255)
                                                                                                                          gTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTgc                        >  2:2076757/1‑150 (MQ=255)
                                                                                                                           ttGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTgcg                       <  2:1822312/150‑1 (MQ=255)
                                                                                                                               acccTATGAGAAGGCCCTTGAAGGCCTGACCGACAACATGACGGTACTTGCCGGTGGTTTCGGCCTGTGCGGCATCCCGAAAAACCTTATCAGTGAAATCAAGCGGCGCGGCGTCAAGGGACTGACCGTGGTTTCCAATAACTGTGGCGt                   <  2:78890/148‑1 (MQ=255)
                                                                                                                               aaCCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGt                   >  1:469041/1‑150 (MQ=255)
                                                                                                                                gccgaTGAGCGGGCGCTCGGAGGGGTGACGGACAATATGACGGTAGTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACAGTGGTTTCCAACAACTGCGGCGTg                  <  2:423451/146‑1 (MQ=255)
                                                                                                                                  cTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGa                >  2:1779004/1‑150 (MQ=255)
                                                                                                                                  cTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGa                >  2:2006592/1‑150 (MQ=255)
                                                                                                                                   tATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGAc               >  1:341754/1‑150 (MQ=255)
                                                                                                                                       aGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGtt           >  2:496644/1‑150 (MQ=255)
                                                                                                                                         cAGGCCCTCGAAGGCCTGACCGAACACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGtttt         >  1:933472/1‑150 (MQ=255)
                                                                                                                                             cccTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAAc                                                                                    <  1:145213/71‑1 (MQ=255)
                                                                                                                                             cccTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAAc                                                                                    >  2:145213/1‑71 (MQ=255)
                                                                                                                                                tCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGTTTTGGTCTgg  >  2:1494820/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGACGAAAACAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGTTTTGGTCTGG  >  NC_002947/3533265‑3533558

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