Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 387,805 C→A 10.5% D315E (GAC→GAA) ‡ glyA‑I → serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947387,8050CA10.5% 168.8 / 9.5 57D315E (GAC→GAA) ‡glyA‑Iserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base A (3/3);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

TGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCGGGCCTGCGTATC  >  NC_002947/387657‑387952
                                                                                                                                                    |                                                                                                                                                   
tGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTTGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGAc                                                                                                                                                     <  1:243404/149‑1 (MQ=255)
    cGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAac                                                                                                                                                  >  4:191705/1‑148 (MQ=255)
         tCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTGTCGGTGCCACCGGAAaccacatc                                                                                                                                            <  1:88765/149‑4 (MQ=255)
            cTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTc                                                                                                                                         <  3:245486/149‑1 (MQ=255)
               ccAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTg                                                                                                                                      <  3:162802/149‑1 (MQ=255)
                cAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTg                                                                                                                                      <  1:3805/148‑1 (MQ=255)
                cAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTg                                                                                                                                      >  2:3805/1‑148 (MQ=255)
                   ggCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGAc                                                                                                                                                     <  4:173080/130‑1 (MQ=255)
                   ggCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGAc                                                                                                                                                     <  4:178757/130‑1 (MQ=255)
                   ggCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGAc                                                                                                                                                     >  3:178757/1‑130 (MQ=255)
                   ggCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGAc                                                                                                                                                     >  3:173080/1‑130 (MQ=255)
                   ggCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCa                                                                                                                                  <  6:312384/149‑1 (MQ=255)
                             ttCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTGTCGGTGCCACCGGAAaccacatc                                                                                                                                            <  4:247067/129‑4 (MQ=255)
                             ttCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTGTCGGTGCCACCGGAAaccacatc                                                                                                                                            >  3:247067/1‑126 (MQ=255)
                               cAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGt                                                                                                                       >  6:123136/1‑148 (MQ=255)
                                aaGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGt                                                                                                                       >  6:129923/1‑147 (MQ=255)
                                  ggAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAgg                                                                                                                   >  7:195867/1‑149 (MQ=255)
                                   gAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAggg                                                                                                                  >  4:295062/1‑149 (MQ=255)
                                   gAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAggg                                                                                                                  <  5:286175/149‑1 (MQ=255)
                                   gAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCAAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAggg                                                                                                                  <  8:182250/149‑1 (MQ=255)
                                      gcgcTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCt                                                                                                               >  2:96428/1‑149 (MQ=255)
                                        gcTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCt                                                                                                                                       <  4:181719/123‑1 (MQ=255)
                                        gcTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCt                                                                                                                                       >  3:181719/1‑123 (MQ=255)
                                                cTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCaaa                                                                                                     <  3:365995/149‑1 (MQ=255)
                                                cTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCaaa                                                                                                     <  3:191705/149‑1 (MQ=255)
                                                  ggATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAga                                                                                                   >  5:211260/1‑149 (MQ=255)
                                                   gATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAg                                                                                                    <  8:195867/147‑1 (MQ=255)
                                                        aaGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgacgccga                                                                                             >  1:75103/1‑149 (MQ=255)
                                                                  agcagcAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTg                                                                                   <  5:123136/149‑1 (MQ=255)
                                                                          gTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTTTGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTgc                                                                           <  1:123230/149‑1 (MQ=255)
                                                                          gTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTgc                                                                           <  2:75103/149‑1 (MQ=255)
                                                                          gTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTgc                                                                           <  5:129923/149‑1 (MQ=255)
                                                                            gATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgacgcc                                                                                               >  5:199059/1‑127 (MQ=255)
                                                                            gATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGGCAGGGCCTGACTGGCAAAgacgcc                                                                                               <  6:199059/127‑1 (MQ=255)
                                                                             aTCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTgcgc                                                                         <  6:211260/148‑1 (MQ=255)
                                                                             aTCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCa                                                                        >  2:322585/1‑149 (MQ=255)
                                                                                    aCGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGTAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCAcc                                                                 <  3:295062/149‑1 (MQ=255)
                                                                                      gCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGt                                                               >  7:52011/1‑149 (MQ=255)
                                                                                        ccAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTc                                                              >  1:203499/1‑148 (MQ=255)
                                                                                          aGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTcaac                                                           <  8:29172/149‑1 (MQ=255)
                                                                                              cATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAg                                                        >  8:124489/1‑148 (MQ=255)
                                                                                                           ttCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGaa                                          >  2:23149/1‑149 (MQ=255)
                                                                                                            tattCGACCGTGGCTAGTTTGGGGTTTCGGGTGGCACGGACAACCACCTGTTCCTGTTCTGCTTGATCCGTCAGGGCCTGACTGGAAAAGAGGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGaa                                          <  1:96428/145‑1 (MQ=255)
                                                                                                              aTCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGa                                       >  8:211291/1‑149 (MQ=255)
                                                                                                              aTCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGa                                       >  8:174443/1‑149 (MQ=255)
                                                                                                                 gACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGAccc                                    >  2:151044/1‑149 (MQ=255)
                                                                                                                          tACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGcc                           >  4:130633/1‑149 (MQ=255)
                                                                                                                            cGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGt                         <  1:144974/149‑1 (MQ=255)
                                                                                                                              aTGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTc                       <  4:361930/149‑1 (MQ=255)
                                                                                                                               ggtgGTTGTCGGTGCCACCGGAAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCg                      >  6:50391/2‑149 (MQ=255)
                                                                                                                                gtgGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGt                     >  6:258201/1‑149 (MQ=255)
                                                                                                                                gtgGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGt                     >  2:176978/1‑149 (MQ=255)
                                                                                                                                gtgGTTGTCGGTGCCACCGGAAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGt                     >  5:12662/1‑149 (MQ=255)
                                                                                                                                        cGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCgg             <  1:23149/149‑1 (MQ=255)
                                                                                                                                          gTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCgg             <  7:314741/147‑1 (MQ=255)
                                                                                                                                          gTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCgg             >  8:314741/1‑147 (MQ=255)
                                                                                                                                            gccACCGGAAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCGGGCCt         <  6:12662/147‑1 (MQ=255)
                                                                                                                                              cACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCGGGCCTGc       <  7:124489/149‑1 (MQ=255)
                                                                                                                                               aCCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGAc                                                                                                            >  6:52035/1‑47 (MQ=255)
                                                                                                                                                 cGACAACCACCTGTTCCCGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGAc                                      >  7:30632/1‑115 (MQ=255)
                                                                                                                                                 cGACAACCACCTGTTCCCGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGAc                                      <  8:30632/115‑1 (MQ=255)
                                                                                                                                                   aCAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCGGGCCTGCGTATc  >  7:124832/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCGGGCCTGCGTATC  >  NC_002947/387657‑387952

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