Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,126,848 G→T 10.0% H193N (CAC→AAC) ‡ PP_2744 ← ribose‑phosphate pyrophosphokinase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,126,8480GT10.0% 181.6 / 9.6 60H193N (CAC→AAC) ‡PP_2744ribose‑phosphate pyrophosphokinase family protein
Reads supporting (aligned to +/- strand):  ref base G (31/23);  new base T (3/3);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

CAGCGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCG  >  NC_002947/3126705‑3126989
                                                                                                                                                |                                                                                                                                             
cagcGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGCCCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGcc                                                                                                                                           >  7:162656/1‑148 (MQ=255)
cagcGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGcc                                                                                                                                           <  8:162656/148‑1 (MQ=255)
       tGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTg                                                                                                                                   >  3:153651/1‑149 (MQ=255)
         ccccGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGc                                                                                                                                 >  4:208613/1‑149 (MQ=255)
         ccccGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGc                                                                                                                                 >  3:282777/1‑149 (MQ=255)
           ccGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCg                                                                                                                                >  1:318110/1‑148 (MQ=255)
               gCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGGCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCGAGCAGTTTCGCCAGGCCCTGGCgcacc                                                                                                                            >  5:199484/1‑145 (MQ=255)
               gCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCACCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAggg                                                                                                                                       >  5:333782/1‑137 (MQ=255)
               gCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCACCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAggg                                                                                                                                       <  6:333782/137‑1 (MQ=255)
                cTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCACCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGTTTCGCCAGGCCCTGGCgcacct                                                                                                                           <  5:3392/148‑6 (MQ=255)
                   aTCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTc                                                                                                                       >  7:321777/1‑149 (MQ=255)
                    tCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTgg                                                                                                                                  <  5:127218/137‑1 (MQ=255)
                    tCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTgg                                                                                                                                  >  6:127218/1‑137 (MQ=255)
                      aGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGGCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCGAGCAGTTTCGCCAGGCCCTGGCGCACCTGCtggc                                                                                                                     <  6:199484/149‑4 (MQ=38)
                              gAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAAc                                                                                                                            <  7:74079/133‑1 (MQ=255)
                              gAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAAc                                                                                                                            >  8:74079/1‑133 (MQ=255)
                                  aCGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGt                                                                                                                <  5:272336/141‑1 (MQ=255)
                                  aCGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGt                                                                                                                >  6:272336/1‑141 (MQ=255)
                                       gACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCAc                                                                                                      >  7:159340/1‑146 (MQ=255)
                                       gACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCAc                                                                                                      <  8:159340/146‑1 (MQ=255)
                                       gACGGTGCTGCCTGCCACGTGACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCTGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGcgcg                                                                                                                 <  6:86111/135‑1 (MQ=255)
                                       gACGGTGCTGCCTGCCACGTAACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGcgcg                                                                                                                 >  5:86111/1‑135 (MQ=255)
                                              ctgcctgcCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCgg                                                                                            <  1:313767/149‑1 (MQ=255)
                                               tgcctgccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGt                                                                                           >  7:94878/1‑149 (MQ=255)
                                                 cctgccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGt                                                                                           <  8:94878/147‑1 (MQ=255)
                                                     ccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACg                                                                                     >  3:301841/1‑149 (MQ=255)
                                                           tACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGt                                                                                >  1:281907/1‑148 (MQ=255)
                                                                 ccAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGcc                                                                         <  4:153651/149‑1 (MQ=255)
                                                                  cAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCt                                                                        <  2:318110/149‑1 (MQ=255)
                                                                   aGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTg                                                                       >  3:308116/1‑149 (MQ=255)
                                                                    ggTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGc                                                                      >  1:262767/1‑149 (MQ=255)
                                                                    ggTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGc                                                                      >  3:85099/1‑149 (MQ=255)
                                                                                tCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGc                                                          >  3:200074/1‑149 (MQ=255)
                                                                                    gcctgcctGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTAGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCgggg                                                      >  7:9043/1‑149 (MQ=255)
                                                                                     cctgcctgCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGc                                                     <  8:321777/149‑1 (MQ=255)
                                                                                         cctgcTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAg                                                  <  7:33895/148‑1 (MQ=255)
                                                                                          ctgctgGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGc                                                <  2:262767/149‑1 (MQ=255)
                                                                                                gCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTg                                            <  4:282777/147‑1 (MQ=255)
                                                                                                  ggTGTTTCTCCATCAAGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGc                                         >  8:297477/1‑148 (MQ=255)
                                                                                                     gtTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGaa                                      >  6:181839/1‑148 (MQ=255)
                                                                                                      tttCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGaa                                      >  3:326783/1‑147 (MQ=255)
                                                                                                      tttCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGaa                                      >  3:317229/1‑147 (MQ=255)
                                                                                                      tttCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCAAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGaa                                      >  4:370357/1‑147 (MQ=255)
                                                                                                            catcatGTCCACCGAGACCGGTCGCGGGAGGAGGTGCGCCAGGGCCTGGCGAAAGTGCTCGGCGCGTTTGACGCGACGGATATCCGGTTACACGACCGTGACCGCCTGCTCACCGAGTTGCTTGGCCAGGCGTTGGGCGAATAGCt                                 <  4:157823/146‑1 (MQ=255)
                                                                                                            catcatGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCt                                 <  4:85099/146‑1 (MQ=255)
                                                                                                            catcatGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTgg                               <  2:281907/148‑1 (MQ=255)
                                                                                                                atGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGCGCCAGGCGTTGGGCGAATAGCTGGGTGCa                          >  4:48786/1‑149 (MQ=255)
                                                                                                                        cGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGt                  >  3:134530/1‑149 (MQ=255)
                                                                                                                         gagaCCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGtt                 <  5:181839/149‑1 (MQ=255)
                                                                                                                              cGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGc            >  4:337271/1‑149 (MQ=255)
                                                                                                                               ggTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCt           >  4:39553/1‑149 (MQ=255)
                                                                                                                                 tCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCg         <  4:134530/149‑1 (MQ=255)
                                                                                                                                 tCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCg         >  3:125955/1‑149 (MQ=255)
                                                                                                                                 tCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCg         <  3:208613/149‑1 (MQ=255)
                                                                                                                                  cgcgCCGAGCAGTTTCGCCAGGCCCTGGCGCACCTGCTGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCa                                                   >  2:74978/1‑106 (MQ=21)
                                                                                                                                  cgcgCCGAGCAGTTTCGCCAGGCCCTGGCGCACCTGCTGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCa                                                   <  1:74978/106‑1 (MQ=21)
                                                                                                                                   gcgccgagcagTTTCGCCAGGCCCTGGCGCACCTGCTGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCgg        >  2:107937/7‑148 (MQ=37)
                                                                                                                                     gcCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTc                >  5:253314/1‑138 (MQ=255)
                                                                                                                                     gcCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTc                <  6:253314/138‑1 (MQ=255)
                                                                                                                                         agcagGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCg  <  4:200074/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CAGCGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCG  >  NC_002947/3126705‑3126989

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: