Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,315,814 C→T 12.6% G90E (GGA→GAA)  phnN ← ribose 1,5‑bisphosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,315,8140CT12.6% 182.2 / 16.1 66G90E (GGA→GAA) phnNribose 1,5‑bisphosphokinase
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base T (4/4);  total (34/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 = 4.29e-01

TGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGC  >  NC_000913/4315671‑4315959
                                                                                                                                               |                                                                                                                                                 
tGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCg                                                                                                                                              >  2:1077146/1‑149 (MQ=255)
             gCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGAt                                                                                                                                 >  2:1212718/1‑149 (MQ=255)
               cACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCa                                                                                                                               >  2:277280/1‑149 (MQ=255)
                         cAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGc                                                                                                                                           >  2:578919/1‑127 (MQ=255)
                         cAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGc                                                                                                                                           <  1:578919/127‑1 (MQ=255)
                          aGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgcc                                                                                                                    >  1:980424/1‑149 (MQ=255)
                           ggCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgccg                                                                                                                   >  3:673089/1‑149 (MQ=255)
                            gcgcTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAAtc                                                                                                                            <  2:540918/139‑1 (MQ=255)
                            gcgcTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAAtc                                                                                                                            >  1:540918/1‑139 (MQ=255)
                             cgcTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgccgac                                                                                                                 <  1:1077146/149‑1 (MQ=255)
                                    ggaggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGc                                                                                                                                               >  3:1248218/1‑112 (MQ=255)
                                    ggaggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCataa                                                                                                          >  4:467790/1‑149 (MQ=255)
                                    ggaggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGAAAGCGCGCCCGCGCCTGCGGCAGATGGGCGAGTGAGCCGTTGACCAGCACGTCGAATCCGGc                                                                                                                                               <  4:1248218/112‑1 (MQ=255)
                                         aTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCgg                                                                                                                                                <  1:245608/106‑1 (MQ=255)
                                         aTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCgg                                                                                                                                                >  2:245608/1‑106 (MQ=255)
                                                gCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTa                                                                                              <  4:673089/149‑1 (MQ=255)
                                                  gAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGc                                                                                            <  1:1212718/149‑1 (MQ=255)
                                                        tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCa                                                                                      >  3:1054567/1‑149 (MQ=255)
                                                        tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCa                                                                                      >  4:1418304/1‑149 (MQ=255)
                                                        tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGc                                                                                                                                           >  4:413351/1‑96 (MQ=255)
                                                        tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGc                                                                                                                                           <  3:413351/96‑1 (MQ=255)
                                                         gTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCg                                                                                           >  4:654385/1‑143 (MQ=255)
                                                         gTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCg                                                                                           <  3:654385/143‑1 (MQ=255)
                                                           aaaCAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacg                                                                                                                      >  3:409601/1‑114 (MQ=255)
                                                           aaaCAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacg                                                                                                                      <  4:409601/114‑1 (MQ=255)
                                                               aGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAg                                                                               <  3:1418304/149‑1 (MQ=255)
                                                                         cagcGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCgacgccgacgcc                                                                                                              <  1:1072429/108‑1 (MQ=255)
                                                                         cagcGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCgacgccgacgcc                                                                                                              >  2:1072429/1‑108 (MQ=255)
                                                                               cGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGcc                                                               <  1:277280/149‑1 (MQ=255)
                                                                                 aTTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGcccc                                                             >  3:216823/1‑149 (MQ=255)
                                                                                     aTAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCa                                                                                                             >  1:1178394/1‑97 (MQ=255)
                                                                                     aTAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCa                                                                                                             <  2:1178394/97‑1 (MQ=255)
                                                                                          gcgcccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGAttt                                                                   >  3:890019/1‑134 (MQ=255)
                                                                                          gcgcccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGAttt                                                                   <  4:890019/134‑1 (MQ=255)
                                                                                              ccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaaa                                                >  4:511406/1‑149 (MQ=255)
                                                                                                 cgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACt                                             <  3:511406/149‑1 (MQ=255)
                                                                                                         gCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgctcg                                     >  1:1209593/1‑149 (MQ=255)
                                                                                                           aGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTc                                            >  4:172956/1‑140 (MQ=255)
                                                                                                           aGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTc                                            <  3:172956/140‑1 (MQ=255)
                                                                                                            gATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgctcgctc                                  >  3:1415503/1‑149 (MQ=255)
                                                                                                                 gcgcGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGc                             <  3:467790/149‑1 (MQ=255)
                                                                                                                     gTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGc                                                                                    >  1:564000/1‑90 (MQ=255)
                                                                                                                     gTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGc                                                                                    <  2:564000/90‑1 (MQ=255)
                                                                                                                        aGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATgg                      >  4:1314646/1‑149 (MQ=255)
                                                                                                                         gCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGt                     <  4:216823/149‑1 (MQ=255)
                                                                                                                            gTTGACCGGCACGTCGACTGCGGCGTGCAGCCAGACACCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAACGGTAACAGATTTTGCCCCGCGCGGGAAAAAAACTCCTGCTCGCTCAGGGCGATATGGttgt                  <  4:1054567/149‑3 (MQ=255)
                                                                                                                                  cAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCg                                                                                           >  3:652868/1‑70 (MQ=255)
                                                                                                                                  cAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCg                                                                                           <  4:652868/70‑1 (MQ=255)
                                                                                                                                   tgCACGCCGGATTCGACGTGCTGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgc                                        >  2:872597/2‑120 (MQ=255)
                                                                                                                                   tgCACGCCGGATTCGACGTGCTGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgc                                        <  3:1418454/119‑1 (MQ=255)
                                                                                                                                   tgCACGCCGGATTCGACGTGCTGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgc                                        <  1:872597/119‑1 (MQ=255)
                                                                                                                                   tgCACGCCGGATTCGACGTGCTGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgc                                        >  4:1418454/2‑120 (MQ=255)
                                                                                                                                   aGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTcc           <  3:1314646/149‑1 (MQ=255)
                                                                                                                                     cACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGc                                                                <  3:1147134/94‑1 (MQ=255)
                                                                                                                                     cACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGc                                                                >  4:1147134/1‑94 (MQ=255)
                                                                                                                                        gTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaaa                                                >  1:814069/1‑107 (MQ=255)
                                                                                                                                        gTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaaa                                                <  2:814069/107‑1 (MQ=255)
                                                                                                                                         tCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaa                                                 <  2:1325618/105‑1 (MQ=255)
                                                                                                                                         tCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaa                                                 >  1:1325618/1‑105 (MQ=255)
                                                                                                                                         tCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCt     <  4:1415503/149‑1 (MQ=255)
                                                                                                                                          cGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTg    >  3:483407/1‑149 (MQ=255)
                                                                                                                                           gAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAgg                               <  2:1091699/121‑1 (MQ=255)
                                                                                                                                           gAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAgg                               >  1:1091699/1‑121 (MQ=255)
                                                                                                                                            aaTCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCa                >  1:127163/1‑135 (MQ=255)
                                                                                                                                            aaTCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCa                <  2:127163/135‑1 (MQ=255)
                                                                                                                                            aaTCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGc  <  2:1209593/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
TGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGC  >  NC_000913/4315671‑4315959

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑