Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,872,212 C→T 10.4% I665I (ATC→ATT xdhB → xanthine dehydrogenase subunit XdhB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,872,2120CT10.4% 590.4 / 44.7 211I665I (ATC→ATTxdhBxanthine dehydrogenase subunit XdhB
Reads supporting (aligned to +/- strand):  ref base C (91/98);  new base T (9/13);  total (100/111)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

AGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAA  >  NC_002947/4872063‑4872361
                                                                                                                                                      |                                                                                                                                                     
aGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATc                                                                                                                                                       >  2:4800035/1‑150 (MQ=255)
  ccAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCgg                                                                                                                                                     <  1:368529/150‑1 (MQ=255)
  ccAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCgg                                                                                                                                                     >  1:1351201/1‑150 (MQ=255)
  ccAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCgg                                                                                                                                                     >  2:1763383/1‑150 (MQ=255)
   cAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGc                                                                                                                                                    <  2:4697594/150‑1 (MQ=255)
   cAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGc                                                                                                                                                    >  2:3757266/1‑150 (MQ=255)
        ccGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGt                                                                                                                                               <  1:1463405/150‑1 (MQ=255)
          gcgGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGACCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTgg                                                                                                                                             <  1:3778666/150‑1 (MQ=255)
           cggccggcCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGa                                                                                                                                            <  1:884320/150‑1 (MQ=255)
              ccggccgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAggg                                                                                                                                         <  2:2127114/150‑1 (MQ=255)
               cggccgTTCTACTACTACGCCTTTGGCGCAGCCTGCGTGGAATTGATCGTCGATACGCTGACCGGCGAATAAAAGATGTTGCGCGCCCACAAACTGCACGACCTGGGCGATTAGCTGAACCCAGCCATCCACATCGGCCAGGTGGAGggt                                                                                                                                        >  1:702913/1‑150 (MQ=255)
                   cgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGtt                                                                                                                                    <  2:1036529/150‑1 (MQ=255)
                   cgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGtt                                                                                                                                    <  2:3975729/150‑1 (MQ=255)
                        tactacTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATc                                                                                                                               <  1:2492418/150‑1 (MQ=255)
                         actactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGggtg                                                                                                                                       >  2:4843757/1‑141 (MQ=255)
                         actactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATcc                                                                                                                              >  2:4139996/1‑150 (MQ=255)
                          ctactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCAGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCa                                                                                                                             >  1:3054976/1‑150 (MQ=255)
                           tactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAg                                                                                                                            <  1:982668/150‑1 (MQ=255)
                            attatgCCTATGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAgg                                                                                                                           <  1:3957722/145‑1 (MQ=255)
                             ctTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAggg                                                                                                                          <  1:4895396/150‑1 (MQ=255)
                              ttCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGc                                                                                                                         >  2:422671/1‑150 (MQ=255)
                                  ccTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATgg                                                                                                                     >  2:4711078/1‑150 (MQ=255)
                                  ccTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATgg                                                                                                                     >  2:3477630/1‑150 (MQ=255)
                                  ccTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATgg                                                                                                                     >  2:325297/1‑150 (MQ=255)
                                  ccTATGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATgg                                                                                                                     <  1:1721397/150‑1 (MQ=255)
                                   cTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATggg                                                                                                                    >  2:2760828/1‑150 (MQ=255)
                                    tttGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggc                                                                                                                   <  1:2350362/150‑1 (MQ=255)
                                    tttGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggc                                                                                                                   >  1:1874692/1‑150 (MQ=255)
                                        gcgcTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctggc                                                                                                               >  1:4476319/1‑150 (MQ=255)
                                        gcgcTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctggc                                                                                                               >  2:502899/1‑150 (MQ=255)
                                           ctgcctgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGa                                                                                                            <  1:1845516/150‑1 (MQ=255)
                                            tgcctgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGac                                                                                                           <  1:323991/150‑1 (MQ=255)
                                            tgcctgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGac                                                                                                           <  1:2218900/150‑1 (MQ=255)
                                            tgcctacGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGac                                                                                                           <  1:531743/150‑1 (MQ=255)
                                             gcctgcTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAaggc                                                                                                                                        >  1:1034351/1‑119 (MQ=255)
                                                tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                                                                                       >  1:3078179/1‑150 (MQ=255)
                                                tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                                                                                       >  2:5400707/1‑150 (MQ=255)
                                                 gcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCg                                                                                                      >  1:3347576/1‑150 (MQ=255)
                                                   gTGGAAGTTATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGaa                                                                                                    >  1:3506670/1‑150 (MQ=255)
                                                   gTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCACGGGCTGGCTGACCACCGaa                                                                                                    >  2:5178857/1‑150 (MQ=255)
                                                    tGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACCTCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAag                                                                                                   >  2:2445551/1‑150 (MQ=255)
                                                       aaGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGc                                                                                                <  2:1843717/150‑1 (MQ=255)
                                                         gTGCTCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTg                                                                                              <  2:4320207/150‑1 (MQ=255)
                                                          tGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTgg                                                                                             <  2:4428262/150‑1 (MQ=255)
                                                          gtttcggtcgATACGCTGACCGGCGAATATGTGGTGGATCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCtt                                                                                              <  2:536648/144‑2 (MQ=255)
                                                            aTCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCCCCGAAGAGCTGgtg                                                                                           <  2:702913/150‑1 (MQ=255)
                                                             tcggtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTgtc                                                                                           >  2:4281445/4‑149 (MQ=255)
                                                              cgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGtgtg                                                                                         >  2:3221827/1‑150 (MQ=255)
                                                              cgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGtgtg                                                                                         >  2:1339306/1‑150 (MQ=255)
                                                              cgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGtgtg                                                                                         <  1:2932265/150‑1 (MQ=255)
                                                              cgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGtgtg                                                                                         <  1:2667110/150‑1 (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:4539750/150‑2 (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)
                                                                             ccGGCGAATACACGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCa                                                                           <  2:3958072/150‑1 (MQ=255)
                                                                             ccGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCa                                                                           >  2:1722048/1‑150 (MQ=255)
                                                                             ccGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCa                                                                           >  1:131094/1‑150 (MQ=255)
                                                                              cGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGt                                                                                                                                               <  1:2397623/81‑1 (MQ=255)
                                                                              cGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCaa                                                                          <  1:1763383/150‑1 (MQ=255)
                                                                              cGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCaa                                                                          <  1:5119700/150‑1 (MQ=255)
                                                                                gCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGc                                                                        >  2:5421815/1‑150 (MQ=255)
                                                                                  gAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtg                                                                      >  2:1259711/1‑150 (MQ=255)
                                                                                       cAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAc                                                                 <  1:210326/150‑1 (MQ=255)
                                                                                        aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGc                                                                            <  1:5352059/138‑1 (MQ=255)
                                                                                        aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                                                >  2:626086/1‑150 (MQ=255)
                                                                                        aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                                                >  2:851437/1‑150 (MQ=255)
                                                                                          gATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGAACGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCaa                                                              >  1:1736442/1‑150 (MQ=255)
                                                                                                      cGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAg                                                  >  1:1473918/1‑150 (MQ=255)
                                                                                                         cATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTa                                               >  2:138359/1‑150 (MQ=255)
                                                                                                            ccTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctgg                                                                                                                >  1:4519697/1‑82 (MQ=255)
                                                                                                            ccTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctgg                                                                                                                <  2:4519107/82‑1 (MQ=255)
                                                                                                             cTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACaat                                           >  1:1644156/1‑149 (MQ=255)
                                                                                                                cACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtgat                                                                    >  1:1188045/1‑122 (MQ=255)
                                                                                                                cACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtgat                                                                    <  2:1188045/122‑1 (MQ=255)
                                                                                                                cACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATc                                        <  1:422671/150‑1 (MQ=255)
                                                                                                                  cgacgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATccc                                      <  2:1658003/150‑1 (MQ=255)
                                                                                                                  cgacgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATccc                                      <  1:4818281/150‑1 (MQ=255)
                                                                                                                  cgacgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATccc                                      <  1:4819450/150‑1 (MQ=255)
                                                                                                                      gTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCTTTGGCTGGCTGACCCCCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGcc                                  <  2:3040891/150‑1 (MQ=255)
                                                                                                                      gTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGcc                                  <  1:3638555/150‑1 (MQ=255)
                                                                                                                        gggCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGt                                >  2:3200258/1‑150 (MQ=255)
                                                                                                                        gggCGAGTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGccac                                <  2:291925/150‑3 (MQ=255)
                                                                                                                         ggCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTg                               <  1:1889746/150‑1 (MQ=255)
                                                                                                                            gATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGcc                            <  2:2877239/150‑1 (MQ=255)
                                                                                                                              ttCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGa                          >  2:3660153/1‑150 (MQ=255)
                                                                                                                                cGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACa                        <  2:131094/150‑1 (MQ=255)
                                                                                                                                 gCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACAt                       >  1:4515582/1‑150 (MQ=255)
                                                                                                                                 gCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACAt                       >  1:4072088/1‑150 (MQ=255)
                                                                                                                                  cTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCACGGCCgggatc                      <  1:4940963/150‑6 (MQ=255)
                                                                                                                                       cccTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTCCAAGATCCCGGCCGTGGCCGACATGCCACt                 >  2:226603/1‑150 (MQ=255)
                                                                                                                                        ccTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACaa                                            >  1:4603284/1‑122 (MQ=255)
                                                                                                                                          tGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGa              <  1:2760828/150‑1 (MQ=255)
                                                                                                                                           gCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGAc             <  2:1558166/150‑1 (MQ=255)
                                                                                                                                                      cGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGaa  <  1:288364/150‑1 (MQ=255)
                                                                                                                                                      cGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGaa  >  1:1115625/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAA  >  NC_002947/4872063‑4872361

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: