Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,297,638 G→A 6.5% R154C (CGC→TGC)  PP_3768 ← shikimate 5‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,297,6380GA 6.5% 180.8 / 2.1 61R154C (CGC→TGC) PP_3768shikimate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (32/25);  new base A (2/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GCCAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGCCAGGTAATCGGTGT  >  NC_002947/4297499‑4297787
                                                                                                                                           |                                                                                                                                                     
gCCAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACgatg                                                                                                                                             >  1:67464/1‑150 (MQ=255)
  cAGCTCTCCTGTCTCCATGCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACgatgat                                                                                                                                           <  2:974593/150‑1 (MQ=255)
        ttctGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCAAACACGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTc                                                                                                                                     >  1:991221/1‑150 (MQ=255)
              cTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCaaa                                                                                                                               >  2:988410/1‑150 (MQ=255)
               tCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAAc                                                                                                                              >  1:632678/1‑150 (MQ=255)
               tCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCACGCAACgagca                                                                                                                                            <  2:287580/136‑4 (MQ=255)
               tCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCACGCAACgagca                                                                                                                                            >  1:287580/1‑133 (MQ=255)
               tCCAGCCCGCCTGCCATACCGATCGGCGACAAGTTCAACAGGATGGGCGTGCAAATTGCGCCTGGCTCGGGCACTAGAAGACGTCCAGGCACACCCGCTAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGCAGTATCCTTCGGCAAAc                                                                                                                              >  1:875889/1‑150 (MQ=255)
                   gCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCgg                                                                                                                          <  2:1047248/150‑1 (MQ=255)
                    cccGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGc                                                                                                                         >  2:517325/1‑150 (MQ=255)
                    cccGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGc                                                                                                                         >  1:613011/1‑150 (MQ=255)
                           gCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCACGCAACgag                                                                                                                                              <  2:1050874/122‑2 (MQ=255)
                           gCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCACGCAACgag                                                                                                                                              >  1:1050874/1‑121 (MQ=255)
                                 ccGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGgcg                                                                                                            >  1:83013/1‑150 (MQ=255)
                                        gCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCa                                                                                                     >  2:420125/1‑150 (MQ=255)
                                           tCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACgg                                                                                                  <  2:875889/150‑1 (MQ=255)
                                           tCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACgg                                                                                                  >  1:1059971/1‑150 (MQ=255)
                                                 tCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTgg                                                                                            <  2:1059971/150‑1 (MQ=255)
                                                      cgcctGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCAAGGCCGTGGCCagc                                                                                       <  2:991221/146‑3 (MQ=255)
                                                        cATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAGTGCGCTGGCCACCGCCTTGGCCATGcc                                                                                     >  2:454629/1‑150 (MQ=255)
                                                                       gTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAg                                                                      <  2:1200835/150‑1 (MQ=255)
                                                                       gTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAg                                                                      <  2:632678/150‑1 (MQ=255)
                                                                        tCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAgg                                                                     >  1:601875/1‑150 (MQ=255)
                                                                          gCCTGGCTCAGGCACCCAGCGATAACCACAGACAACCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGc                                                                   >  1:533247/1‑150 (MQ=255)
                                                                               gCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAgg                                                              >  1:184219/1‑150 (MQ=255)
                                                                                cTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGc                                                             <  2:768303/150‑1 (MQ=255)
                                                                                         ccAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGTAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCgg                                                    >  1:543142/1‑150 (MQ=255)
                                                                                          cAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGa                                                   >  2:115544/1‑150 (MQ=255)
                                                                                           aGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATTATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTg                                                                                             <  1:667741/107‑1 (MQ=255)
                                                                                           aGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATTATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTg                                                                                             >  2:667741/1‑107 (MQ=255)
                                                                                            gCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCg                                                                                                                                    <  2:294219/67‑1 (MQ=255)
                                                                                            gCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCg                                                                                                                                    >  1:294219/1‑67 (MQ=255)
                                                                                              gATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGa                                               >  2:664279/1‑150 (MQ=255)
                                                                                               aTAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGAc                                              <  1:805084/150‑1 (MQ=255)
                                                                                                    cacaGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGAt                                         <  1:988410/150‑1 (MQ=255)
                                                                                                      caGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGc                                       >  1:262639/1‑150 (MQ=255)
                                                                                                         aCATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGa                                    <  2:67464/150‑1 (MQ=255)
                                                                                                         aCATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGa                                    <  1:642069/150‑1 (MQ=255)
                                                                                                         aCATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGa                                    <  2:613011/150‑1 (MQ=255)
                                                                                                         aCATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGa                                    <  2:510701/150‑1 (MQ=255)
                                                                                                           aTCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGc                                  >  2:710750/1‑150 (MQ=255)
                                                                                                            tCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGcc                                 <  2:601875/150‑1 (MQ=255)
                                                                                                              cGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAg                               >  1:573202/1‑150 (MQ=255)
                                                                                                               gCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGc                                       >  2:1074969/1‑141 (MQ=255)
                                                                                                               gCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGc                                       <  1:1074969/141‑1 (MQ=255)
                                                                                                                ccAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCa                             >  1:851270/1‑150 (MQ=255)
                                                                                                                    cTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTg                         >  1:1037478/1‑150 (MQ=255)
                                                                                                                     tGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTgg                        <  2:1037478/150‑1 (MQ=255)
                                                                                                                      gACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGgc                       <  2:851270/150‑1 (MQ=255)
                                                                                                                        cGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGgcgc                     >  1:315252/1‑150 (MQ=255)
                                                                                                                           ccGGTCTGCTCGTTGCGTGCAACCATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCAACGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACg                  <  1:421922/150‑1 (MQ=255)
                                                                                                                           ccGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAAt                           >  1:729360/1‑141 (MQ=255)
                                                                                                                           ccGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAAt                           <  2:729360/141‑1 (MQ=255)
                                                                                                                           ccGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACg                  >  1:996470/1‑150 (MQ=255)
                                                                                                                              gCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGcc               >  1:477477/1‑150 (MQ=255)
                                                                                                                                cTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGcc                                                       <  2:1014733/108‑1 (MQ=255)
                                                                                                                                cTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGcc                                                       >  1:1014733/1‑108 (MQ=255)
                                                                                                                                   cTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCa                                                                       >  2:394261/1‑89 (MQ=255)
                                                                                                                                   cTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCACCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCa                                                                       <  1:394261/89‑1 (MQ=255)
                                                                                                                                    tCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGCCAGGTaa         >  2:788805/1‑150 (MQ=255)
                                                                                                                                         gCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGCCAGGTAATCGgt    <  2:262639/150‑1 (MQ=255)
                                                                                                                                           gTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGCCAGGTAATCGgtgt  >  2:75807/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
GCCAGCTCTTCTGCCTCCAGCCCGCCTGCCATACCGATCGGCGTCACGTTCACCAGCATGGGCGGGCAAATGTCGCCTGGCTCAGGCACCCAGCGATAACCACAGACATCCGCCAACTGACGCCCGGCCTGCTCGTTGCGTGCAACGATGATGCCTTCGGCAAACCCGGCATCGCGCAAGGCGCTGGCCACGGCCTTGGCCATGCCACCGCTGCCTCGCAGGGCAAAGGCCGTGCCCGGATCGACCTGATGCTGAGCCAGCAATTGGCGCACGGCCAGGTAATCGGTGT  >  NC_002947/4297499‑4297787

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