Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,870,087 G→A 6.7% V598M (GTG→ATG)  xdhB → xanthine dehydrogenase subunit XdhB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,870,0870GA 6.7% 168.8 / 2.9 60V598M (GTG→ATG) xdhBxanthine dehydrogenase subunit XdhB
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base A (2/2);  total (28/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.62e-01

TCCTCAAGAAACGCCTGACCGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAA  >  NC_002947/4869944‑4870229
                                                                                                                                                |                                                                                                                                              
tCCTCAAGAAACGCCTGACCGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGAACAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGc                                                                                                                                          >  1:214235/1‑150 (MQ=255)
 ccTCAAGAAACGCCTGACCGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCt                                                                                                                                         <  2:662960/150‑1 (MQ=255)
   tCAAGAAACGCCTGACCGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGt                                                                                                                                       >  1:302910/1‑150 (MQ=255)
   tCAAGAAACGCCTGACCGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGt                                                                                                                                       >  1:148162/1‑150 (MQ=255)
   tCAAGAAACGCCTGACCGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGATGTCGCTGt                                                                                                                                       >  2:154458/1‑150 (MQ=255)
   tCAAGAAACGCCTGACCGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGATGTCGCTGt                                                                                                                                       >  2:154509/1‑150 (MQ=255)
     aaGAAACGCCTGACCGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTcc                                                                                                                                     >  2:413254/1‑150 (MQ=255)
             ccTGACCGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGc                                                                                                                                  >  2:667736/1‑145 (MQ=255)
             ccTGACCGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGc                                                                                                                                  <  1:667736/145‑1 (MQ=255)
                 aCCGAGTTTGCCGCACGCCACTAGCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGACATGTGCGGGTGCGCGACCAGATCGTCAGATTCGAATAACTGTGGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTc                                                                                                                         <  1:349782/150‑1 (MQ=255)
                   cGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTa                                                                                                                       <  2:834395/150‑1 (MQ=255)
                   cGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTa                                                                                                                       >  1:1123344/1‑150 (MQ=255)
                    gAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGCACGTCAGCTTCGATCAGCGTGGTCGTCAGGCCCACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTAc                                                                                                                      <  1:24125/150‑1 (MQ=255)
                       tttGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGc                                                                                                                   >  2:133047/1‑150 (MQ=255)
                            cGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGcc                                                                                                              >  1:721697/1‑150 (MQ=255)
                               aCGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGaa                                                                                                           <  2:3013/150‑1 (MQ=255)
                                    aCTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGAtct                                                                                                      >  1:581599/1‑150 (MQ=255)
                                        tCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGATGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTa                                                                                                  <  1:154509/150‑1 (MQ=255)
                                        tCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGATGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTa                                                                                                  <  1:154458/150‑1 (MQ=255)
                                                 cGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGt                                                                                                                           <  1:1030538/116‑1 (MQ=255)
                                                 cGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGt                                                                                                                           >  2:1030538/1‑116 (MQ=255)
                                                 cGAGGAAGACGTCGAATCTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAg                                                                                         >  2:761411/1‑150 (MQ=255)
                                                         aCGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCAc                                                                                                                 >  1:1177397/1‑118 (MQ=255)
                                                         aCGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCAc                                                                                                                 <  2:1177397/118‑1 (MQ=255)
                                                         aCGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGccc                                                                                 >  2:1093159/1‑150 (MQ=255)
                                                                       aaCGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTc                                                                   <  2:1123344/150‑1 (MQ=255)
                                                                       aaCGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTc                                                                   >  2:83343/1‑150 (MQ=255)
                                                                             cATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAAATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCtactac                                                             >  1:276863/1‑150 (MQ=255)
                                                                              aTGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCtactacg                                                            <  2:880712/150‑2 (MQ=255)
                                                                               tgtgCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCTGGTGTCGCAGTCCACTTCCGGGTTCTACCGCACGCCGAACACCTTCTACGTCCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACtt                                                           <  1:1253725/150‑1 (MQ=255)
                                                                                gtgCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTc                                                          <  2:302910/150‑1 (MQ=255)
                                                                                gtgCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTc                                                          >  2:728829/1‑150 (MQ=255)
                                                                                      gTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCttt                                                    <  2:409125/150‑1 (MQ=255)
                                                                                       tGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTg                                                   <  2:1274600/150‑1 (MQ=255)
                                                                                         cgcgACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGgc                                                 <  1:1070678/150‑1 (MQ=255)
                                                                                          gcgACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCGTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGgcg                                                >  2:604012/1‑150 (MQ=255)
                                                                                              ccAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGctgc                                            <  1:761411/150‑1 (MQ=255)
                                                                                                 gATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGAtcttc                                                                                                    <  1:670226/91‑1 (MQ=255)
                                                                                                 gATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGAtcttc                                                                                                    >  2:670226/1‑91 (MQ=255)
                                                                                                    catccgCTTCGAGCAACTGGCGCAGCGCGCCTACGTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGTCCGCGGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGt                                      <  2:276863/145‑1 (MQ=255)
                                                                                                     gTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCTCTGTCTGCGTg                                     <  2:148162/150‑1 (MQ=255)
                                                                                                       cAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGgcg                                                <  2:131758/137‑1 (MQ=255)
                                                                                                       cAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGgcg                                                >  1:131758/1‑137 (MQ=255)
                                                                                                           ttCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTg                               >  2:493801/1‑150 (MQ=255)
                                                                                                           ttCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCGTTGGCGCTGCCTGCGTGGAAGTg                               >  2:1184097/1‑150 (MQ=255)
                                                                                                              gAGCAACTGGTGCAGCAGGCCTAGTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGAtc                            <  2:196525/150‑1 (MQ=255)
                                                                                                                   aCTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCcgcg                                                                              >  2:491949/1‑95 (MQ=255)
                                                                                                                   aCTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCcgcg                                                                              <  1:491949/95‑1 (MQ=255)
                                                                                                                   aCTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGcc                                                       <  2:652610/118‑1 (MQ=255)
                                                                                                                   aCTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGcc                                                       >  1:652610/1‑118 (MQ=255)
                                                                                                                   aCTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGa                       <  2:581599/150‑1 (MQ=255)
                                                                                                                        gTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACg                   >  2:439551/1‑150 (MQ=255)
                                                                                                                            agcagGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGa               <  2:641191/150‑1 (MQ=255)
                                                                                                                            agcagGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGa               <  2:721697/150‑1 (MQ=255)
                                                                                                                            agcagGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGa               >  1:194715/1‑150 (MQ=255)
                                                                                                                            agcagGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGa               <  1:1093159/150‑1 (MQ=255)
                                                                                                                             gcagGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCtt                                                     >  2:735317/1‑111 (MQ=255)
                                                                                                                             gcagGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCtt                                                     <  1:735317/111‑1 (MQ=255)
                                                                                                                                 gCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTCGAAGTGATCGTCGATACGCTGACCGGc          <  2:194715/150‑1 (MQ=255)
                                                                                                                                         cGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACaa  <  1:439551/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TCCTCAAGAAACGCCTGACCGAGTTTGCCGCACGGCACTATCAAGTGACCGAGGAAGACGTCGAATTTCGCAACGGCCATGTGCGGGTGCGCGACCAGATCGTCAGCTTCGAGCAACTGGTGCAGCAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAA  >  NC_002947/4869944‑4870229

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