Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,872,260 C→G 7.5% T681T (ACC→ACG xdhB → xanthine dehydrogenase subunit XdhB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,872,2600CG 7.5% 598.3 / 27.7 183T681T (ACC→ACGxdhBxanthine dehydrogenase subunit XdhB
Reads supporting (aligned to +/- strand):  ref base C (76/93);  new base G (7/7);  total (83/100)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

TGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCACTCCAA  >  NC_002947/4872111‑4872409
                                                                                                                                                     |                                                                                                                                                     
tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                                                                                                                                       >  2:5400787/1‑150 (MQ=255)
 gcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCg                                                                                                                                                      >  1:3347576/1‑150 (MQ=255)
   gTGGAAGTTATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGaa                                                                                                                                                    >  1:3506670/1‑150 (MQ=255)
   gTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGaa                                                                                                                                                    <  1:647358/150‑1 (MQ=255)
   gTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCACGGGCTGGCTGACCACCGaa                                                                                                                                                    >  2:5178857/1‑150 (MQ=255)
    tGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACCTCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAag                                                                                                                                                   >  2:2445551/1‑150 (MQ=255)
      gAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAagag                                                                                                                                                 >  2:1878672/1‑150 (MQ=255)
       aaGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGc                                                                                                                                                <  2:1843717/150‑1 (MQ=255)
        aGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGTAGAGCt                                                                                                                                               <  2:3209710/150‑1 (MQ=255)
         gTGCTCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTg                                                                                                                                              <  2:4320207/150‑1 (MQ=255)
          tGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTgg                                                                                                                                             <  2:4428262/150‑1 (MQ=255)
           gATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACGTCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgt                                                                                                                                            <  2:375454/150‑1 (MQ=255)
           gATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACGTCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgt                                                                                                                                            <  1:4893383/150‑1 (MQ=255)
            tcggtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTgtc                                                                                                                                           >  2:4281445/4‑149 (MQ=255)
              cgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGtgtg                                                                                                                                         >  2:3221827/1‑150 (MQ=255)
               gtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTgg                                                                                                                                        >  1:4043178/1‑150 (MQ=255)
               gtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTgg                                                                                                                                        >  1:1539893/1‑150 (MQ=255)
               gtcgATACGCGGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTgg                                                                                                                                        <  1:485293/150‑1 (MQ=255)
                tcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTggg                                                                                                                                       <  2:1540181/150‑2 (MQ=255)
                tcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCCGGCTGACCACCGAAGAGCTGGTGTGGa                                                                                                                                       >  1:2437808/1‑150 (MQ=255)
                  gATACGCTGACCGGCGAATACACGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAAc                                                                                                                                     <  2:3577580/150‑1 (MQ=255)
                       gCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaa                                                                                                                                >  2:1535666/1‑150 (MQ=255)
                        cTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGc                                                                                                                                                >  2:5263806/1‑133 (MQ=255)
                        cTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaaa                                                                                                                               <  2:585470/150‑1 (MQ=255)
                         tGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTGGGATGACAAAg                                                                                                                              <  1:2221360/150‑1 (MQ=255)
                            ccGGCGAATACACGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCa                                                                                                                           <  2:3958072/150‑1 (MQ=255)
                             cGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCaa                                                                                                                          <  1:1763383/150‑1 (MQ=255)
                             cGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCaa                                                                                                                          <  1:161093/150‑1 (MQ=255)
                              ggCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAg                                                                                                                         <  2:1754001/150‑1 (MQ=255)
                                 gAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtg                                                                                                                      >  2:1259711/1‑150 (MQ=255)
                                    tACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtgatg                                                                                                                   >  2:4907702/1‑150 (MQ=255)
                                      cAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAc                                                                                                                 <  1:210326/150‑1 (MQ=255)
                                       aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGc                                                                                                                            >  2:5352059/1‑138 (MQ=255)
                                       aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                                                                                                >  2:851437/1‑150 (MQ=255)
                                       aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                                                                                                >  2:626086/1‑150 (MQ=255)
                                       aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                                                                                                >  1:2723974/1‑150 (MQ=255)
                                         gATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCaa                                                                                                              >  1:1736417/1‑150 (MQ=255)
                                         gATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCaa                                                                                                              >  2:1697557/1‑150 (MQ=255)
                                             ttGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggc                                                                                                          >  1:4303835/1‑150 (MQ=255)
                                               gcgcgcCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggccc                                                                                                        >  1:1955983/1‑150 (MQ=255)
                                                   accGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTgtct                                                                                                                                          <  2:2942511/111‑3 (MQ=255)
                                                   accGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTgtct                                                                                                                                          >  1:2942511/2‑110 (MQ=255)
                                                         aTCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTAc                                                                                              <  1:1320114/150‑1 (MQ=255)
                                                           ccTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACaa                                                                                            <  1:10655/150‑1 (MQ=255)
                                                               cACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtgat                                                                                                                    <  2:1188045/122‑1 (MQ=255)
                                                               cACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATc                                                                                        <  1:422680/150‑1 (MQ=255)
                                                                 cgacgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATccc                                                                                      <  1:4819450/150‑1 (MQ=255)
                                                                 cgacgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATccc                                                                                      <  1:1866042/150‑1 (MQ=255)
                                                                 cgacgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATccc                                                                                      <  1:4818281/150‑1 (MQ=255)
                                                                     gTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGcc                                                                                  <  1:3638555/150‑1 (MQ=255)
                                                                     gTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCGGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGcc                                                                                  <  2:3040966/150‑1 (MQ=255)
                                                                       gggCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGt                                                                                <  1:1647274/150‑1 (MQ=255)
                                                                       gggCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGt                                                                                <  2:2860731/150‑1 (MQ=255)
                                                                       gggCGAGTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGccac                                                                                <  2:291925/150‑3 (MQ=255)
                                                                        ggCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTg                                                                               <  1:2892459/150‑1 (MQ=255)
                                                                        ggCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTg                                                                               <  1:1889746/150‑1 (MQ=255)
                                                                          cGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGc                                                                             <  1:2394037/150‑1 (MQ=255)
                                                                             ttCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGa                                                                          <  1:2857440/150‑1 (MQ=255)
                                                                               cGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACa                                                                        <  2:131094/150‑1 (MQ=255)
                                                                                gCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCTAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCACGGCCgggat                                                                       <  2:4559022/150‑5 (MQ=255)
                                                                                gCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACAt                                                                       >  1:4072088/1‑150 (MQ=255)
                                                                                gCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACAt                                                                       >  1:5337846/1‑150 (MQ=255)
                                                                                 cTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCACGGCCgggatc                                                                      <  1:4940963/150‑6 (MQ=255)
                                                                                       ccTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACaa                                                                                            >  1:4603284/1‑122 (MQ=255)
                                                                                       ccTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTg                                                                <  2:4515582/150‑1 (MQ=255)
                                                                                         tGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGa                                                              <  1:2760828/150‑1 (MQ=255)
                                                                                          gccaacGACATCCGCCAGGAGGAGGCTGGTAACATCCAGGGCAAGGGCTGCCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGAc                                                             <  2:5370309/145‑1 (MQ=255)
                                                                                          gCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGAc                                                             <  2:1351201/150‑1 (MQ=255)
                                                                                           cGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACt                                                            >  2:1597512/1‑150 (MQ=255)
                                                                                                gACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGCTCCCGGCCGTGGCCGACATGCCACTGGACTTGCgg                                                       <  1:1224469/150‑1 (MQ=255)
                                                                                                gACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCgg                                                       <  2:2163641/150‑1 (MQ=255)
                                                                                                 aCATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAag                                                                                                                                                   >  2:3512784/1‑57 (MQ=255)
                                                                                                 aCATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCggg                                                      >  1:1468479/1‑150 (MQ=255)
                                                                                                 aCATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCggg                                                      >  1:4509808/1‑150 (MQ=255)
                                                                                                   aTCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTg                                                    >  1:4912951/1‑150 (MQ=255)
                                                                                                   aTCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTg                                                    >  2:2651877/1‑150 (MQ=255)
                                                                                                   aTCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTg                                                    >  2:117666/1‑150 (MQ=255)
                                                                                                     cGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGaa                                                  >  1:1115625/1‑150 (MQ=255)
                                                                                                      ggCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAg                                                 >  2:4605673/1‑150 (MQ=255)
                                                                                                       gCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGACGc                                                >  2:419234/1‑150 (MQ=255)
                                                                                                        ccAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCt                                               >  2:4079947/1‑150 (MQ=255)
                                                                                                            gTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCtggtg                                           >  2:2312718/1‑150 (MQ=255)
                                                                                                            gTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCtggtg                                           <  1:5116762/150‑1 (MQ=255)
                                                                                                              ggAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGga                                         <  1:1441431/150‑1 (MQ=255)
                                                                                                               gAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtgatga                                                                                                                  <  2:228958/76‑1 (MQ=255)
                                                                                                                 gggTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGaa                                      <  1:3660153/150‑1 (MQ=255)
                                                                                                                        ttCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAg                               >  2:2009347/1‑150 (MQ=255)
                                                                                                                           aTCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAAc                            >  1:1879556/1‑150 (MQ=255)
                                                                                                                             ccAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGATCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAAccc                          >  2:552702/1‑150 (MQ=255)
                                                                                                                             ccAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAAccc                          <  1:4711078/150‑1 (MQ=255)
                                                                                                                              cAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCg                         <  1:5174492/150‑1 (MQ=255)
                                                                                                                              cAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCg                         <  1:3221827/150‑1 (MQ=255)
                                                                                                                                             ctgACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTc          <  2:948084/150‑1 (MQ=255)
                                                                                                                                             ctgACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTc          >  2:791486/1‑150 (MQ=255)
                                                                                                                                             ctgACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTc          <  1:5421815/150‑1 (MQ=255)
                                                                                                                                               gACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCa        <  2:3575692/150‑1 (MQ=255)
                                                                                                                                                accaccGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCAc       >  2:3678598/1‑150 (MQ=255)
                                                                                                                                                  caccGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCACTc     >  1:850166/1‑150 (MQ=255)
                                                                                                                                                   accGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCACTcc    >  1:3873625/1‑150 (MQ=255)
                                                                                                                                                    ccGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCACTCCa   >  1:263039/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCACTCCAA  >  NC_002947/4872111‑4872409

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