Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,872,260 C→G 5.2% 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 5.2% 942.7 / 28.5 279T681T (ACC→ACGxdhBxanthine dehydrogenase subunit XdhB
Reads supporting (aligned to +/- strand):  ref base C (134/130);  new base G (3/12);  total (137/142)
Fisher's exact test for biased strand distribution p-value = 3.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

TGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCACTCCAA  >  NC_002947/4872111‑4872409
                                                                                                                                                     |                                                                                                                                                     
tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                                                                                                                                       >  2:3387591/1‑150 (MQ=255)
     ggAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAaga                                                                                                                                                  <  2:1245194/150‑1 (MQ=255)
       aaGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGc                                                                                                                                                <  2:6542291/150‑1 (MQ=255)
        aGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCt                                                                                                                                               <  1:3795133/150‑1 (MQ=255)
         gtatcggtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCtt                                                                                                                                              <  1:6393457/144‑2 (MQ=255)
          tGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTgg                                                                                                                                             <  1:7905075/150‑1 (MQ=255)
           atcggtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTgt                                                                                                                                            <  2:9075707/146‑1 (MQ=255)
           atcggtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTgt                                                                                                                                            <  2:5170993/146‑1 (MQ=255)
            aTCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtg                                                                                                                                           >  1:4836630/1‑150 (MQ=255)
             tcgtcgATACGCTGCCCGGCGAATACCAGATGTTGCGCGCCGACCTCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtgt                                                                                                                                          <  2:2202642/150‑1 (MQ=255)
             tcgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtgt                                                                                                                                          <  1:100998/150‑1 (MQ=255)
             tcgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtgt                                                                                                                                          >  2:695668/1‑150 (MQ=255)
             tcgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtgt                                                                                                                                          <  1:2124810/150‑1 (MQ=255)
               gtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTgg                                                                                                                                        <  1:5706503/150‑1 (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)
                       gCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaa                                                                                                                                >  2:2950593/1‑150 (MQ=255)
                       gCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaa                                                                                                                                >  1:7057959/1‑150 (MQ=255)
                        cTGACCGGCGAATACACGATGTTGCGCGCCGACATCCTGCACGACGGGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaaa                                                                                                                               <  2:8139252/150‑1 (MQ=255)
                        cTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaaa                                                                                                                               <  2:4058501/150‑1 (MQ=255)
                        cTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaaa                                                                                                                               <  1:2363578/150‑1 (MQ=255)
                         tGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTGGGATGACAAAg                                                                                                                              >  1:265239/1‑150 (MQ=255)
                         tGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAg                                                                                                                              <  1:8636188/150‑1 (MQ=255)
                              ggCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAg                                                                                                                         <  1:1232349/150‑1 (MQ=255)
                                  aaTACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtga                                                                                                                     <  2:3275570/150‑1 (MQ=255)
                                   aTACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtgat                                                                                                                    <  2:5821820/150‑1 (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)
                                       aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                                                                                                >  1:1893113/1‑150 (MQ=255)
                                       aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                                                                                                >  1:3114530/1‑150 (MQ=255)
                                       aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                                                                                                >  1:2546466/1‑150 (MQ=255)
                                       aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                                                                                                >  2:846161/1‑150 (MQ=255)
                                       aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCAGATGAcc                                                                                                                >  1:7483725/1‑150 (MQ=255)
                                        aGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCa                                                                                                               >  2:1138337/1‑150 (MQ=255)
                                        aGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCa                                                                                                               >  1:2242984/1‑150 (MQ=255)
                                        aGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCa                                                                                                               >  1:7629960/1‑150 (MQ=255)
                                           tGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcg                                                                                                            >  2:8552975/1‑150 (MQ=255)
                                               gcgcgcCGACATCCTGCACTACGTGGGCGATTCGCTTAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggccc                                                                                                        <  1:3494526/150‑1 (MQ=255)
                                               gcgcgcCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggccc                                                                                                        <  2:7405986/150‑1 (MQ=255)
                                                 gcgcCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggcccgg                                                                                                      >  2:3436673/1‑150 (MQ=255)
                                                 gcgcCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggcccgg                                                                                                      >  2:2315653/1‑150 (MQ=255)
                                                 gcgcCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggcccgg                                                                                                      >  2:3535228/1‑150 (MQ=255)
                                                 gcgcCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggcccgg                                                                                                      >  2:3607225/1‑150 (MQ=255)
                                                             tGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGa                                                                                          >  1:7694251/1‑150 (MQ=255)
                                                                  gacgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggcccggc                                                                                                     >  1:5499910/1‑134 (MQ=255)
                                                                  gacgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCg                                                                                     >  1:4721602/1‑150 (MQ=255)
                                                                  gacgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCg                                                                                     <  2:470207/150‑1 (MQ=255)
                                                                    cgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGc                                                                                   >  2:4114283/1‑150 (MQ=255)
                                                                      tGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCg                                                                                 >  2:6812771/1‑150 (MQ=255)
                                                                      tGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCg                                                                                 >  2:5328363/1‑150 (MQ=255)
                                                                             ttCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCggccgggccgtt                                                                         <  1:5200969/150‑3 (MQ=255)
                                                                                gCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACAt                                                                       <  1:8450239/150‑1 (MQ=255)
                                                                                gCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACAt                                                                       >  2:8394169/1‑150 (MQ=255)
                                                                                     aCCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACa                                                                        <  2:8654084/144‑1 (MQ=255)
                                                                                      cccTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACt                                                                 <  1:4763423/150‑1 (MQ=255)
                                                                                       ccTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTg                                                                >  2:4519556/1‑150 (MQ=255)
                                                                                       ccTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTg                                                                <  2:2887433/150‑1 (MQ=255)
                                                                                         tGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGa                                                              >  1:3634287/1‑150 (MQ=255)
                                                                                         tGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGa                                                              >  2:2548212/1‑150 (MQ=255)
                                                                                          gCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGAc                                                             <  1:8347698/150‑1 (MQ=255)
                                                                                            gATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACtt                                                           >  1:6634211/1‑150 (MQ=255)
                                                                                            gATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACtt                                                           >  2:4005316/1‑150 (MQ=255)
                                                                                              tCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGc                                                         <  1:872596/150‑1 (MQ=255)
                                                                                                gACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCgg                                                       <  1:4339687/150‑1 (MQ=255)
                                                                                                    tCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGa                                                   >  2:3486852/1‑150 (MQ=255)
                                                                                                      ggCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAg                                                 >  1:2923423/1‑150 (MQ=255)
                                                                                                         cAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCtg                                              <  2:4745954/150‑1 (MQ=255)
                                                                                                          aGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCtgg                                             <  1:6422770/150‑1 (MQ=255)
                                                                                                           ggTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCtggt                                            <  1:9049233/150‑1 (MQ=255)
                                                                                                            gTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCACCGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCtggtg                                           <  1:837909/150‑1 (MQ=255)
                                                                                                            gTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCtggtg                                           >  2:5449347/1‑150 (MQ=255)
                                                                                                            gTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCtggtg                                           <  2:4213887/150‑1 (MQ=255)
                                                                                                            gTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCtggtg                                           <  2:3194184/150‑1 (MQ=255)
                                                                                                            gTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCtggtg                                           <  2:2460374/150‑1 (MQ=255)
                                                                                                             tGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCtggtgg                                          <  1:3193469/150‑1 (MQ=255)
                                                                                                                  ggtggtTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAAc                                     <  1:3033908/150‑1 (MQ=255)
                                                                                                                       tttCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCaa                                >  1:5067349/1‑150 (MQ=255)
                                                                                                                         tCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGa                              >  1:8058089/1‑150 (MQ=255)
                                                                                                                          cATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGaa                             >  1:4498511/1‑150 (MQ=255)
                                                                                                                           aTCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAAc                            <  2:1073023/150‑1 (MQ=255)
                                                                                                                             ccAGGGCATGGGCTGGCTGACCACCGAATAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAAccc                          <  2:266837/150‑1 (MQ=255)
                                                                                                                               aGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAAccc                          <  2:8906354/148‑1 (MQ=255)
                                                                                                                               aGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAAccc                          >  1:8906354/1‑148 (MQ=255)
                                                                                                                               aGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCgg                        >  1:8937389/1‑150 (MQ=255)
                                                                                                                                   cATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACAt                                                                       <  1:7155676/99‑1 (MQ=255)
                                                                                                                                   cATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGa                    >  2:2078515/1‑150 (MQ=255)
                                                                                                                                        gctggctgACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGt                                                                                <  1:6802733/85‑1 (MQ=255)
                                                                                                                                        gctggctgACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCg               <  1:2993028/150‑1 (MQ=255)
                                                                                                                                          tggctgACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCgtg             <  1:1292049/150‑1 (MQ=255)
                                                                                                                                           ggctgACCACCGAAGAGCTGGTGTGGAACGCCTAAGGCAAGCTGAGGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCgtgt            <  2:5841226/150‑1 (MQ=255)
                                                                                                                                                  caccGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCACTc     >  1:1611919/1‑150 (MQ=255)
                                                                                                                                                  caccGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCACTc     >  2:1738100/1‑150 (MQ=255)
                                                                                                                                                  caccGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGAGCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCACTc     <  1:7264681/150‑1 (MQ=255)
                                                                                                                                                   accGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCACTcc    >  2:6323305/1‑150 (MQ=255)
                                                                                                                                                     cGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCCGACATGCCACTGGACTTGCGGGTGAAGCTGGTGGAGAACCGCAAGAACCCGGAAGACACCGTGTTCCACTCCaa  >  2:7189445/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: