Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,978,708 C→G 5.9% F45L (TTC→TTG)  cmk → cytidylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,978,7080CG 5.9% 210.9 / 2.1 65F45L (TTC→TTG) cmkcytidylate kinase
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base G (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCC  >  NC_002947/1978562‑1978853
                                                                                                                                                  |                                                                                                                                                 
agagGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCaa                                                                                                                                                 <  7:116081/149‑1 (MQ=255)
 gagGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGCTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAAt                                                                                                                                                <  1:113678/149‑1 (MQ=255)
    ggCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAAt                                                                                                                                                <  5:371972/146‑1 (MQ=255)
    ggCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAAt                                                                                                                                                >  6:371972/1‑146 (MQ=255)
     gCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGc                                                                                                                                              >  7:254715/1‑147 (MQ=255)
     gCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCa                                                                                                                                            <  3:305143/149‑1 (MQ=255)
       aaGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAg                                                                                                                                           >  8:344857/1‑148 (MQ=255)
         gTCGTTAATTCAATGGCACCGGTTCTGACCATCGACGGACCATGCGGCTCGGGCAAGGGAACGGTCGGCGTGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaa                                                                                                                                        <  3:250512/149‑1 (MQ=255)
         gTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaa                                                                                                                                        <  8:254715/149‑1 (MQ=255)
            gTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAAcc                                                                                                                                      <  2:60512/148‑1 (MQ=255)
            gTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCa                                                                                                                                     >  1:227927/1‑149 (MQ=255)
                 ttCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCg                                                                                                                                >  6:237423/1‑149 (MQ=255)
                 ttCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCg                                                                                                                                >  5:299163/1‑149 (MQ=255)
                    aaTGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTc                                                                                                                              >  1:223328/1‑148 (MQ=255)
                      tGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGAc                                                                                                                           <  2:223328/149‑1 (MQ=255)
                        gCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGAcc                                                                                                                          <  6:118679/148‑1 (MQ=255)
                            cGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGAc                                                                                                                      >  6:333625/1‑148 (MQ=255)
                                 aTCACCATCGACGGACCAAGCGGCTCTGGCAAGGGCACGGTCGCCGGGATGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGCAATCCGGCGCGCTGAATAGCTTGCTGGCATAGAATCCCAGCACCCACGGCGTCGACCTGACCAACg                                                                                                                 >  4:52479/1‑148 (MQ=255)
                                  tCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACg                                                                                                                 >  7:241574/1‑147 (MQ=255)
                                     ccATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCAc                                                                                                                                    >  6:64416/1‑125 (MQ=255)
                                     ccATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCAc                                                                                                                                    <  5:64416/125‑1 (MQ=255)
                                        tCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgctg                                                                                                         <  7:92008/149‑1 (MQ=255)
                                          gACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTc                                                                                                                              >  8:65144/1‑126 (MQ=255)
                                           aCGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgctgct                                                                                                       >  5:26913/1‑148 (MQ=255)
                                            cGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGaa                                                                                                               >  6:568/1‑139 (MQ=255)
                                            cGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGaa                                                                                                               <  5:568/139‑1 (MQ=255)
                                                  aaGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGGTGGACTCCGGCGCGGTGTATCGCTTGCTGGCATTCAATGCGAGCTCGCACGGCGTCGACCGGGCCAACGAAGAGCTGCTGACGAAGc                                                                                               >  8:341618/1‑149 (MQ=255)
                                                   aGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGCCCAAGCt                                                                                              >  4:2319/1‑149 (MQ=255)
                                                           gggCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGccc                                                                                      <  8:241574/149‑1 (MQ=255)
                                                                 gggCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTgg                                                                                <  6:299163/149‑1 (MQ=255)
                                                                       ggTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTc                                                                                                                              >  1:172440/1‑97 (MQ=255)
                                                                       ggTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTc                                                                                                                              <  2:172440/97‑1 (MQ=255)
                                                                       ggTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGc                                                                          <  7:344857/149‑1 (MQ=255)
                                                                              gggCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCAt                                                                   >  4:166864/1‑149 (MQ=255)
                                                                                    aTTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTg                                                              <  2:199252/148‑1 (MQ=255)
                                                                                    aTTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTg                                                              <  6:26913/148‑1 (MQ=255)
                                                                                      tGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGcc                                                            <  3:2319/148‑1 (MQ=255)
                                                                                              aGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGt                                                   <  5:237423/149‑1 (MQ=255)
                                                                                                 tgggctggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAg                                                >  8:144169/1‑149 (MQ=255)
                                                                                                   ggctggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGc                                               >  2:268388/1‑148 (MQ=255)
                                                                                                     ctggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTac                                            <  5:333625/149‑1 (MQ=255)
                                                                                                        gAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTacaac                                         <  3:166864/149‑1 (MQ=255)
                                                                                                         aaGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGcc                                                            >  4:285659/1‑129 (MQ=255)
                                                                                                         aaGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGcc                                                            <  3:285659/129‑1 (MQ=255)
                                                                                                            ctgctgGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGat                                     <  1:268388/149‑1 (MQ=255)
                                                                                                               cgggACTCCGGCGCGCTGTATCGCTTGCTTGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGatcat                                  <  8:336458/147‑1 (MQ=255)
                                                                                                                  gACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATcc                                <  2:227927/148‑1 (MQ=255)
                                                                                                                     tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTacaaca                                        >  5:151970/1‑137 (MQ=255)
                                                                                                                     tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTacaaca                                        <  6:151970/137‑1 (MQ=255)
                                                                                                                     tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCt                               <  6:43774/146‑1 (MQ=255)
                                                                                                                     tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCt                               >  5:43774/1‑146 (MQ=255)
                                                                                                                       cGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAgg                          >  4:8201/1‑149 (MQ=255)
                                                                                                                               tGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGAc                  >  3:35121/1‑149 (MQ=255)
                                                                                                                                   tCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGt                >  6:264147/1‑147 (MQ=255)
                                                                                                                                    cGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAg             >  4:59108/1‑149 (MQ=255)
                                                                                                                                             gCATTCAATGCCAGCAACCGCGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGa                            >  7:129710/1‑125 (MQ=255)
                                                                                                                                             gCATTCAATGCCAGCAACCGCGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGa                            <  8:129710/125‑1 (MQ=255)
                                                                                                                                             gCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAt                                                                              >  1:63174/1‑75 (MQ=255)
                                                                                                                                             gCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAt                                                                              <  2:63174/75‑1 (MQ=255)
                                                                                                                                             gCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCAt    >  4:92103/1‑149 (MQ=255)
                                                                                                                                              cATTCAATGCCAGCTACCAAGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTCCATCGCTGCCGAGCCCTGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATc   <  3:52479/149‑1 (MQ=255)
                                                                                                                                               aTTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGAc                  <  3:334949/133‑1 (MQ=255)
                                                                                                                                               aTTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGAc                  >  4:334949/1‑133 (MQ=255)
                                                                                                                                               aTTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATcc  >  5:79570/1‑149 (MQ=255)
                                                                                                                                               aTTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATcc  >  7:43394/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGAGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCC  >  NC_002947/1978562‑1978853

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