Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,318,912 T→G 6.7% M385R (ATG→AGG)  PP_2039 → acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,318,9120TG 6.7% 229.8 / 6.8 75M385R (ATG→AGG) PP_2039acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (39/31);  new base G (3/2);  total (42/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 = 7.25e-01

GCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCACGGCCTGGATGATTGCGGTACTGCTGGCGCTGTTGAT  >  NC_002947/2318767‑2319058
                                                                                                                                                 |                                                                                                                                                  
gCCGATCGCGCCTTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGc                                                                                                                                                  <  5:258454/148‑1 (MQ=255)
  cggcaccgcCGTGTATCTGCACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCt                                                                                                                                                >  2:258853/7‑148 (MQ=255)
   ggcaccgcCGTGTATCTGCACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCt                                                                                                                                                >  5:313091/6‑147 (MQ=255)
    ttCGCGGCGTGGAGATACACGGTGGTGCCGGGCTGCGGCGTGGCAACCCGGTGGAGTTGATTTTCCTTGCGTCACGTTTTATGCGCATTCTCACCGGATCCAGCGTGATCCAGCTCAAAACCATCGCCAAGGTCGTGCTCATGCCTGAg                                                                                                                                             <  8:97822/148‑1 (MQ=255)
      cgcgCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGc                                                                                                                                            >  1:155478/1‑148 (MQ=255)
       gcgcCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCAGGCATGAGCac                                                                                                                                          <  8:273809/149‑2 (MQ=255)
        cgcCGTGTATCTGCACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAAt                                                                                                                                         >  5:261574/1‑149 (MQ=255)
         gcCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTTCTCCTGAATGAGCAAt                                                                                                                                         >  5:117037/1‑148 (MQ=255)
          ccGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGa                                                                                                                                       <  8:237366/149‑1 (MQ=255)
          ccGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGa                                                                                                                                       <  8:288793/149‑1 (MQ=255)
                      acGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATgccc                                                                                                                           >  5:19928/1‑149 (MQ=255)
                        ggcggTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATgcccgg                                                                                                                         >  8:180768/1‑149 (MQ=255)
                            gTGCCGGGCTGCTGCGGGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATgcccggcccg                                                                                                                     <  6:313091/149‑1 (MQ=255)
                               ccGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCgtgg                                                                                                                  <  8:65425/149‑1 (MQ=255)
                                 gggCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCgtggtg                                                                                                                >  6:271746/1‑149 (MQ=255)
                                     tgctgcGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCg                                                                                                            >  2:304785/1‑149 (MQ=255)
                                       ctgcGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCgg                                                                                                           <  5:271746/148‑1 (MQ=255)
                                        tgcGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGt                                                                                                         >  4:252191/1‑149 (MQ=255)
                                         gcGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGtt                                                                                                        >  2:48521/1‑149 (MQ=255)
                                            tGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTtctc                                                                                                     >  6:60300/1‑149 (MQ=255)
                                            tGGAAACCCGGTGGAGTGGATATTCCGTGCGACACGTAAGATGCGCATTCTAAACGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTCAGCAATGACGCGGTATGCCCGGCCCGTCGTGGCCGGGTtctc                                                                                                     >  4:99747/1‑149 (MQ=255)
                                               cAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTtctc                                                                                                     >  8:75454/1‑146 (MQ=255)
                                                    cGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGAc                                                                                                >  4:201735/1‑146 (MQ=255)
                                                    cGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGAc                                                                                                >  2:190271/1‑146 (MQ=255)
                                                            tggaTATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATaa                                                                                     <  7:180768/149‑1 (MQ=255)
                                                                 atTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACga                                                                                  >  6:325417/1‑147 (MQ=255)
                                                                  ttCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACgagag                                                                               <  2:41671/149‑1 (MQ=255)
                                                                  ttCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACgaga                                                                                <  2:155478/148‑1 (MQ=255)
                                                                      gTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGatca                                                                           <  3:252191/149‑1 (MQ=255)
                                                                          gACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAAt                                                                                       >  6:193185/1‑133 (MQ=255)
                                                                          gACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAAt                                                                                       <  5:193185/133‑1 (MQ=255)
                                                                           acacGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATcatcat                                                                       >  7:121260/1‑148 (MQ=255)
                                                                            cacGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATg                                                                      >  6:230267/1‑148 (MQ=255)
                                                                             acGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGAc                                                                    <  1:304785/149‑1 (MQ=255)
                                                                              cGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACt                                                                   <  8:138153/149‑1 (MQ=255)
                                                                               gTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTc                                                                  <  1:190852/149‑1 (MQ=255)
                                                                                 atgatgCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCa                                                                 >  6:192769/1‑148 (MQ=255)
                                                                                  tgatgCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATa                                                               <  8:121260/149‑1 (MQ=255)
                                                                                    atgCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACt                                                             >  5:190354/1‑149 (MQ=255)
                                                                                      gcgcATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATCACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACtca                                                           >  7:122413/1‑147 (MQ=255)
                                                                                       cgcATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGa                                                           <  8:45994/148‑1 (MQ=255)
                                                                                          aTTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACAcc                                                       >  2:253364/1‑149 (MQ=255)
                                                                                           ttCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGATATCATCATGACTCATACTGACAccc                                                      >  8:103665/1‑149 (MQ=255)
                                                                                              tcACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCg                                                   >  8:186281/1‑149 (MQ=255)
                                                                                               cACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGac                                                          >  2:96108/1‑141 (MQ=255)
                                                                                               cACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGac                                                          <  1:96108/141‑1 (MQ=255)
                                                                                                 ccGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGGGGCCGGGTTCTCCAGACAAGAACAATAACGAGATCTCATCATGACTCATACTGACACCCACGc                                                  >  1:77772/1‑147 (MQ=255)
                                                                                                   ggAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGt                                              <  7:174010/149‑1 (MQ=255)
                                                                                                      agcAGCTAGATCCAGCGGAACTCCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGATCAATAATGAGAGATCATCATGACTCATACTGACACCCCGGCACGTGCa                                           <  8:122413/147‑1 (MQ=255)
                                                                                                      aCCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCa                                           <  1:190271/149‑1 (MQ=255)
                                                                                                       ccAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAg                                          <  1:253364/149‑1 (MQ=255)
                                                                                                           cGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGatca                                                                           <  5:96679/112‑1 (MQ=255)
                                                                                                           cGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGatca                                                                           >  6:96679/1‑112 (MQ=255)
                                                                                                           cGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCAc                                      >  5:27957/1‑149 (MQ=255)
                                                                                                           cGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCAc                                      >  5:59838/1‑149 (MQ=255)
                                                                                                               aTCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACgaga                                                                                <  1:160425/103‑1 (MQ=255)
                                                                                                               aTCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACgaga                                                                                >  2:160425/1‑103 (MQ=255)
                                                                                                                    gcgcAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCACGGCCTGgat                             <  5:60300/149‑1 (MQ=255)
                                                                                                                    gcgcAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCACGGCCTGgat                             <  5:230267/149‑1 (MQ=255)
                                                                                                                       cAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCACGGCCTGgatgat                          >  5:189802/1‑149 (MQ=255)
                                                                                                                       cAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCGCACGTGCAGGCACGGCCTGGATGAtt                         <  1:48521/149‑1 (MQ=255)
                                                                                                                             cATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATAc                                                              >  1:5595/1‑107 (MQ=255)
                                                                                                                             cATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATAc                                                              <  2:5595/107‑1 (MQ=255)
                                                                                                                             cATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCACGGCCTGGATGATTGCgg                     >  1:5242/1‑148 (MQ=255)
                                                                                                                                  ccAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACAccc                                                      <  7:235983/110‑1 (MQ=255)
                                                                                                                                  ccAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACAccc                                                      >  8:235983/1‑110 (MQ=255)
                                                                                                                                     cgtcattGCTCAGGCATGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCACGGCCTGGATGATTGCGGTACTGCTgg            <  2:69717/143‑1 (MQ=255)
                                                                                                                                      gggCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCACTGCCTGGATGATTTCGGTActgatg             >  3:166107/1‑147 (MQ=255)
                                                                                                                                         cGTGCTCATGCCTGAGCAATGACGCGGTATGCACGGCCCGTGGAGGCCGGGTTCACCAGACAACAACAATAACGAGAGATCATCATGACTCATACTGACACCACCGCACGTGCAGGCACGCCCTGGATGATTGCGGTACTGCTGGCGCt        >  4:269854/1‑149 (MQ=255)
                                                                                                                                         attGCTCAGGCATGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCACGGCCTGGATGATTGCGGTACTGCTGGCGCt        >  7:262867/3‑149 (MQ=255)
                                                                                                                                          ttGCTCAGGCATGAGCAATGCCGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCACGGCCTGGATGATTGCGGTACTGCTGGCGCTg       >  6:58129/2‑149 (MQ=255)
                                                                                                                                          gTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCACGGCCTGGATGATTGCGGTACTGCTGGCGCTg       <  5:192769/149‑1 (MQ=255)
                                                                                                                                            gCTCAGGCATGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGCCTCATACTGACACCCACGCACGTGAAGGCACGGCCTGGATGATTGCGGTACTGCTGGCGCTGtt     >  7:80248/1‑149 (MQ=255)
                                                                                                                                               cATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACTTGCAGGCACGGCATGGATGATTGCGGTAATGCTGGCGCTGTTGAt  >  3:314461/1‑149 (MQ=255)
                                                                                                                                                aTGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCACGGCCTGGATGATTGCGGTACTGCTGGCGCTGTTGa   <  7:103665/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGTGGCAACCCGGTGGAGTGGATATTCCGTGCGACACGTATGATGCGCATTCTCACCGGAACCAGCGAGATCCAGCGCAACACCATCGCCAAGGGCGTGCTCATGCCTGAGCAATGACGCGGTATGCCCGGCCCGTGGTGGCCGGGTTCTCCAGACAAGAACAATAACGAGAGATCATCATGACTCATACTGACACCCCCGCACGTGCAGGCACGGCCTGGATGATTGCGGTACTGCTGGCGCTGTTGAT  >  NC_002947/2318767‑2319058

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