Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,807,254 G→A 6.9% S213F (TCC→TTC)  PP_4234 ← bifunctional aldehyde oxidase/xanthine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,807,2540GA 6.9% 213.2 / 3.9 72S213F (TCC→TTC) PP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (37/30);  new base A (3/2);  total (40/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 = 3.67e-01

GCGGGTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGCCTGGCGGAAG  >  NC_002947/4807114‑4807400
                                                                                                                                             |                                                                                                                                                  
gcgggTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCa                                                                                                                                            <  6:73627/149‑1 (MQ=255)
  gggTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTCCAGCCTGAAGAACTccat                                                                                                                                          <  2:6346/149‑4 (MQ=255)
   ggTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCAtcg                                                                                                                                         >  5:148115/1‑149 (MQ=255)
    gTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCAtcgt                                                                                                                                        >  2:14910/1‑149 (MQ=255)
     ttGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCAtcgtc                                                                                                                                       <  3:174584/149‑1 (MQ=255)
        tCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCAtcgtcgtc                                                                                                                                    <  2:226525/149‑1 (MQ=255)
          cAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCAtcgtcgtcg                                                                                                                                   >  5:99916/1‑148 (MQ=255)
             ctgctCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGGGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGGCGAACg                                                                                                                               <  8:178417/149‑1 (MQ=255)
                ctCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGc                                                                                                                             <  7:187379/148‑1 (MQ=255)
                ctCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCg                                                                                                                            >  6:70990/1‑149 (MQ=255)
                       aTGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGcc                                                                                                                     <  1:63288/149‑1 (MQ=255)
                          cgcTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGTCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGccacc                                                                                                                  <  2:257451/149‑1 (MQ=255)
                                  gTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCa                                                                                                          >  8:145215/1‑149 (MQ=255)
                                  gTAGAACAACCCCTGCTCGGTCAGGTAGCCGTGCTGGTTGGGCTGCTTCGACCACACCGCCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCa                                                                                                          >  7:189856/1‑149 (MQ=255)
                                   tAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCAc                                                                                                         >  8:148731/1‑149 (MQ=255)
                                      aacaacCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACg                                                                                                        >  2:208122/1‑147 (MQ=255)
                                        caacCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGc                                                                                                    >  7:8891/1‑149 (MQ=255)
                                         aacCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGgc                                                                                                                        >  8:73900/1‑128 (MQ=255)
                                         aacCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGgc                                                                                                                        <  7:73900/128‑1 (MQ=255)
                                             acTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGGTGAGAATGGTGGGGACCAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCCCCCGCACGAAGCGGacg                                                                                               >  8:176798/2‑147 (MQ=255)
                                               tGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGccac                                                                                                                   >  3:108705/1‑127 (MQ=255)
                                               tGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGccac                                                                                                                   <  4:108705/127‑1 (MQ=255)
                                                 cTCGGGCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTg                                                                                            <  8:346290/148‑1 (MQ=255)
                                                  tCGGTCATCTCCCCGTCCTGGATCATCTGCTTCGACCACACCGGCTTGCCGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGc                                                                                           >  4:139052/1‑148 (MQ=255)
                                                    ggTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGATCAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCa                                                                                         >  5:217891/1‑148 (MQ=255)
                                                        aGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGGGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAg                                                                                     <  6:66499/148‑1 (MQ=255)
                                                          gTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGt                                                                                  >  7:60292/1‑149 (MQ=255)
                                                           tCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTc                                                                                 <  6:99916/149‑1 (MQ=255)
                                                             gCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAggg                                                                                   <  8:110413/145‑1 (MQ=255)
                                                             gCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAggg                                                                                   >  7:110413/1‑145 (MQ=255)
                                                                 tGCTGGGTGGGCTGCCTCCACCACACCGGCTTGCGGTTGAGAATGGTGGGGATCAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTg                                                                           <  6:217891/149‑1 (MQ=255)
                                                                     ggTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGtt                                                                        >  5:369507/1‑147 (MQ=255)
                                                                       ttGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAg                                                                     <  7:148731/149‑1 (MQ=255)
                                                                         gggCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAgcg                                                                   >  4:266822/1‑149 (MQ=255)
                                                                           gctgctTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAgcggc                                                                 >  6:14089/1‑149 (MQ=255)
                                                                               ctTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCt                                                             >  5:67348/1‑149 (MQ=255)
                                                                                 tCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTGCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGATGCGGATGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTgg                                                           >  4:19776/1‑149 (MQ=255)
                                                                                  cGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTgg                                                           <  8:8891/148‑1 (MQ=255)
                                                                                  cGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTg                                                            <  1:14910/147‑1 (MQ=255)
                                                                                  cGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTg                                                            <  5:14089/147‑1 (MQ=255)
                                                                                   gACCACACCGGCTTGCGGTTGAGAATGGTGGCGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCAACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTggg                                                          >  7:349089/1‑148 (MQ=255)
                                                                                       acacCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCg                                                     >  3:98536/1‑149 (MQ=255)
                                                                                          ccGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTaa                                                  <  6:148115/149‑1 (MQ=255)
                                                                                                    gTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCct                                        <  5:70990/149‑1 (MQ=255)
                                                                                                    gTTCAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACCGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGACATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCct                                        >  2:69341/1‑149 (MQ=255)
                                                                                                     ttGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACACACCCGCACGAAGCGGAAGCTGCCATATGGGTCGCGCTGGTACAGCGGCTCCTGGGCACCGTAACGGATCGCctg                                       >  6:131517/1‑149 (MQ=255)
                                                                                                        agaATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCctgct                                     <  8:60292/148‑1 (MQ=255)
                                                                                                        agaATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGAAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCctgct                                     >  1:252629/1‑148 (MQ=255)
                                                                                                         gaATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCctgctg                                    <  3:139052/148‑1 (MQ=255)
                                                                                                         gaATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGa                                   >  6:317350/1‑149 (MQ=255)
                                                                                                           aTGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGATAACTCATGGTGGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGGGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGAtt                                 <  3:19776/149‑1 (MQ=255)
                                                                                                               tggGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAAc                                                 >  3:291168/1‑129 (MQ=255)
                                                                                                               tggGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAAc                                                 <  4:291168/129‑1 (MQ=255)
                                                                                                                  ggAACAGCATGGAGTTCTTCAGGCTGGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTg                           >  1:210024/1‑149 (MQ=255)
                                                                                                                  ggAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATc                                            >  5:323302/1‑131 (MQ=255)
                                                                                                                  ggAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATc                                            <  6:323302/131‑1 (MQ=255)
                                                                                                                   gAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCa                         <  8:163827/149‑1 (MQ=255)
                                                                                                                      cAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCt                      >  5:200390/1‑149 (MQ=255)
                                                                                                                      cAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCt                      >  4:282806/1‑149 (MQ=255)
                                                                                                                         cgatgAGTTCTCCAGCCTGAAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTgc                   >  5:287789/5‑149 (MQ=255)
                                                                                                                           tGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGGGGCGCCACCCACCCGCACGAAGCGGAAGCTGCAATAAGGGTCGCGCTGGTACATCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTgcgc                 <  1:290563/149‑1 (MQ=255)
                                                                                                                            tgAGTTCTCCAGCCTGAAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCt                >  1:348812/2‑149 (MQ=255)
                                                                                                                               gTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTgg             >  8:338564/1‑149 (MQ=255)
                                                                                                                                  cttcAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGcc           <  6:200390/148‑1 (MQ=255)
                                                                                                                                  cttcAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGCCt          >  4:168684/1‑149 (MQ=255)
                                                                                                                                   ttcAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGCCTg         >  8:289302/1‑149 (MQ=255)
                                                                                                                                      aGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGCCTgtc       <  8:189856/148‑3 (MQ=255)
                                                                                                                                       ggCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCa                         >  7:165144/1‑129 (MQ=255)
                                                                                                                                       ggCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCa                         <  8:165144/129‑1 (MQ=255)
                                                                                                                                       ggCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTgcgc                 >  4:327318/1‑137 (MQ=255)
                                                                                                                                         cTGAAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGCCTGGCGGaa   <  6:287789/149‑1 (MQ=255)
                                                                                                                                          tGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGCCTGGCGGAAg  >  6:302153/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GCGGGTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTGCGGTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGCCTGGCGGAAG  >  NC_002947/4807114‑4807400

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