Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,872,149 C→T 8.7% Y644Y (TAC→TAT xdhB → xanthine dehydrogenase subunit XdhB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,872,1490CT 8.7% 909.1 / 54.7 309Y644Y (TAC→TATxdhBxanthine dehydrogenase subunit XdhB
Reads supporting (aligned to +/- strand):  ref base C (126/156);  new base T (10/17);  total (136/173)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 321 total aligned reads displayed.

TTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAC  >  NC_002947/4872000‑4872298
                                                                                                                                                      |                                                                                                                                                     
ttCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATAc                                                                                                                                                       <  1:1241981/150‑1 (MQ=255)
ttCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATAc                                                                                                                                                       >  2:3795455/1‑150 (MQ=255)
 tCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACa                                                                                                                                                      <  2:2056852/150‑1 (MQ=255)
  cGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACaa                                                                                                                                                     <  2:4510967/150‑1 (MQ=255)
          tgtCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTgc                                                                                                                                             <  2:532301/150‑1 (MQ=255)
             cGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAATTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTgcgcg                                                                                                                                          <  1:1946969/150‑1 (MQ=255)
              gCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTgcgcgc                                                                                                                                         >  2:2236826/1‑150 (MQ=255)
              gCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTgcgcgc                                                                                                                                         >  1:5012411/1‑150 (MQ=255)
               cTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGcc                                                                                                                                        <  1:2649225/150‑1 (MQ=255)
                tGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCg                                                                                                                                       >  2:3041539/1‑150 (MQ=255)
                    cAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACAt                                                                                                                                   <  1:9124099/150‑1 (MQ=255)
                    cAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACAt                                                                                                                                   >  1:3422141/1‑150 (MQ=255)
                         ccGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGc                                                                                                                              >  2:678385/1‑150 (MQ=255)
                             gTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCacga                                                                                                                          >  2:5783444/1‑150 (MQ=255)
                                 tACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTg                                                                                                                      <  1:1418681/150‑1 (MQ=255)
                                   ccGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTggg                                                                                                                    >  1:537264/1‑150 (MQ=255)
                                     gCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCg                                                                                                                  >  2:8347698/1‑150 (MQ=255)
                                          ccGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCg                                                                                                             <  1:9010109/150‑1 (MQ=255)
                                            gAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCt                                                                                                           >  2:4339687/1‑150 (MQ=255)
                                            gAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCt                                                                                                           >  1:5821820/1‑150 (MQ=255)
                                                 tcttctACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAAcc                                                                                                      >  2:1777984/1‑150 (MQ=255)
                                                   ttctACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCt                                                                                                    >  2:4763394/1‑150 (MQ=255)
                                                            cGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGAc                                                                                           <  1:3152599/150‑1 (MQ=255)
                                                            cGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGAc                                                                                           <  1:5144699/150‑1 (MQ=255)
                                                              cAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGACGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACAt                                                                                         >  2:9049233/1‑150 (MQ=255)
                                                                 ccAGGCCCGCGGCCGGCCGTTCTACTCCTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCgg                                                                                      >  1:6975209/1‑150 (MQ=255)
                                                                     gCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGTGATTCGCTGAACCCTGCGATCGACATCGGCCAg                                                                                  >  2:2995121/1‑150 (MQ=255)
                                                                       ccGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGt                                                                                >  2:3171517/1‑150 (MQ=255)
                                                                       ccGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGt                                                                                <  1:5957378/150‑1 (MQ=255)
                                                                           ggccggccGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAg                                                                            <  1:4304933/150‑1 (MQ=255)
                                                                            gccggccgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAgg                                                                           <  2:3604322/150‑1 (MQ=255)
                                                                                  cgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGtt                                                                     <  2:3219360/150‑1 (MQ=255)
                                                                                  cgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGtt                                                                     <  1:7986712/150‑1 (MQ=255)
                                                                                    ttCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTc                                                                   <  1:2905837/150‑1 (MQ=255)
                                                                                       tactacTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATc                                                                >  1:4745954/1‑150 (MQ=255)
                                                                                       tactacTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATc                                                                >  1:5615977/1‑150 (MQ=255)
                                                                                          tattatgCCTATGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAg                                                             <  1:1079026/144‑1 (MQ=255)
                                                                                          tactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAg                                                             <  1:7332230/150‑1 (MQ=255)
                                                                                          tactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAg                                                             >  2:1590624/1‑150 (MQ=255)
                                                                                             ttCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGc                                                          >  1:788570/1‑150 (MQ=255)
                                                                                                gCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACAGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATg                                                       >  2:7181788/1‑150 (MQ=255)
                                                                                                  ctatGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATggg                                                     >  2:8578027/4‑150 (MQ=255)
                                                                                                  ctatGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATggg                                                     >  1:6259150/4‑150 (MQ=255)
                                                                                                   tttGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggc                                                    >  2:3725967/1‑150 (MQ=255)
                                                                                                    ttGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggct                                                   <  2:6196948/150‑1 (MQ=255)
                                                                                                      ggCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctgg                                                 <  1:6892900/150‑1 (MQ=255)
                                                                                                      ggCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGggctgg                                                 <  2:6770610/150‑1 (MQ=255)
                                                                                                       gcgcTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctggc                                                <  2:4912886/150‑1 (MQ=255)
                                                                                                         gcTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGgctggctg                                              >  2:7264681/1‑150 (MQ=255)
                                                                                                          ctgcctgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGa                                             <  1:4093617/150‑1 (MQ=255)
                                                                                                           tgcctgcGTGGAAGTGATCGTCGATACGCTGACCGTCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGac                                            <  1:4960629/150‑1 (MQ=255)
                                                                                                           tgcctgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGac                                            <  1:7573469/150‑1 (MQ=255)
                                                                                                             cctgcTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGacca                                          >  1:6636680/1‑150 (MQ=255)
                                                                                                              ctgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccac                                         >  2:2502554/1‑150 (MQ=255)
                                                                                                              ctgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccac                                         <  1:2690368/150‑1 (MQ=255)
                                                                                                               tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                        >  1:8972235/1‑150 (MQ=255)
                                                                                                               tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                        <  2:4097920/150‑1 (MQ=255)
                                                                                                               tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                        >  2:3455477/1‑150 (MQ=255)
                                                                                                               tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                        >  1:3275570/1‑150 (MQ=255)
                                                                                                               tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                        >  1:2985206/1‑150 (MQ=255)
                                                                                                               tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                        >  2:3387591/1‑150 (MQ=255)
                                                                                                                gcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCg                                       <  2:7562559/150‑1 (MQ=255)
                                                                                                                  gTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGaa                                     >  1:2062911/1‑150 (MQ=255)
                                                                                                                  gTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGaa                                     >  1:2174833/1‑150 (MQ=255)
                                                                                                                  gTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGaa                                     <  1:4569948/150‑1 (MQ=255)
                                                                                                                      aaGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGc                                 <  2:6542291/150‑1 (MQ=255)
                                                                                                                       aGTGATCGTCGATACGCTGACCGGCGAATACAAGATGGTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCt                                <  1:3795455/150‑1 (MQ=255)
                                                                                                                         gtatcggtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCtt                               <  1:7092374/144‑2 (MQ=255)
                                                                                                                          tatcggtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCttg                              <  2:466748/145‑3 (MQ=255)
                                                                                                                          gATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgt                             <  2:6391534/150‑1 (MQ=255)
                                                                                                                          gATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgt                             >  2:1137413/1‑150 (MQ=255)
                                                                                                                           atcggtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTgt                             <  2:5170993/146‑1 (MQ=255)
                                                                                                                           aTCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtg                            >  1:232075/1‑150 (MQ=255)
                                                                                                                            tcgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGa                                             >  1:1183611/1‑132 (MQ=255)
                                                                                                                            tcgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtgt                           >  2:695668/1‑150 (MQ=255)
                                                                                                                            tcgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtgt                           <  1:100998/150‑1 (MQ=255)
                                                                                                                            tcgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtgt                           <  1:8888110/150‑1 (MQ=255)
                                                                                                                            tcgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtgt                           <  2:2338632/150‑1 (MQ=255)
                                                                                                                            tcgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGAGTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTTGCTGACCACCGAAGAGCTGgtgt                           <  2:8089834/150‑1 (MQ=255)
                                                                                                                               gtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTgg                         >  1:6407132/1‑150 (MQ=255)
                                                                                                                                tcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGa                        <  2:7102245/150‑1 (MQ=255)
                                                                                                                                 cgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTggga                       <  1:8121797/150‑3 (MQ=255)
                                                                                                                                 cgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTggga                       <  1:1806969/150‑3 (MQ=255)
                                                                                                                                  caTACGCTGACCCGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCAGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAAc                      <  2:6081012/149‑1 (MQ=255)
                                                                                                                                    tACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGc                    <  2:1532399/150‑1 (MQ=255)
                                                                                                                                    tACGCTGACCGGCGAATACAAGATGGTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGc                    <  1:4845865/150‑1 (MQ=255)
                                                                                                                                       gCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaa                 >  2:1273190/1‑150 (MQ=255)
                                                                                                                                       gCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaa                 <  2:5601442/150‑1 (MQ=255)
                                                                                                                                       gCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaa                 <  2:4028045/150‑1 (MQ=255)
                                                                                                                                        cTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaaa                <  1:2363773/150‑1 (MQ=255)
                                                                                                                                        cTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaaa                <  2:4058501/150‑1 (MQ=255)
                                                                                                                                         tGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAg               <  1:8636188/150‑1 (MQ=255)
                                                                                                                                              ggCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAg          >  2:8644130/1‑150 (MQ=255)
                                                                                                                                                 gAATACAAGATGTTGCGCGCCGACAGCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtg       <  2:8866830/150‑1 (MQ=255)
                                                                                                                                                  aaTACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtga      <  2:3275570/150‑1 (MQ=255)
                                                                                                                                                   aTACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtgat     >  2:2678656/1‑150 (MQ=255)
                                                                                                                                                    tACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtgatg    <  1:1562975/150‑1 (MQ=255)
                                                                                                                                                      cAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAc  >  2:6250641/1‑150 (MQ=255)
                                                                                                                                                      cAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAc  >  1:2460374/1‑150 (MQ=255)
                                                                                                                                                      cAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAc  >  2:796633/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAC  >  NC_002947/4872000‑4872298

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: